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_ from Donald W. Douglas, Jr. 


In your field, as most of you well know, it’s easy 
to be complicated ... it’s hard to be simple. At 
Douglas, I’m happy to say, we do things the 
“hard” way. This matter of simplicity is vitally 
important. We work intensive hours, days and 
months to achieve it. 

Why this extra effort? Well, simple things work 
easier, last longer, are more easily maintained 
and are lots more reliable. We are rewarded for 
our greater engineering effort with a product 
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President, Douglas Aircraft Company 


that performs better for our various customers. 

We know that good engineers, working in an > 
atmosphere which stimulates them to do their 
best, have been largely responsible for our suc- 
cess. If you enjoy solving challenging problems in | 
the simplest manner, we’d like to talk with you. 
about joining us. 

Please write to Mr. C. C. LaVene 
Douglas Aircraft Company, Box J-620 
Santa Monica, California 
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Scope of JET PROPULSION 


This Journal is a publication of the Ameri. 
can Rocket Society devoted to the advange. 
ment of the field of jet propulsion throug) 
the dissemination of original papers disclosing 
new knowledge or new developments, gs 
used herein, the term ‘“‘jet propulsion” 
embraces all engines that develop thrust by 
rearward discharge of a jet through a nozzle 
or duct; and thus it includes air-consuming 
engines and underwater systems as well as 
rockets. JET PROPULSION is open to con- 
tributions dealing not only with propulsion 
but with other aspects of jet-propelled flight 
such as flight mechanics, guidance, telemeter. 
ing, and research instrumentation. [p-. 
creasing emphasis will be given to the scientific 
problems of extraterrestrial flight. 


Information for Authors 


Manuscripts must be as brief as tie proper 
presentation of the ideas will alluw. Ex. 
clusion of dispensable material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references 
1 page of abstract, and 12 illustrations. 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
in exceptional cases. 

hort manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes or Technical Comments, 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are usually published within two 
months of the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to tile indus- 
try. A manuscript that does not qualify for 
publication according to the above-stated 
requirements as to subject, scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as & 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 

word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows: for journal articles: authors first 
then title, journal, volume, year, page numbers; 
for books: authors first, then title, publisher 
city, edition, and page or chapter numbers 
Line drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. F 

Manuscripts must be accompanied by 
written assurance as to security clearance 1D 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. at 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
Professor of Aeronautical Engineering, Prince- 
ton University, Princeton, N. J. Preprints 
of papers presented at ARS national meetings 
are automatically considered for publication. 
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The author was born in 1922 in Italy. Following a clas- 
ical education, he attended the University of Rome and was 
warded a doctorate in Aeronautical Engineering in 1946. Dr. 
Miele has held positions in engineering, teaching, research 
ad consulting in Italy, Argentina and the United States. Heis 
wwa professor at Purdue University, where he has taught for 
fe last three years. He has acted as a consultant to Doug- 
js Aircraft Company and Allison Division of General Motors 
(poration; and has directed research projects sponsored by 
he National Advisory Committee for Aeronautics, Office of 
iientific Research of the United States Air Force and Army 
jullistic Missile Agency. Dr. Miele’s research has been 
minly concerned with the mechanics of terrestrial and extra- 
errestrial flight and with the application of the Calculus of Vari- 
tions to the general theory of aircraft, missile and satellite per- 
mance. He is the author of approximately forty scientific 


pers, published in American and European journals. 


Introduction 
HE mechanics of flight of an aircraft is an ‘important 
branch of the aeronautical sciences. The scope of this 
ubject—a final synthesis of classical mechanics, theoretical 
ud experimental aerodynamics and theory of jet engines— 
sthe determination of the performance of a flying machine. 
lsparticular objective is to provide a solution for the follow- 
ig fundamental and complementary problems: 
| To define the operational utility of an aircraft or missile 
thich has already been designed. 


ussile for a given mission to be accomplished. 
In the years preceding World War II it was a common — 
lief among engineers that the mechanics of flight had — 
miched a conclusive and rather stagnant stage of develop- 
ett. The progress achieved in the last fifteen years has 
sproved this belief. In fact, the advent of jet engines as 


rraft propulsion systems and the parallel increase in flight 
locities have generated a wealth of important and unsolved 
~wblems of applied mathematics. Some of these problems, 
lich cannot be handled by conventional methods of per- 
imance analysis, are mainly responsible for the engineering 
wival of the branch of mathematics known as the Calculus 
i Variations. 


Object of the Present Paper 


This article offers a critical sketch of some recent advances 
ithe mechanics of terrestrial flight. In spite of the con- 
werable amount of bibliographical information included, 
welaim for completeness is made. Consequently, the author 
ologizes to those original writers who, because of oversight 
‘space considerations, have been excluded from the present 
arvey. 

In this paper the aircraft is conceived as a particle. The 
“dynamic lag is disregarded and attention is focused on 


— 
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Some Recent Advances in the Mechanics of 
_ 


two categories of questions associated with motion in a 
vertical plane:? Problems of steady or quasi-steady flight, 


and problems of nonsteady flight. 5 
Problems of Steady or Quasi-Steady Flight _ wal 


The quasi-steady approach is useful in the analysis of those 
conditions of flight where the acceleration can be regarded as 
negligible. Initially developed for conventional aircraft, 
the quasi-steady approach has also been of interest for air- 
craft propelled by air-breathing jet engines. As an example, 
the cruising flight of a turbojet powered commercial aircraft 
can be well analyzed by means of this approach. The latter 
has the merit of leading to simple analytic relationships of 
great usefulness for design purposes. 

Inherent in the quasi-steady approach is the fact that each 
point of the flight path of an aircraft can be investigated inde- 
pendently of the preceding one and of the following one. 
Thus, the so-called point performance can be studied within 
the framework of elementary algebra. Furthermore, as- 
sociated with the optimum flight conditions are analytical 
problems amenable to treatment by the Ordinary Theory of 


Maxima and Minima. 
Let 
T = thrust 
D = drag 
‘@ L = lift 
= weight 


= inclination of the vector velocity with respect to a hori- 
zontal plane 


a = inclination of the vector thrust with respect to the vec- 
tor velocity 
= altitude 
V = flight velocity 


Consider a drag function of the form D = D(h, V, L) and a 
thrust function 7 = T(h, V, a@) where a is a variable* con- 
trolling the engine performance and having a function analo- 
gous to that of the accelerator pedal in a typical automobile. 
With the above considerations in mind, and assuming 
that the inertia terms are negligible, the equations of motion 
in a vertical plane are written as 


T(h, V, a) cos w — Dih, V, L) — W sin 6 = 
T(h, V, a) sinw + L — Weos 6 = 0 


For a given weight W, Equations [1] represent a set of two 

2 Due to space limitations, the problems of motion in a hori- 
zontal plane are not reviewed here. The interested reader is 
referred to References (1, 19, 20, 21, 69, 71). 

3 Within the context of this paper the a-variable is termed en- 
gine control parameter or thrust control parameter or power 
setting. Nevertheless, its actual physical meaning is not speci- 
fied, both because of the great variety of existing types of engines 
and for the purpose of developing a general theory. As an ex- 
ample, however, a can be identified with the number of revolu- 
tions of the turbine-compressor group for a turbojet engine hav- 
ing fixed geometry, or with the position of the fuel control lever, 
ete. 


Note: 
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This is the first of two articles on the mechanics of 5 ong” 
‘restrial flight and is entitled “Recent Advances in Astrodynamics’ 


The other to follow next month is on the subject of extra- 


’ by S. Herrick and R.M.L. Baker Jr. 
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algebraic relations in the six variables h, V, L, 0, a and w. 
It follows that four degrees of freedom are left, thus imply- 
ing that a quadruple infinity of solutions is mathematically 


possible. 
Additional Constraints 4 


In many engineering applications it is of paramount interest 
to study particular solutions of the algebraic system formed by 
Equations [1], more specifically, those solutions which 
simultaneously satisfy a set of p < 4 additional constraints.* 
The latter are symbolically indicated as 


V, L, 8, a, w) = 0 (i= 


Clearly, the number (n) of degrees of freedom associated 
with the system of Equations [1 and 2] is now 


As a consequence, for p < 3, n > 1, so that some optimum 
requirement can be imposed on the flight performance, specifi- 
cally, the one of extremizing (i.e., maximizing or minimizing) 
an arbitrarily specified function of the form 


In practical cases the W-function can be, for instance, the 
range* per unit fuel consumed (V/cT7'), the endurance per unit 
fuel consumed (1/c7'), the altitude (h), the flight velocity (V), 
the path inclination (6), the rate of climb (V sin 6), ete. 

, Within the frame of expressions [1 to 4] many of the recent 
advances in quasi-steady flight mechanics of jet-propelled 
aircraft can be interpreted. First, the comprehensive papers 
by Lippisch (1,2)® deserve particular mention. They were 
varried out as part of the pioneering development of jet- 
propelled vehicles in Germany during World War II. 

In the years following World War II, the performance 
problem attracted considerable interest in several countries. 
Under the assumption that the drag polar is parabolic with 
constant coefficients, analytical results were derived in (3 to 
10). 

The important range problem for a turbojet aircraft was 
the object of the investigations (11 to 13). In particular, the 
economy in fuel, which may result by flying at variable alti- 
tude (cruise-climb technique), was pointed out in (11,12). 

The case of a drag polar embodying—in addition to a con- 
stant term and a term quadratic in the lift coefficient—a 
term linear in the lift coefficient was considered by Ludwig 
(14), who extended to a jet-propelled aircraft the analysis 
previously carried out by Wendt (15) for a piston-engined air- 
craft. Finally, the computation of the flight performance 
with compressibility effects accounted for was stressed in (3, 
16 to 18). 


Optimum Conditions for a Given Configuration and Basic 


Determinantal Equation 


Consider Equations [1 and 2], and assume p = 8, i.e., 
that three additional constraints ®; = 0, , = 0 and 4; = 0 
are imposed. Assume also that the only remaining degree 
of freedom is to be saturated in such a way that the W-func- 
tion (Eq. [4]) is stationary. By the method of the Lagrange 
multipliers—and within the framework of the Ordinary 
Theory of Maxima and Minima—it can be shown (19) that 


4 Within the context of the present paper the constraints [1] 
are a fundamental; the constraints [2] are termed addi- 
tional. 

_ * The symbol c = c(h, V, «) denotes specific fuel consumption, 
i.e., weight of fuel consumed per unit time and unit thrust. 
® Numbers in parentheses indicate References at end of paper. 


the following optimizing condition is to be satisfied 
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The above determinantal equation is fundameital for 
flight mechanics applications. Actually, all the optim:m flight 
conditions of articles (1 to 18) are only particular cases oj 
Equation [5]. The latter, therefore, unifies the best ;erforn- 
ances in a vertical plane into a single generalized st:tement 
Thus, apparently diversified problems, such as maximun 
range, maximum endurance, ceiling, maximum spee«|, steep- 
est ascent, fastest ascent, flattest descent, etc., are all coy- 
ered’ by the same determinantal equation. 


Optimum Aerodynamic Configurations 


The method outlined in (19) not only yields the best flight 
performance for a given aircraft configuration, but also cai 
be employed to solve design questions. 

Assume for instance that the configuration of an aircraft’ 
depends on m parameters 7, (k = 1...m), so that the drag 
function becomes D = D(h, V, L, m, . . ., tm) and the weight 
W = W(m,...,7m). Assume also that the above m param- 
eters must be determined so as to extremize the functio: 
[4] subjeet to the equations of motion and to the constraints 
[2], where p = 3. The number of degrees of freedom becomes 
n = m + 1, so that—after eliminating the Lagrange multi 
pliers—(m + 1) optimum relationships must be writter 
Under the hypothesis that the W-function and the ,-con- 
straints are independent of both the lift and the z,-parameters, 
it can be shown (19) that the optimizing condition is repre- 
sented by the combination of Equation [5] and the following 
m additional relations 

ow 


D D 
(«in 6+ cos 0) =0 (k= 


Om; 

In spite of the generality of these results, a word of caution is 
necessary. Unavoidably, an aircraft design is always ‘ 
compromise between many contrasting requirements. As 3 
consequence, the previously indicated optimum aerodynamic 
configurations must be taken cum grano salis. 


Problems of Nonsteady Flight 


The nonsteady approach to the mechanics of flight is 0! 
interest (for instance) in the analysis of the flight paths 0! 
rocket-powered aircraft, missiles, satellite carriers, skip ant 
glide bombers and in the study of the transient behavior’ 0 
aircraft propelled by air-breathing engines. For flight in 4 
vertical plane over a flat Earth"? the following set of differen 


7 The only difference between one problem and another lies ! 
the analytical nature of the function W to be extremized and 0! 
the additional constraints #; = 0 (7 = 1, 2, 3) to be satisfied. 

5 In practical cases the +; parameters can be, for instance, thé 
wing surface, the aspect ratio, the apex angle, ete. 

* Typical, among these problems, is the investigation of the 
brachistocronic climbing technique for a turbojet aircraft. ' 

1” Problems of motion over a spherical, rotating Earth are di 
cussed at the end of this paper. 
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jequations is to be considered. ; : The Problems of Bolza, Mayer and Lagrange 


With regard to the optimization of the flight trajectories 
_ the substantial difference between the quasi-steady approach 
_ and the nonsteady approach is the following: Because of the 


X — Vecos 6 = 0 
h-Vsin@ =0 
D(h, V, L) — T(h, V, a) 


+ g sin @ + — “ =e inertia terms, the minimal problems can no longer be handled 
cos6 L + Th, V, a) sinw from a local point of view; considerations of an integral na- 
@+4a ert _ = @ ture are in order, in the sense that the flight path is to be 
b treated as a whole. Thus, the nucleus of the problems of 
‘ (7] performance optimization shifts from the domain of the 
Ordinary Theory of Maxima and Minima into the realm of the 
Calculus of Variations. This fact, in turn, signifies a con- 
; : siderable increase in analytical difficulties, insofar as the 
\ = horizontal distance : optimum paths are no longer governed by equations in finite 
acceleration of gravity ; terms, but rather controlled by nonlinear systems of dif- 
= Seen ae of the vehicle ferential equations which are integrable, in general, only by 
Notice that, if the additional constraints are such that p < 
ae Notic that Equations [7] form a system of five differential 2, the number of degrees of freedom isn > 1. As a conse- 
verform.(q uations in the eight unknown functions X(t), A(t), V(t), quence, some ae requirement can be imposed on the 
tement L(t), a(t) and w(t). Thus, n= 3 degrees of free- flight trajectory. In this connection, the following problem 
adie jon are left, as is logical in view of the possibility of con- (of Bolza type) is formulated (22): “In the class of functions 
L, steep. tolling the time history of: The elevator position, the thrust x (t), h(t), V(t), A(t), m(t), Lit), a(t) and w(t) satisfying Equa- 
pid nagnitude, and the thrust direction.!* tions [7 and 8] and certain prescribed end conditions, find 
that special set such that the functional form? 
W = (G(X, h, V, 0, m, 4+ VC my 
7 ti 
st flight 
also cat is extremized.”’ 

Additional Constraints It should be remarked that, if the G-function is iden- 
aircraft’ tically zero, the Bolza problem reduces to the problem of 
he drag In many ehginccring applications It Is of interest to study Lagrange. If, on the other hand, the F-function is identically 
» weight particular solutions of the set of differential Equations [7]; zero, the Bolza problem reduces to the problem of Mayer. 
param: fl '®s those solutions which simultaneously satisfy a set of It is also stressed that, while the problem of Bolza appears to 
junction ? <3 additional constraints be—at first glance—more general than the other two, the 
straints 0(X, hb, V, 6, m; X,h, V, 0, m, L, a, w, t) = 0 three problems of Mayer, Lagrange and Bolza have, actually, 
becomes (@=1,...,p)..{8] the same degree of generality. As a matter of fact, by simple 
e multi transformations (22), one problem can be taken into another 
written. Their effect is to reduce the number of degrees of freedom to and vice versa. A dials 1 

-COll- 
General Considerations 


s repre-§ Within the framework of expressions [7 to 9] many of the 
lowing rent contributions to nonsteady flight mechanics can be 
interpreted. Quite appropriately, these contributions can be 


With the theory of the first variation it can be shown that 
the special trajectory which solves the proposed Bolza prob- 


livided into papers dealing with optimum flight trajectories, lem must satisfy, as a necessary condition, the set of Euler- 
m)..{6)§ ad papers dealing with arbitrary (i.e., non-optimum) flight Lagrange equations (22), omitted here for brevity. The for- 
trajectories. mulation of the latter becomes particularly simple if use is 
anil Because of space limitations, only a survey of papers of made of a special set of auxiliary variables, the Lagrange 
the first type is presented here. Furthermore, they are gener- multipliers. The boundary conditions (partly of the fixed 


4s Uy more significant from a mathematical standpoint. Con- 


iibutions of the second type are therefore excluded, even 


end points type, partly of the natural type) must be con- 


wees sistent with the general transversality condition of the Calculus 


though there is no question that several of them are cer- of Variations (22). ; : 
tainly interesting from an engineering point of view. It is to be noted that the analytical nature of the set of 
: = Euler equations and constraining equations is such that the 
"It should be emphasized that the phrase “degree of freedom’’ attainment of solutions in : closed form is the reaper 
_ [i sbeing used in this section with a significance different from the rather than the rule. An important additional difficulty 
ht is o! ile used in the preceding section on Problems of Steady or Quasi- stems from the circumstance that the variational problems of 
yaths ol re Flight. The symbol n of the preceding section denotes the mechanics of flight are boundary value problems, i.e., 
<ip and ner of coordinates to be given in order to supply one par- problems with conditions prescribed in part at the initial 
of ar solution of the set of algebraic Equations [1 and 2]. On 
v10r' "Bf te other hand, the symbol n of this section indicates the number point and in part at the final point. As a consequence, the 
ht i “fiffunctions to be specified in order to supply—with the appro- use of trial and error techniques is, except for particular cases, 
|ifferen-§ tate boundary conditions—one particular solution of the set of an unavoidable necessity. 
iiferential Equations [7 and 8]. 
._.f# "In severa! papers of the variational literature the angle of 
7 — y tack (€) has been used in place of the lift (Z) as a dependent 13 Tt should be emphasized that, for n = 0, it is still possible to 
ri “Bitiable. Analytically, Equations [7] are modified insofar as the optimize the boundary conditions. Problems of this type have 
ed. h irag becomes D = D(h, V, e) and the lift L = L(h, V, «). The been treated in (58). 
ice, WH minimal problem, therefore, is no longer studied in the (X, h, V, 14 The subscript 7 denotes initial point; the subscript f refers 
+ the 4m, L, a, w, t) space but in the (X, h, V, 6, m, €, a, w, £) space. to final point. The symbols G and F denote arbitrarily specified 
| Ol WE The latter coordinate system offers—perhaps—a more imme- functions of the coordinates of the vehicle and their derivatives. 
ode te understanding of the analytical nature of the problem in the 16 Typical G-functions, having engineering interest, are: The 
ire GF se where particular hypotheses are introduced in the equations range (G = X), the time (@ = 2), the flight velocity (G = V), the 
(motion (see following sections). mass (G = m), ete. 
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Review of the Variational Literature 


Generally speaking, the analytical difficulties of a Bolza 
problem increase as the number of degrees of freedom in- 
creases. Thus, the simplest problems are those in which two 
degrees of freedom have been conveniently frozen, so that 
attention is focused on the optimization of the remaining one. 
Typical problems arising from the optimization of only the 
lift distribution (or the angle of attack) or only the thrust 
modulus or only the thrust direction are now reviewed. 
More general problems are also considered.'® 


Programming. of the Lift Distribution or the Angle 
of Attack 


. Consider the case p = 2, and assume that the two con- 
straints @, = 0 and & = 0 are such that they prescribe a 
priori both the mode of operation of the engine and the thrust 
a direction (relative to the aircraft or to the vector velocity). 
One degree of freedom is consequently left, namely the de- 
gree of freedom associated with the programming of either 


the lift distribution or the angle of attack. Some typical 
problems are now presented (23 to 39). 


The initial investigations of the optimum climbing pro- 
gram for turbojet aircraft were based on a somewhat simpli- 
fied form of Equations [7]. The assumption was made that 
the mass of the aircraft is ideally a constant, and the centri- 

- petal component of the acceleration is disregarded, with only 

the tangential component being accounted for. Thus, ap- 
proximate solutions were detected by the writer (23,24) whose 
analytical method, based on Green’s theorem, avoids the 
use of indirect variational procedures. The Calculus of 
Variations was also avoided in a paper by Lush (25), with an 
intuitive graphical method based on the concept of energy- 
height. 

On the same problem, an attack with the indirect methods 

_ of the Calculus of Variations was attempted in (26 to 28) 

leading, however, to incomplete results. In fact, the true 
Eulerian solution is discontinuous and generally formed of 
three different branches (dive, zoom, variable path inclina- 
tion sub-are); nevertheless, only one of the three sub-arcs 
_ forming the extremal are was detected in (26 to 28), namely, 
the sub-are flown with variable path inclination. Complete 
variational solutions, on the other hand, were found in 

(29,30). 

Recently, the results indicated in (23,24) were rederived by 
Cartaino and Dreyfus (31) along the lines cf the promising 
theory of dynamic programming, as developed by Bellman 

_ (70). Furthermore, solutions in a closed form for the alti- 

_ tude-Mach number relationship were indicated in (32), thus 

bypassing the so-called energy-height method commonly 
used by aircraft manufacturers. 


Simplified Climbing Program for Rocket-Powered Aircraft 


For a rocket-powered aircraft, the timewise variation of 
_ Mass is important and must be accounted for in the study of 
the climbing performance. Historically, it is worth noticing 


16 Tt should be emphasized that the really important charac- 
_ teristic of a variational problem is the number of degrees of free- 
dom associated with the set of fundamental constraints and addi- 
tional constraints being considered. Any other form of classifica- 
tion is—to a certain extent—arbitrary. The terminology “‘ift 
: programming problems” or “thrust programming problems’’ or 
“problems where the thrust direction is programmed” has only 
_ been used with the purpose of rendering a more vivid exposition 
of these problems. The reader, therefore, should accept the pres- 
. ent terminology with a note of caution. 
This problem, as well as the one of the following section, is 
better understood if the angle of attack is used in place of the lift 
as a dependent variable. 
_ 18 With regard to either time or fuel consumed. 
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that the theoretical development of new climbing techniques 
for rocket-powered aircraft preceded that for turbojet pow. 
ered aircraft. 

The two papers (1,33) are of particular interest. Even 
though they left the bulk of variational questions associated 
with the climbing flight substantially unsolved, nevertheless 
they shed considerable light on a new class of problems of the 
mechanics of flight. 

It appears that in the paper by Kaiser (33) the concept of 
energy-height was first employed for flight mechanics «pplica- 
tion. On the other hand, the accelerated climbing flight was 
investigated in (1) under the assumption that the total drag 
coefficient is constant. As in the case of (26 to 2s), the 
mathematical techniques employed in (1) lead to an incom- 
plete solution, insofar as only one of the three sub-arcs form- 
ing the extremal arc associated with the problem under con- 
sideration was detected. 

Recently, the brachistocronic problem with fixed end points 
was attacked in (34) with Green’s theorem and in (35) with 
the indirect methods of the Calculus of Variations, uncer the 
assumption that the induced drag is negligible with respect 
to the zero-lift drag. In particular, analytical solutious were 
detected in (35) for the altitude-Mach number relationship, 
for both tropospheric and stratospheric flight. 


More General Investigations 


The investigations of the preceding sections were carried 
out under particular hypotheses, of which the essential an- 
alytical objective was to simplify the equation of motion on 
the normal to the flight path. By lifting the above limita- 
tions, a more general category of variational problems is 
originated. These problems attracted the attention of 
Cicala and the writer, who stressed the advantages of the 
Mayer formulation in an initial note (36) dealing with mini- 
mum time flight paths. 

On the other hand, the Lagrange formulation was used by 
Behrbohm (37) who investigated—among other problems— 
trajectories of minimum time and of minimum fuel con- 
sumption, by referring to a parabolic polar with coefficients 
depending upon the Mach number. 

In a recent Readers’ Forum note (38) the minimum time 
problem was once more attacked essentially along the lines of 
(36), but using the altitude in place of the time as the inde- 
pendent variable. 

Finally, a recent paper by Fraeijs de Veubeke (39) is to be 
mentioned. It refers particularly to the problem of maxi- 
mum range and formulates the equations of motion in para- 
metric form. 


Thrust Programming Problems 


Consider again the case p = 2, and assume that the con- 
straints 6; = ®, = 0 prescribe a priori the direction of the 
thrust and the distribution of lift (or some equivalent geo- 
metric condition). One degree of freedom is left, the one 
associated with the programming of the modulus of thrust. 
Papers in this area of problems (40 to 53) are now reviewed. 

Concerning the level flight of a rocket-powered aircraft, 
the early analysis by Hibbs (40) is mentioned in connection 
with the problem of maximizing the range, the flying time 
being free of choice. An interesting characteristic of this 
problem is that—in addition to the constraints [7 and 8]—it 
involves a side condition expressed in the form of an inequality 


specifying that the engine mass flow 8 is bounded between al 
upper value and a lower value. This circumstance, in tur, 
leads to a discontinuous extremal arc, composed of sub-ares 
of three different kinds: Sub-ares flown with maximum thrust, 
sub-ares flown by coasting and sub-ares flown with variable 
thrust. Methods for handling inequality [11] were de 
veloped in (41) by the writer, who made use of Green’s the- 
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orem and by Cicala, who introduced the concept of inder 
ylue in combination with the Weierstrass device. More 
general Mayer problems of the form G = G(X, V, m, t) were 
investigated in (42). 

With regard to vertical flight, it is of historical interest to 
notice that Goddard (438) clearly recognized that the Calculus 
of Variations is an important tool in the performance analy- 
sis of rockets; even though—following an oversimplified pro- 
cedure—he dealt with questions of the Ordinary Theory of 
Maxima and Minima. Following a brief paper by Hamel 
(44), the variational problems of minimum propellant con- 
sumption attracted the attention of Tsien and Evans (45) 
and, subsequently, of Ward (46), Leitmann (47,48) and Faulk- 
ner (49). In addition, the brachistocronic problem, as well 
as Mayer problems of the form G = G(h, V, m, t) were stud- 
ied in (50,51). 

Arbitrarily inclined rectilinear paths were investigated in 
(52), while nonlifting paths flown with minimum propellant 
consumption were analyzed by Bryson and Ross (53). 


Programming of the Thrust Direction 


If the two constraints ©; = ®, = 0 prescribe a priori both 
the distribution of lift and the modulus of the thrust, the 
trajectory can still be optimized with regard to the thrust 
direction (54 to 56). 

For the special condition of vacuum flight (L = D = 0), 
simple and intriguing results are obtained. It . been shown 
by Lawden (54) that, for a Mayer problem of the form G = 
G(X, h, V, 0, m, t), the absolute inclination in space of the vec- 
tor thrust is a bilinear function of the time. 

On the other hand, for the particular case of Mayer prob- 
lems of the form G = G(h, V, 8, m, t) not involving range, the 
absolute inclination of the vector thrust becomes a linear 
function of the time, as shown by Fried (55). Finally, for 
the other particular case of Mayer problems of the form G = 
G(X,h,m,t) where the vector velocity is free at the final point, 
the optimum inclination in space of the vector thrust is a 
constant, a result proved by Fried and Richardson (56). 


Simultaneous Optimization of Several Degrees of Freedom 


The investigations of the preceding sections were con- 
cerned with problems involving only one degree of freedom, 
ie., the one associated with either the lift distribution or the 
thrust modulus or the thrust direction. In the present para- 
graph questions of a more general nature, involving either 
two or three degrees of freedom, are analyzed (57 to 60). 

First, the contributions due to Cicala (57,58) must be men- 
tioned in connection with the programming of both the lift 
distribution and the modulus of the thrust. The latter 
problem was also studied in an interesting report by Break- 
well (59). Finally, a very general theory by the writer (60) 
is mentioned in connection with the simultaneous optimiza- 
tion of lift distribution, thrust modulus and thrust direction. 


Legendre-Clebsch, Weierstrass and Jacobi-Mayer 
Conditions 


As is known from the theory of the first variation, the opti- 
mum paths are solutions of the Euler-Lagrange differential 
equations. The fact that a trajectory satisfies the Euler 
equations, however, only guarantees its stationary character. 
To determine whether the functional VY assumes a maximum 
or minimum value along the Eulerian path, a further investi- 
gation is in order. In particular, the so-called Legendre- 
Clebsch, Weierstrass and Jacobi-Mayer conditions must be 
analyzed. A general discussion of these conditions can be 
found in (57,58). Particular aspects of them are dealt with 
in (39,55,61). 

It is of interest to note that in many problems of the me- 
chanics of flight the conditions due to Legendre-Clebsch and 
Weierstrass fail to supply any information on the minimal or 
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maximal character of the solution along a major portion of — 
the extremal arc. This important difficulty can be overcome 
with some of the techniques developed in (58). Another de-— 
vice is based on Green’s theorem, which has a particularly | 
useful application in the study of linear variational problems _ 
of either the simplest type or the isoperimetric type (24, 42, — 
62). 

J 


Effect of Inequalities on the Solutions 


There are innumerable problems of the mechanics of flight — 
where not only the constraints [7 and 8] must be satisfied, — 
but also a number of additional constraints expressed in the | 
form of inequalities. 

As an example, the mass flow of a rocket engine must satisfy — 
the inequality [11] at all points of the flight path. Other — 
possible inequalities may involve the longitudinal accelera- — 
tion, the normal acceleration, the lift, the flight altitude, 
the maneuverability, etc. In this connection, the following 
statement is pertinent: The effect of each inequality is to— 
introduce one additional (potential) branch into the solution. 
Discontinuous solutions are generally found and the Erdmann- 
Weierstrass corner conditions (22) must be applied in order — 
to join the different sub-arcs forming one extremal arc. 

As an additional example, consider a high speed glider, and 
assume that the lift Z is to be consistent with the constraint 


at all points of the flight path, W being the weight and K a 
nondimensional constant. By means of analytical methods 
similar to those developed in (42,51,53) it can be shown that, 
because of the inequality [12], the additional sub-arc 


L—KW= [13] 


may appear in the composition of the extremal are. In this 
connection, methods for treating problems involving dis- 
continuous solutions can be found in several of the papers of 
Cicala and the writer. 


On the Relationship Between Calculus of Variations and 
Ordinary Theory of Maxima and Minima 


An interesting link can be established between Calculus of 
Variations and Ordinary Theory of Maxima and Minima by 
means of the following procedure. 

Let the acceleration terms appearing in Equations [7] de- 
crease, tending to zero in the limit. The associated set of 
Euler-Lagrange equations modifies and, in the limit, yields a 
new set of Euler equations. The latter set supplies optimiz- 
ing conditions identical with those given by the Ordinary 
Theory of Maxima and Minima for problems of quasi-steady 
flight. This important property of the Euler equations has 
been proved in (17,19,63) for particular problems, specifically, 
the maximum range and the maximum endurance of a turbo- 
jet powered aircraft. 


Problems of a Thermal Nature 


The variational techniques surveyed in this article can also 
be applied to problems of a thermal nature, for instance the 
problem of the trajectory which minimizes the increase in 
skin temperature for a rocket-powered missile (64). 

The analytical character of these problems is not different 
from those previously examined. Nevertheless, the set of 
constraining equations is to be completed by a further equa- 
tion, the one relative to the heat balance between missile 
skin and surrounding medium (64). 


Variational Problems with Consideration of the Variable 
Gravitational Field and Rotation of the Earth 


For problems where either great ranges or great velocities 
are involved, or for problems where a gret at precision is nee ded, 


| 
Rix 
| 
a 
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it becomes important to account for two factors: 

1 The gravitational field of the Earth is net uniform but 
central. 

2 Newton’s law is rigorously true only for a reference 
system rigidly associated with the so-called Fixed Stars. 
Thus, when the equations of motion are referred to a coor- 
dinate system rigidly associated with the Earth, the correc- 
tive terms due to the so-called transport acceleration and 
Coriolis acceleration must be considered. The latter is of more 
relevant importance than the former for ballistic missiles 
and satellite vehicles. 

When items 1 and 2 are accounted for, Equations [7] are 
modified and become more complicated. Nevertheless, the 
analytical nature of the minimal problems does not change; 
they still fall under the formulations of either Bolza, Mayer 
or Lagrange, even though the attainment of solutions in a 
closed form becomes more difficult. 

Problems involving either one or several degrees of free- 
dom (with sphericity and/or rotation of the Earth considered) 
have been investigated by Fried and Richardson (56), 

- Breakwell (59), Miele (60), Leitmann (65), Okhotsimskii 
and Eneev (66) and Lawden (67,68). The latter is also the 
author of many contributions to the mechanics of extra- 
terrestrial flight, excluded however from the present survey 
which is confined to problems of terrestrial flight. 


waa Summary and Conclusions 


A general survey is presented of the problems of the me- 

- chanics of terrestrial flight, with particular regard to the 
optimization of the trajectories of both quasi-steady and non- 
steady flight. 

The optimum problems of quasi-steady flight are amenable 
if as to treatment by the Ordinary Theory of Maxima and Min- 
ima. For the case where the aircraft configuration i is given, a 
= new and general determinantal equation is presented. The 

_ latter covers a wide variety of problems, such as maximum 

range, maximum endurance, ceiling, maximum speed, steep- 
est ascent, fastest ascent, flattest descent, etc., thus unifying 
all the mechanics of flight in a vertical plane into a single 
statement. A general solution is also presented for the 
- aerodynamic configuration which is the best for an arbitrary 
performance to be optimized. 

On the other hand, nonsteady flight mechanics leads to 
minimal problems falling in the domain of the problems 
_ of Bolza, Mayer and Lagrange of the Calculus of Variations. 
_ The recent variational literature is reviewed with particular 
— regard to the programming of the lift distribution, thrust 
-modulus and thrust direction. The following topics are also 
_ discussed: Legendre-Clebsch and Weierstrass conditions, 

effect of inequalities on the solutions, problems of a thermal 

nature, relationships between Calculus of Variations and 

Ordinary Theory of Maxima and Minima, and variational 

_ problems with consideration of the variable gravitational 
— and rotation of the Earth. 

It is clear from the present survey that much has been 
~ achieved in recent years in the field of terrestrial flight me- 
_ chanies. Many problems have been conquered. Neverthe- 
4 less, an even larger domain is still unexplored, both from a 
a theoretical standpoint and with regard to practical engineering 

applications. 
_ There is an immediate need for improved methods for inte- 
Mi bse the system of Euler equations and constraining equa- 
fp tions and for solving the associated boundary value problems. 
_ An extension of the available closed form solutions would be of 
great value with regard to engineering applications. In view 
of the rather weak character of the maxima and minima of the 
mechanics of flight, the finding of short cuts and simplifica- 
tions applicable to particular problems would also be valuable. 
At the present time, the work in the area of sufficient condi- 


tions for an extremum lags far behind the work accom- 
These sufficiency 


plished in obtaining necessary conditions. 
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conditions have given rise to questions, the answers to some 
of which are still incomplete or unknown, especially in eon. 
nection with discontinuous extremal solutions. 

In the era of supersonic interceptors, intercontinental mis. 
siles, satellites and interplanetary vehicles, the variational 
methods constitute a much needed and important step for. 
ward in advanced performance calculations. The vista is now 
beginning to open on this application. It is the opinion of the 
writer that, as the industry progresses toward faster and 
faster vehicles, the Calculus of Variations will become the 
standard, rather than the specialized, tool of optimum per. 
formance analysis of aircraft and missiles. 
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‘ HIS i is a preliminary report of results obtained concerning 
; ; radiation intensities measured with a single geiger tube 
earried by the artificial earth satellites 1958 a and 1958 y.® 
a The counting rate of the counter in 1958 a was transmitted 
-_ gontinuously, and the data were recorded only when the 
z- oe satellite was quite near one of the 16 receiving stations dis- 
“es tributed over the earth. 

The data collected by 1958 y were also telemetered con- 
tinuously. In addition, a small magnetic tape recorder stored 
oy the data obtained during each entire orbit. Then, as the 
satellite passed near one of the receiving stations, a radio 
command from the ground caused these data to be read out. 

A preliminary study of the data obtained from 1958 a@ and 
several interrogations of 1958 y has been carried out, with the 
following results. 

Reasonable cosmic ray counting rates have been obtained 
for altitudes below about 1000 km. In particular, we have 
obtained a plot of omnidirectional intensity vs. height in 
the vicinity of California for the first two weeks in February. 
. This curve, extrapolated down to altitudes previously 
reached by rockets, agrees with earlier data. 
altitudes greater than about 1100 km, v ery high counting 
rates were obtained. This conclusion is the result of a some- 
a what lengthy analysis. Geiger tube output rates up to about 
140/sec have actually been observed. In addition, periods 
have been found during which the geiger tube put out less 
than 128 pulses in 15 min. (We have a scaling factor of 128.) 
The considerations detailed in section 3 cause us to conclude 
that this is not due to equipment malfunction, but is caused 
by a blanking of the geiger tube by an intense radiation field. 
We estimate that if the geiger tube had had zero dead time, it 
would on these occasions have been producing at least 35,000 
counts/sec. 
: ey We surmise that the radiation we have found is closely re- 
e to the soft radiation previously detected during rocket 

flights in the auroral zone.” 

The radiation intensity necessary just to blank the geiger 
_ tube is equivalent to 60 mr/hr. In this connection the recom- 
_ mended permissible dose for human beings is 0.3 r/week.® 
The present radiation is 0.3 r in 5 hr or less. 
Several geophysical effects of this radiation seem possible. 

It is very likely closely related to aurorae and geomagnetic 
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storms. In addition, a rough calculation suggests that th 
radiation may be sufficiently intense to contribute ini portant 
heating to the upper atmosphere. It will be importaut to in- 
vestigate the amount of atmospheric ionization, licht and 
radio noise which would be produced, under various ::ssump- 
tions as to the nature of the radiation. 


a 1 Instrumentation for 1958 a and 


1958 vy 


The instrumentation for 1958 @ consisted essentially of a 
single Geiger Mueller tube, a scaling circuit for reducing the 
number of pulses to be worked with, and telemetry systems 
for transmitting the scaler output to the ground receiving sta- 
tions. The system contained in 1958 y was identical, with the 
addition of a miniature tape recorder for storing the data for 
the duration of each orbit and a command system to cause the 
telemetry of the stored information over a ground receiving 
station (Fig. 1). 

Identical G.M. counters, scaler input circuits and scaling 
circuitry were used in the two cases. The G.M. counters 
were Anton halogen quenched counters having approximately 
0.050 in. thick stainless steel walls. In addition, the counters 
were surrounded by the stainless steel cases of the payload, 
which were 0.023 in. thick. Thus the total absorption was 
approximately 1.5 gm em~? of stainless steel (approximately 
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weuT FROM FORK SCALER were reset, the transmitter turned off, and the next recording 
begun. 

et ial ; — The information thus telemetered to the ground was the 
- train of pulses emanating from the inhibitor circuit, except 

ee ee ee | that it was much compressed in time. It can be seen then that 
oe eee 4. scaler input pulse rates between 0 and 128 per sec were prop- 
erly passed on, and that all rates above 128 per sec appeared 
as a rate of 128 per sec, that is, all pulses missing. ns 

Fig. 2 Illustration of the function of the inhibitor circuit Table I is a list of the stations receiving data and reporting 


them to us. The stations labeled JPL are operated under the 
auspices of the Jet Propulsion Laboratory at Pasadena, Calif. 
Those labeled NRL are operated by the Naval Research 
Laboratory in Washington, D.C. Data were obtained from 
1958 a only when it was reasonably near one of these stations, 
since it had no provision for storing data for a later readout. 
We have already analyzed most of the data from the JPL 
stations, and some of that from the NRL stations as well. 
This work is continuing. 

A small magnetic tape recorder in 1958 y stored the cosmic 
ray information for an entire oribit, and then played it into 
a transmitter on command from the ground. Data from nine 
of these orbits have been reduced in a preliminary way. We 


iat the B75 per cent iron, 25 per cent chromium). The G.M. tubes had 
ortant sentially infinite lives, small variation in counting efficiency 
ttoin- § over the range —55 to 175 C, approximately 85 per cent 
it and # counting efficiency for cosmic rays, and about 0.3 per cent 
‘SuMp- F counting efficiency for photons of energy 660 kev. The dead 
time of the counters was approximately 100 microsec. The 
length of the counter wire was 4 in.; the inside diameter of 
the counter was 0.781 in. 

Following the counters were current amplifiers, which 
directly fed the first scaler stages. The scalers were bistable 
transistor multivibrators, which operated over a wide range of 


b's supply voltage and over a temperature range of —15 to 85 C. already have on hand many more of these passes, and are 
stems 4 tis limitation was caused by the supply batteries. The reducing the data from them in a routine way. 

ost. | Stler resolving time was 250 microsec. If input pulses It is evident from the above summary that the present re- 
+ ie at higher rates than 4000 per sec periodic were received, port is a very preliminary one. The nine cases from 1958 ¥ 
ta for | the Scaler simply indicated a constant rate of 4000 per sec. occur during the last four days of March, and we expect ulti- 
se the | That is to say, the sealer would not go out of operation if this mately to have data obtained during several weeks after those 
rim rate was exceeded. It did, however, have an input pulse days. In addition, we have so far reduced the data from 1958 
"6 f amplitude discrimination level, so that counter pulses of 7 only in a rather rough way, as explained in the following 
sali amplitude less than approximately one eighth normal were not paragraphs. Finally, we do not yet have highly accurate data 
mor counted. on the satellites’ orbits. We do have the position of 1958 aas 
ara In each of the satellites, the output of a scale of 32 was a function of time tabulated in 1 min interv als as supplied by 
hten 4 lemetered directly by the low power transmitter. In ad- the Vanguard computing center for the month of February. 
load dition, it was transmitted by the high power transmitter in These data seem to be in error by several minutes in time, but 
was | 1958 a. In all cases, the shift of state of the output scaler apparently are sufficiently accurate for the purposes of the 
ately | stages caused a discontinuous shift in the frequency of the present report. So far for 1958 y, we have only a set of orbital 
_— elements for March 26 and position vs. time for one orbit on 


subcarrier oscillators, of which the outputs were transmitted 
by the appropriate transmitters. The data telemetered in 
this manner have been readable when the rates of input pulses 
to the sealers were between 0.14 pulses per sec (16 pulses or 
one change of state per 2 min pass) and 80 pulses per sec 
(limited by the bandwidth of the receiving and data reduction 
systems). 

In 1958 y additional scaling circuits were included to pro- 
= §f vide a total scaling factor of 128 for the data to be stored. 
It was also necessary to include a time base, in order that a 
\RZED proper correlation could be established between the data and 
the satellite position. These two bits of information were 
combined in such a way that they could both be stored and 


April 1, together with estimates of the various perturbations. 
This information, supplied to us by the Vanguard Computing 
Center, has made it possible for us to estimate the orbit 
during the last days of March with reasonable accuracies. 
In particular, we estimate that our error in determining the 
time of passage through perigee is not more than about 5 
min on March 31, and is less on earlier dates. Our errors in 
estimating latitude and longitude may amount to 10 deg in 
some cases. 

Accurate orbital data will ultimately be supplied to us by 
the Vanguard Computing Center. 


telemetered on a single channel. Fig. 2 indicates the manner 
in which an inhibitor circuit effected this combination. The a : 7 
7 §f time base input was a train of pulses at the rate of one each Tebtet Receiving stations 
R sec. These pulses appeared at the output of the inhibitor, Blossom Point, Md. NRL 
D and were recorded, unless one was preceded by an output from Fort Stewart, Ga. be oe NRL 
=i the scale of 128, in which case it was suppressed. Antigua, Br. W. Ind. NRL 
The tape recorder was advanced in a discontinuous manner Havana, Cuba NRL 
; at the rate of one step per sec. As the tape advanced, it San Diego, Calif. NRL 
wound a spring for the eventual return of the tape to the Quito, Ecua. aeny 
startir : Lima, Peru NRL 
: 1g point. Antof, ta. Chil NRL 
| Upon receipt of a properly coded interrogation signal by coca. Chie 7” NRI 
the command receiver in the satellite, a relay system was Woomer, promt NRL 
" activated which caused the higher power transmitter to be Pudah Ale Poses Base, Fla. “JPL 
| turned on and the tape to be released, so that the spring was Earthquake Valley, Calif. JPL 
4 free to return it to zero. The return tape speed was con- Singapore = JPL 
$ trolled by an eddy current damping system, so that the play- Ibadan, Nigeria a) - JPL 
00, J back time was approximately 5 sec. As the tape returned, Temple City, Calif. : a PL 
the information was read off the tape, telemetered, and the Pasadena, Calif. JPL 
tape was erased. Upon completion of the cycle, the relays 
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T T T T T | t T We discuss first the data obtained from 1958 a avel 
105 «140 Fig. 3 is a plot of height against counting rate near the pare 
‘ California coast. All of the passes recorded by JPL stations ] sort 
aie ai pore a Ban in California are included in this graph. There is some § fron 
an Wiakteed ee variation in latitude, which presumably accounts for some of J mar 
‘ the seatter of the points. In addition, as explained above, the pow 
pi at orbital data are not yet known with good accuracy, and § ond: 
5 this presumably contributes significantly to the scatter. A § this 
8 | linear extrapolation down to a height of 100 km yields a value A 
of omnidirectional intensity of 1.22 (em*-sec)~, in adequate J we 
agreement with values we have previously obtained from rocket J mee 
flights, considering the crudity of the extrapolation. The data § it d 
Z 2 shown in this figure were nearly all taken before Feb. 11. rapi 
5 The data obtained by the NRL stations in South Ainerica rate 
is} | during the first two weeks of February are altogeth:r dif- § this 
ferent from those just discussed. The passes fall into two §& fal 
“f classes. In the first case, one obtains a counting rate of about | the 
30/sec, a roughly reasonable value. In the second case, the Fi 
telemetered signal fails to show a single scaler output pulse § tude 
ee 200 400 600 800 1000 00 ~~ 1400 1600 during the approximately 2 min of clean signal. This repre- hay 
. HEIGHT (km) sents an input rate to the scaler of less than about 0.1 /sec. & tran 
s Fig. 3 Counting rate vs. height near California for 1958 a There are, in addition, a few cases showing a strong change in & tuni 
' counting rate during the pass. miss 
| « For reasons discussed in section 3, we believe that the ex- § thre 
eee . z = tremely low output rate of the scaler is caused by very intense ff >15 
: radiation which “j jams” the geiger tube so that it puts out § dash 
re 6 pulses of such low i sight that they are below the threshold of § thei 
gue =F i the counting circuits. Laboratory tests show that this first J thes 
happens for the present equipment when the radiation the 
oF © ~, 30 COUNTS/SEC reaches such an intensity that a counter of the same effective peri 
dimensions and efficiency as the present geiger counter but corn 
‘ > .*° with a zero dead time would produce 35,000 counts/sec. high 
“ae Fig. 4 is a plot of height vs. geographic latitude in the § long 
bes vicinity of 75 W longitude. The positions of 1958 @ during § cour 
| reception of its telemetering signal by various of the NRL § 
; Li > x stations are marked. A code designates the kind of informa- § jn se 

rs tion received. It is at once evident that the extremely low 

1 — counting rates observed all occur at a high altitude, while 

(°) the more or less normal rates occur at a low altitude. Transi- 
4: fer of date tional cases occur at a low altitude. Transitional cases occur W 
from 1958 a over South America intermediat altitudes. : no 
Quite similar behavior is observed near Singapore, and pose 
probably also Ibadan. In these two cases no thorough study TI 


has been made, mostly because of the lack of trajectory data J may 
for the dates on which extremely low telemetered counting 
rates occur. In the one case at Singapore where such a rate reige 
occurred on a date for which orbital data were available, the 


= NO COUNTS OCCUR | 
— COUNTS OCCUR 


extremely low counting rate observed occurred at an altitude 
: of about 2000 km. able 

aan AD Fig. 5 is a plot of geographic latitude vs. geographic longi- J neq; 
tude for various orbits. Only the high altitude cases are oper: 
plotted on this figure. The fact that the segments of data do fF from 

— not correspond to positions of closest approach to the in- oper 

terrogating stations is due to our so far inaccurate knowledge tion, 

of the trajectory. whic! 

These data already suggest a picture of the geophysical temp 

phenomenon being measured. The data from 1958 y are § perat 

a much more explicit. Fig. 6 is a plot of the scaler output asa & carrii 

function of time as given by the tape recorder readout for Ar 

the pass ending near San Diego on March 28, 1748 UT. § thro: 


Since the tape recorder can only record one scaler output § This 
pulse each second (see section 1) themaximum indicationon the J one g 


x tape recorder output corresponds to 128 counts/sec for the § j 
ANTOFAGASTA | geiger tube output rate. (Our scaling factor is 128 in this A 
case.) It is evident from the figure that reasonable counting Th 
rates occur near the two ends of the pass. These ends cor- § geige 
| respond to the most northern latitudes and the lowest heights § effect 
an ons <. Toye) 90 above the earth. The section where the counting rate indi- § to ex 
cation is zero corresponds to a portion of the magnetic tape exper 
GEOGRAPHIC ig (°w where no tuning fork pulses were missing, and hence no scaler ‘¥ 
Fig. 5 Positions in latitude vs. longitude for telemetry of data output pulses occurred. This condition lasted 15 min, and 128 § whic} 
from 1958 a over South America pulses were fed to the scaler during this time. This is a0 § volta 
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average counting rate for the interval of 0.14/sec, to be com- 
pared with the usual cosmic ray rate for a geiger tube of this 
sort of about 50/sec. The counter goes through the transition 
from putting out essentially no counts to putting out a great 
many very quickly, and we presume that most of the 128 
counts observed during this 15 min interval occurred near the 
ends of the interval. There is, of course, no real evidence for 
this. 

As discussed in detail in the next section, we believe that if 
we had had a detector with zero dead time, and a storage 
mechanism of unlimited capacity, Fig. 6 would begin where 
it does now, and at about 13 min would have begun rising 
rapidly to a peak near 25 min at which point the counting 
rate would have been greater than 35,000 counts/sec. After 
this time, the rate would gradually have subsided, returning 
finally tv about the value actually recorded near the end of 
the pass. 

Fig. 7 is a plot of geographic latitude vs. geographic longi- 
tude of those orbits for which the tape recorder readout data 
have so far been analyzed. We have simply identified the 
transition points between portions of the record where no 
tuning fork pulses are missing, all tuning fork pulses are 
missing, or some tuning fork pulses are missing. These 
three different kinds of regions are identified on the graph as 
>15,000 ‘see, 128 to 15,000, and <128, respectively. The 
dashed portions of the various curves represent regions where — 
the identification as to counting rate range is uncertain. Since 
these passes all occurred during March 28 through March 31, 
the orbit did not have time to precess appreciably. Since 
perigee was near the most northern latitude, a given latitude 
corresponds closely to a given altitude. It is evident that at 
high altitudes and low latitudes, mostly in a certain range of 
longitude, the counting rate is very high. Near perigee the 
counting rate is low. Elsewhere intermediate counting rates 
occur. Possible interpretations of this result will be discussed 
in section 3. 4 


3 Interpretation of Observed Data 


We now propose to justify our claim that when essentially 
no sealer output pulses occur, the apparatus is, in fact, ex- 
posed to very intense radiation. 

Three possibilities are immediately evident. The apparatus 
may have some simple malfunction. This possibility can im- 
mediately be rejected except for the scalers, geiger tubes, and 
geiger tube voltage supplies, since the subsequent treatment 
of the information is completely different in the 1958 @ and 
1958 y. Some effect of temperature seems the only reason- 
able possibility here. The temperature of the geiger tube was 
measured in 1958 y and telemetered on the continuously 
operating transmitter. The observed temperatures range 
from zero to about 15 C. As discussed in section 1, the 
operating range of the circuitry is —15 to 85 C. In addi- 
tion, the frequencies of the continuously telemetering channels 
which carried the cosmic ray information are significantly 
temperature sensitive. These showed that no extreme tem- 
peratures occurred at the location of the corresponding sub- 
carrier generators. 

Another possibility might be that the satellite passed 
through regions which very few cosmic rays could reach. 
This is extremely unlikely. A magnetic field of the order of 
ohe gauss extending over thousands of kilometers and remain- 
ing unbelievably free of local irregularities would be required 
to exclude a sufficient fraction of the cosmic radiation. 

The possibility that we firmly believe is correct is that the 
geiger tube encountered such intense radiation that dead time 
fects reduced the counting rate essentially to zero. In order 
to explore this possibility, we have carried out the following 
experiments. 

A spare flight unit for 1958 @ was placed in an X-ray beam 
which was hardened by a 2 in. thick, brass absorber. The 
voltage on the X-ray tube was varied between 50 and 90 
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Fig. 6 A sample of the results of a tape recorder readout near 
San Diego on March 28, 1748 UT 
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Fig. 7 A plot of various orbits of 1958 ~ showing the range of 
counting rates as a function of position 
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key to vary the flux over a wide range. The counting rate 
was measured with and without lead shields which permitted 
only part of the beam to reach the geiger tube. In this man- 
ner the counting rates with and without the dead time effects 
were determined. As shown in Fig. 8, the dead time effects 
are negligible up to highest rates which can be handled by 
the telemetering systems. At high fluxes few of the pulses 
from the geiger tube have sufficient amplitude to operate the 
sealing circuit, and the counting rate returns to the range 
which can be telemetered. At very high fluxes no pulses have 
sufficient amplitude, and the counting rate is zero. 

An ion chamber placed in the position of the satellite ap- 
paratus measured an intensity of 60 milliroentgens per hr at 
the minimum flux required to reduce the counting rate to 
zero. The ionization produced by different energy X-rays or 
by charged particles producing this effect would of course be 
different from this measurement. The X-rays used for this 
measurement had energies in the range 50 to 90 kev. 

We have little concrete evidence concerning the nature of 
this radiation. Apparently, however, it is not electromag- 
netic. It makes its effects felt through the 1.5 g/cm? of ab- 
sorber which constitute the hull of the satellite and the walls 
of the counter. Photons with such energy should then be 
seen down to the lowest altitudes our equipment reaches. 
The radiation can presumably be either protons or electrons. 
If it is electrons, we then are probably detecting bremsstra- 
lung formed in the satellite shell. 


4 Implications 


Any reasonable identification of this radiation strongly sug- 
gests several geophysical consequences. It is unlikely that 
the particles have several Bev of energy each. Then in order 
to reach such low heights through the geomagnetic field they 
must at least initially be associated with plasmas which seri- 
ously perturb the magnetic field at an earth radius or so. 
We presume that this plasma is closely related to geomagnetic 
storms and aurorae. 


Secondly, at heights only a little above 1000 km, there jg 
still some atmosphere. Crude quantitative estimates suggest 
that the energy loss in this residual atmosphere of the radia- 
tion we detect may contribute significantly, if not domi- 
nantly, to the heating of the high atmosphere. In addition to 
considering this heating effect, it will be important to caley- 
late, on various assumptions as to the nature of the radiation, 
the amount of visible light, radio noise, and ionization pro- 
duced. 

Finally, there are obvious biological implications 0: these 
results. As discussed in section 3, if photons are being de- 
tected directly by the geiger tube, and if these photons are in 
the energy range 50 to 90 kev, then the radiation field ins ile the 
satellite corresponds to about 0.06 r/hr. The maximum 
permissible dose for human beings is 0.3 r/week. Other 
assumptions as to the nature of the radiation would 
obviously lead to different results. 
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Flame Stabilization in the Boundary Layer of Heated Plates 


This paper describes studies conducted with flames sta- 
bilized in the boundary layer of a heated flat plate. A 
graphical analysis of boundary layer flame stabilization is 
proposed and is verified by an experimental investigation of 
propane-air flames. The free stream flow was laminar 
with a turbulence intensity of 0.4 per cent, and the ap- 
proach mixture velocity was 25 and 50 ft per sec. The 
maximum Reynolds number encountered, based on the 
length of the heated plate, was 122,000. Surface temper- 
atures up to 2025 R were investigated. 


Introduction 


ig TURBOJET afterburners and ramijets, it is necessary to 
anchor the flame and prevent it from being blown out by 


the high velocity gases. One method which has been sug- 
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gested is the stabilization of the flame on a high tem- 
perature surface. Such a surface may have the form of a 
circular cylinder or a flat plate, aligned parallel with the flow. 
The phenomenon whereby a flame anchors itself within the 
boundary layer of the surface as a result of boundary layer 
configuration and surface temperature is called boundary layer 
flame stabilization. No flow separation occurs, and, in addi- 
tion, the surface upon which the stabilizing boundary layer is 
created need not obstruct the flow stream. No turbulence is 
created by the body other than that from the transition from 
the laminar to the turbulent regime. ae 


Previous Investigations 


Probably the use of heated bodies to stabilize flames was 
first proposed in 1950 by Tsien (1)* who suggested a flame- 
holder in the shape of a streamlined body. Such an airfoil 
would first be internally heated until combustion was initi- 
ated. Once the flame was started, it would be stabilized on 


’ Numbers in parentheses indicate References at end of p:per- 
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the surface of the streamlined body. Thus, the portion of the 
airfoil downstream of the flame front would be heated by the 
hot gases. By lengthening the portion of the airfoil exposed 
to the hot gases, the temperature at the front of the airfoil 
could be raised, and the approach mixture velocity could be 
increased delaying blow-off. 

Gross (2) studied boundary layer flame stabilization of pre- 
mixed propane-air flames on thin flat plates as well as on 
symmetrical airfoils. He determined maximum blowout 
velocities for flames stabilized in the boundary layer or in the 
wake at the trailing edge of the plate. 

A theoretical analysis of flames stabilized on a flat plate of 
constant temperature was made by Dooley (3) who formu- 
lated thie differential equations for laminar boundary layer flow 
in which heat is released by chemical reaction. These equa- 
tions were solved by mathematical iteration for any nonfrac- 
tional order of reaction. Applying the solution of a first order 
combustion reaction, the stabilization performance for azo- 
methaie flames was calculated. 

Gilbert, Haddock and Metzler (4) studied high heat release 
rates from the combustion of propane-air mixtures in a com- 
bustio. chamber consisting of an electrical heating element 
mounted concentrically inside a steel tube. They observed 
that a flame front became stabilized on the heated rod and 
that its position in the chamber shifted as the mixture veloc- 
ity, the equivalence ratio and the energy input to the heating 
element were varied. 

Hottel, Toong and Martin (5) examined the stability of 
lean propane-air flames in the boundary layer on a water 
cooled rod. For a given coolant temperature, they correlated 
the Blasius velocity gradient at the rod surface with the pro- 
pane-air ratio. 

Toong (6) studied the incipient ignition of ethanol-air 
mixtures flowing over a hot surface. He made qualitative 
theoretical predictions of incipient ignition from boundary 
layer theory and correlated these experimentally. Finally, 
ina more recent study, Toong (7) studied the mechanism of 
flame stabilization by boundary layer theory, 


Proposed Model 


In this paper boundary layer flame stabilization by hot sur- 
faces is explained by the hypothetical model shown in Fig. 1. 
The proposed model consists of a flat plate over which a com- 
bustible mixture flows steadily, and a laminar boundary layer 
isereated. A flame propagates upstream through the bound- 
ary layer on the plate until the upstream position is reached 
where the boundary layer is sufficiently thin to quench any 
further penetration of the flame, and the flame front becomes 
stabilized at this position. 

If the flowing fluid is a combustible mixture, it has associ- 
ated with it a finite velocity of laminar flame propagation. 
Since the fluid velocity tends to vanish as the surface is ap- 
proached, there will be a region near the surface where the 
flame velocity exceeds the mixture flow velocity. If a flame 
is initiated in this region, the flame can propagate upstream 
through the combustible mixture. On the other hand, the 
surface of the plate acts as a heat sink as well as an absorber 
of active chemical radicals, thereby inhibiting the combustion 
reaction in a layer of the mixture adjacent to the surface. In 
this layer, near solid surface, there exists a zone devoid of 
flame which is called the dead space.4 Smaller in depth, the 
dead space lies within the boundary layer. A flame prop- 
agating in the boundary layer region can occur no closer 


_ ‘The dead space is defined as the depth within which the sur- 
lace extinguishes the flame (8). Immediately above the dead 
space is a zone in which the combustion reaction is inhibited, re- 
sulting in a decrease of the laminar flame velocity. The depth 
at which the velocity of flame propagation is 99 per cent of the 
velocity at an infinite depth of combustible mixture and which 
includes the dead space, is defined here, as quenching distance. 
Although this definition is an arbitrary one, it is based on the 
precedent used in defining boundary layer thickness. 
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to the surface than the dead space distance. As the flame 
propagates upstream, the boundary layer through which it 
travels becomes thinner as the leading edge of the plate 
is approached. The flame traveling through the bound- 
ary layer can move upstream only as long as the flame veloc- 
ity exceeds the fluid velocity at a distance from the surface 
greater than the dead space. The distance normal to the 
surface at which the mixture velocity is equal to the flame 
velocity of the mixture is an important parameter, because 
above this locus of equal velocities the fluid velocity is greater 
than the flame velocity, and hence any flame in this region 
would be carried downstream. Below the locus of equal 
velocities, the flame velocity exceeds the fluid velocity, and 
any flame in this region would propagate upstream parallel to 
the direction of flow. From the physics of the phenomena, it 
may be predicted that the limiting upstream position of the 
flame front should correspond to the location along the plate 
where the locus of equal velocity intersects the dead space, 
because at this position quenching will prevent the flame from 
penetrating any further upstream. The flame front should 
stabilize here and propagate into the fresh mixture. 

When the surface temperature is raised by external or in- 
ternal heat addition, several important changes take place. 
The flame velocity in the boundary layer will increase because 
of the higher local mixture temperature. The boundary layer 
thickness will increase because, with increased heat transfer 
from the surface to the boundary layer, the rate of boundary 
layer growth is more rapid. The dead space thickness will 
decrease, because as the surface temperature increases less 
quenching occurs. All of these changes associated with 
higher surface temperatures allow stabilization at a higher 
free stream velocity. For a given surface temperature, the 
boundary layer thickness may also be increased by lengthen- 
ing the heated plate, thereby permitting flames to be stabilized 
at higher mixture velocities. 


Graphical Solution 


It will be shown in this paper that it is possible to determine 
the position of the stabilized flame on a flat plate if the locus 
of equal velocities and the depth of the dead space can be 
specified for the desired conditions of surface temperature 
distribution and free stream velocity. 

In summary, the graphical solution is based upon two basic 
combustion properties, the laminar flame velocity and the 
dead space distance. The analysis does not require a knowl- 
edge of the steps of the chemical reaction, the order of the 
reaction, or values of the activation energy, the diffusion co- 
efficient or the frequency factor, which often are not known. 
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Fig. 2 Boundary layer velocity profiles for a flat plate of uniform 
temperature, 7 = 0.050 sec 


In developing the graphical solution, the following as- 
sumptions are made: 

1 No appreciable chemical reaction occurs upstream of the 
flame front. 

2 The boundary layer temperature profile and the bound- 
ary layer velocity profile upstream of the flame front are not 
appreciably affected by the flame. 

3 The mixture viscosity varies with the static temperature 
as defined in (8). 

4 The curvature of the flame front within the quenching 
zone can be neglected. 

5 The specific heat of the mixture is a constant. 

The advantage of the proposed method is that the following 
assumptions, common to numerical solutions, need not be 
made: 

1 The surface has a particular temperature distribution. 

2 The Prandtl number equals the Schmidt number, and 
both are equal to unity. 

3 The chemical kinetics is expressed by the Arrhenius re- 
action rate equation or some similar relation. 

4 The order of the combustion reaction is nonfractional. 

5 The average molecular weight of the mixture is constant. 

The proposed graphical construction requires knowledge of 
the following functions: 

1 The variation of temperature and velocity in the fluid 
boundary layer along the plate surface. 

2 The variation of the laminar flame velocity with mixture 
temperature. 

3 The variation of the depth of the dead space with plate 
surface temperature. 

The methods of developing the graphical functions and the 
manipulations required to obtain the final solution of the 
analysis will be exemplified through consideration of a stoi- 
chiometric propane-air mixture for various plate temperatures, 
and for an initial mixture temperature of 300 K. A uniform 
temperature distribution along the plate will be assumed, al- 
though this assumption is not necessary to obtain a solution 


of the problem. First, each of the three previously listed a 
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graphical functions will be determined, and a discussion of the 
necessary assumptions and approximations will be under- 
taken. Then, the technique of utilizing the functions to pre 
dict the flame stabilization performance will be presented. 


Velocity and Temperature Profiles 


The first functions of the illustrative solution to be calcu- 
lated are the fluid boundary layer velocity and temperature pro- 
files. These may be determined analytically by the general 
boundary layer solution of Chapman and Rubesin (9). 
Their solution applies to a laminar compressible boundary 
layer on a flat plate having an arbitrary temperature distribu- 
tion. For a stoichiometric mixture of propane and air, the 
fuel represents only 4 per cent of the mixture by volume, 
and hence it is assumed that the fluid mechanic properties of 
the combustible mixture are the same as for air. The deter- 
mination of the fluid boundary layer profiles by the above 
method does not take into consideration any chemical change 
in the fluid, and it is assumed that the degree of chemical 
reaction within the combustible mixture is negligible between 
the plate leading edge and the flame front. This is a valid 
assumption, because Linnett (10) has shown that the thickness 
of the reaction zone in laminar combustion is of the order of 
one millimeter, and hence is small compared to the length o! 
surface being considered. 

The boundary layer temperature and velocity profiles were 
calculated, using the solution of Chapman and Rubesin, for 
plates of uniform temperature distribution in the range of 
temperatures from 500 to 1100 K at increments of 100 deg. 
The profiles were computed for several values of the ratio 
7 = 2/U,. It is then necessary to plot the distance y versus 
the dimensionless velocity u/U, and the fluid temperature 7, 
respectively, for different values of r and 7». In this study, 
values of + between 0.001 and 0.100 sec were selected because 
they represent practical values of plate lengths for the pat 
ticular mixture velocities encountered. In Figs. 2 and 3, 8 
typical set is shown for r = 0.050 sec. Similar curves were 
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jotted for other values of r but they are not shown. The use 
of the computed boundary layer profiles will be described 


later. 


Laminar Flame Velocity 


The second graphical function, namely the variation of the 
laminar flame velocity with mixture temperature is required 
next. At a given pressure and for a given combustible mix- 
ture, the velocity of laminar flame propagation is a function 
of the mixture temperature, and this relation is readily ob- 
tained from experimental investigations. 

In addition to being a function of the mixture temperature, 
the local laminar flame velocity is influenced by the quenching 
effect of a solid surface (8). However, it can be shown that 
this effect may be neglected in the analysis. Consider, for 
example, the flame velocity profile for a laminar flame front 
traveling through a quiescent isothermal mixture over a solid 
surface. The reduction of the flame velocity in the region 
within the quenching distance but above the dead space 
causes thle moving flame front to be curved rearward. Since 
the fam: front travels as a whole and since the displacement, 
parallel to the surface, of the flame front within the quenching 
distance is small compared to the length of the flat plate, the 
flame front can be considered to have the same effect as a 
plane flume front extending toward the surface to the border 
of the d-ad space. The resulting uniform flame velocity pro- 
files are assumed to vanish discontinuously at the border of 
the dead space. A similar situation exists if the mixture is 
flowing across the surface and the flame front is fixed relative 
tothe plate. Consequently, the flame velocity can be consid- 
ered as 2 function of the local mixture temperature only, and 
the inhibiting effect of quenching above the dead space can be 
neglected. The latter assumption greatly simplifies the solu- 
tion of the theoretical analysis. 

For the illustrative solution, the variation of flame velocity 
with mixture temperature for a stoichiometric propane-air 
mixture is required. Experimental data from Dugger (11, 12) 
are plotted in Fig. 4. Since both sets of flame velocity data 
coincide up to a temperature of 600 K, it is reasonable to as- 
sume that the curve for the maximum flame velocity ade- 
quately represents the flame velocity variation for an equiva- 
lence ratio of unity in the range of 600 to 1000 K. This rela- 
tion is the second graphical function which is required for the 
proposed solution 


Déad Space 

The third function required for the proposed analysis is 
the relation between the depth of the dead space and the plate 
surface temperature. Values for the dead space distance for 
various surface temperatures are best obtained from experi- 
mental measurements. The dead space distance is closely 
related to the quenching distance, and, according to Lewis 
and von Elbe (8) and Wohl et al. (13), the dead space depth 
should be about one half of the quenching distance. More- 
over, it has been shown by Friedman and Johnston (16) that 
the composition of the quenching surface has no measurable 
efiect. Because in this study the quenching distance above a 
single flat plate is required, the values of quenching distance 
4 commonly measured between parallel plates are halved, 
and hence the dead space distance above a single surface is 
approximately one fourth of the quenching distance for paral- 
lel plates. 

When experimental data are lacking or available data must 
de extrapolated to a higher range of temperature, semithe- 
oretical or empirical equations may be used. Potter and 
Berlad (14, 15) have derived equations for the flame quench- 
ing distance d in terms of the temperature. Their empiri- 
tal equation reads ( 


sod = 1.206 X 10-5 em?/sec 
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where 
8s. = burning velocity 
7, = mixture temperature and also the surface temperature © 
T; = adiabatic flame temperature 


The equation applies to stoichiometric hydrocarbon-oxygen- 
nitrogen flames at one atmosphere pressure. Since burning 
velocities are generally known for a wide range of tempera- 
tures, the quenching distances may be calculated for inter- 
mediate and high temperature surfaces. 

For the problem under consideration, the variation of the 
dead space distance with surface temperature for a stoichio- 
metric propane-air mixture is required. Such a relation may 
be partially derived from the experimental values of quench- 
ing distance for a stoichiometric propane-air mixture reported 
by Friedman and Johnston (16). The dead space distance, 
however, must be known for temperatures considerably above 
the range tested by Friedman and Johnston. Therefore, 
Potter and Berlad’s expression for quenching distance based 
on the burning velocity is used to obtain the dead space dis- 
tance in the range of surface temperatures from 560 to 1000 
K. Thisis shown in Fig. 5. This curve is the third graphical 
function required for the graphical solution. 

The basic graphical functions have been developed and the 
graphical solution may now be obtained. The next step of 
the solution is to determine the locus of equal fluid velocity 
and flame velocity. The boundary layer fluid temperature 
profiles have been plotted in Fig. 3, and since the relation 
between the fluid temperature and the laminar flame velocity 
is known from Fig. 4, a flame velocity profile can be con- 
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Fig.6 Boundary layer flame velocity profiles for a stoichiometric 
_propane-air mixture and for a flat plate of uniform temperature, 
t = 0.050 sec, U, = 25 ft/sec 


0.40 


ry corresponding to each particular fluid temperature 
profile. In other words, at every point in the fluid boundary 
. layer, the local fluid temperature defines a corresponding local 
flame velocity. A plot of typical flame velocity profiles is 
shown in Fig. 6. 

After the flame velocity profiles throughout the boundary 
ayer are known, each is superimposed graphically upon the 
- corresponding fluid velocity profile. The intersection of the 
two profiles defines the point of equal velocities. This is 
shown in Fig. 7 for one particular value of r. However, 

similar curves for the other values of 7 must be constructed so 

that the locus of equal velocities throughout the boundary 

layer can be determined. When this has been done for a 

sufficient number of points to define a curve, the resulting 

curve is plotted, as in Fig. 8, and the graphical relation be- 
tween the dead space depth and the surface temperature is 
superimposed. The intersection of the locus of equal veloc- 
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ities for each value of 7, with the curve of the dead spac 
distance, represents the location in the boundary layer wher 
conditions are satisfactory for a flame front to become sts- 
bilized. 

The resulting performance map, defined by the graphical 
solution, is plotted in Fig. 9 for two mixture velocities. 


In 
stabil 
heate 
exper 
proac 
plate 
tion 
tance 
front 
ture 
propa 
upon 
graph 
The 
heatec 
upon 
A des 
follow: 


The 
2'/ in 
the tes 
shown 
makin, 
chamb 
plate 
The he 
flat pla 
The 
Heat t 
diagon 
by pas 
were ol 
in an 
chambr 
£ases 
directic 
the hor 
propan 
egulat 
The 
removii 
stream 
a tr: 
The bo 


In or 


2200 
2000 
1800 o| 
cre 
> 
& 1600 
: Ng 
1200 rw 
8 
= 
| 
1000 ° ~ 
< 900 | 
® 800 | ch 
| 
© Up = 25 FT/SEC. H 
700 U.=50 FT/SEC. 
woe | 
3 45678910 2 3 456789 


DISTANCE TO FLAME FRONT, INCHES 


Fig. 9 Comparison of theory and experiment for stoichiometti 
propane-air flames stabilized on a heated flat plate of uniform 
temperature 
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— om proximately 0.4 per cent. The mixture was not preheated 
oa TEST CHAMBER and was at a temperature of approximately 535 R for all 
— Flow tests. Commercial grade propane having an analysis of 95 
lee 3 PE aud per cent to 98 per cent propane with 2 per cent to 5 per cent 
— UNIOK other saturated hydrocarbons was used in the experimenta- 
iV RECORDING tion. The velocity profiles in the test chamber were found 
K De aren —" to be flat to within } in. from the floor of the chamber at 
uct co velocities of 26.6 50.7 ft per sec. At in. from the 
7 aie surface, the velocity was 95 per cent or greater of the free 
8 = stream value. The distance from the leading edge of the 
24 TO ORIFICE {7 plate to the flame front was determined photographically by 
> TO ROTAMETER onmen scaling. The photographs were taken with an exposure of 
SWITCH 3 sec at f 4.5. 
Fig. 10 Schematic diagram of test apparatus 
= 
: - — The tests of flame stabilization characteristics on a heated 
Work flat plate were performed according to the following procedure. 
The main air flow rate was first adjusted so that, when the 
—7 In order to verify the solution of the graphical analysis, the proper amount of fuel was added to the air a predetermined 
stabilization of propane-air flames in the boundary layer of a mixture ratio and the desired free stream velocity would occur 
heated flat plate was investigated experimentally. All the in the test chamber. The boundary layer suction was initiated 
aperimental data were obtained for constant mixture ap- and set to the required flow rate for the mixture velocity. 
proach velocities of 25 and 50 ft per sec and for a range of The air and fuel for the heating combustion chamber were 
plate temperatures from 1000 to 2000 R. The flame stabiliza- turned on and ignited by a high voltage spark. The second- 
— tion performance was measured as the variation of the dis- arv air was then turned on to dilute and reduce the tempera- 
tance from the flat plate leading edge to the stabilized flame ture of the combustion gases. The fuel and secondary air 
— Ff front with change in the plate temperature for a constant mix- were adjusted to give the desired uniform plate temperature. 
ture velocity. The experimental results for stoichiometric The next steps of the procedure were carried out in rapid suc- 
~ propane-air mixtures are plotted in Fig. 9, and superimposed cession. The propane was turned on to a predetermined 
050 upon the theoreti “al performance curves predicted by the flow rate which would give the desired equivalence ratio. The 
graphical analysis. : ; mixture was then ignited by a high voltage spark jumping 
y pte The experimental apparatus was designed to provide & between the chamber floor and a stainless steel rod which was 
heated flat plate with a uniform temperature distribution positioned just inside the chamber exit. If conditions were 
upon which premixed propane-air flames could be stabilized. favorable for the stabilization of a flame on the heated sur- 
A description of the experimental equipment and procedure face, the flame would travel upstream from the spark rod 
space follows. near the chamber exit to a position on the heated plate. If 
vhere conditions were not satisfactory for stabilization of a flame 
ste Equipment on the plate, the flame would be stabilized only by the spark 
discharge. If the flame became stabilized on the heated 
hia The sith imental investigation was perfo rmed = 21/2 by plate, (which would occur within a fraction of a second after 
2/, in. cross-sectional test chamber. A schematic diagram of the igniting spark was turned on) the flame was quickly pho- 
the test equipment ’ shown in Fig. 10. The test chamber, tographed through the glass walls of the test chamber. The 
shown ™ Fig. 11, had vertical glass walls of Vycor plate glass, fuel was then shut off, extinguishing the flame..-The flames 
making it possible to observe the combustion process in the were allowed to remain in the combustion chamber only long 
chamber. The horizontal chamber floor consisted of a heated enough to be photographed. This procedure prevented the 
plate of 5 in. length, upon which the flames were stabilized. flame from changing the plate temperature appreciably. 7 
The heated section of the floor of the test chamber simulated a This test procedure was repeated for a series of tests within a 
fat plate in a semi-infinite flow field. ; ; : range of plate temperatures of 1010 to 2025 R, for approach 
The heated plate was made from } in. thick stainless steel. mixture velocities of 25 or 50 ft per sec and for various mixture 
Heat transfer fins were formed on the bottom side by milling ratios within the stable combestion Hiaiie. 
diagonal grooves to a depth of } in., and the plate was heated The time rate of change of the plate temperature, at a posi-_ 


ty passing hot gases along the underside. The hot gases 
were obtained by the combustion of premixed propane and air 
ian auxilliary vertical combustor mounted below the test 
chamber. A horizontal rectangular duct carried the hot 
gases along the bottom of the heated plate parallel to the 


: (rection of the mixture flow in the test chamber. Control of 


the hot gas temperature was obtained by regulation of the 
propane and primary air flow rates and by the addition and 
iegulation of secondary air. 

The boundary layer profile of a flat plate was obtained by 
rmoving the boundary layer previously formed on the up- 
steam portion of the test chamber floor by applying suction 
0a transverse slot at the leading edge of the heated plate. 
The boundary thickness was measured with a probe. 

In order to reduce the turbulence level in the test chamber, 
the propane-air mixture was passed through straightening 
vanes, a diffuser, six turbulence screens of 200 mesh wire cloth, 
‘calming chamber and then a converging nozzle with a 23 to 


stritf | contraction ratio before entering the test chamber. The “Ss 
omf {tee stream turbulence intensity in the test chamber was de- Fig. 11 Test chamber in operation showing flame stabilized on _ 
‘mined by hot film anemometry and was found to be ap- heated flat plate a : 
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Zz 50 rived from the quenching distance measured by Friedman and U1 
4 Johnston (16) exhibits a minimum value at an equivalence § ynbu 
& ratio of approximately 1.1. Since a shorter heated length js prac 
3 7” required to stabilize a flame as the quenching distance is mr. ff srea 
oe oe © N:., duced, the heated length would tend to be a minimum for g § ther 
mixture having the minimum quenching distance. In addi. § toin 
7 tion, Zukoski (17) has found that for propane-air mixtures socia 
under laminar conditions preferential diffusion causes the § sttai 
st mixture at the flame interface to be leaner than the free 
on | ag ot stream. Thus, the mixture in the flame region is «ctually 
leaner than the overall equivalence ratio, causing }erform- Th 
ance to peak at an even richer mixture. All these the ( 
gs factors in turn tend to allow the flame to be stabilized with a § 
7 4 minimum heated length with fuel rich mixtures. a 
Observation of the stabilized flame front showed oscil- 
latory motion to exist. The flame front appeared as a blur yom 


ALENC 
to the eye, but the downstream position was somewh:it more 


discernible, because here the flame front mom: ntarily 
stopped before reversing the direction of travel. Such oscilla- 
tion of the flame front is characteristic of boundary layer 
flame stabilization and had already been observed by Dooley , : 
and by Gross. Dooley reported that the oscillation is more § °° 
2 no. 81 
pronounced if the heated surface has a slightly negative tem- 9 
perature gradient, i.e., the temperature decreases dow) stream p 


Fig. 12 Effect of propane-air equivalence ratio on distance from 
plate leading edge to stabilized flame front 


along the body. In the experiments described in this paper, 3 
it was observed that flames stabilized on the portion of the § Disse 
plate having a small positive temperature gradient were stable 4 
and exhibited little oscillation. Comt 

In an attempt to gain an insight inte this oscillatory phe- tabi 


nomenon, high speed schlieren motion pictures were taken, and J) ine 
frame by frame analysis of the film was made. For a free J ''¢ 
stream velocity of 32 ft per sec and an equivalence ratio of J Roun, 
0.902, the frequency of individual cycles varied from 14 to § patio 
30 cycles per sec and appeared to be a random variation. 
The amplitude of the oscillation varied from 3 in. to 2 in. 
and did not seem to correlate in any way with the frequency 
of the cycle. The ratio of the average upstream velocity of 
the flame front to the average of the downstream velocity of 
the flame front varied from 0.5 to 2.0 and did not seem to be 
influenced by the frequency. The magnitude of this ratio 


Fig. 13 Typical schlieren picture of propane-air flame stabilized 
on a heated flat plate, U, = 34 ft/sec, 7,, = 1020 F, ¢ = 0.902. 
Flow from left to right 


tion approximately an inch downstream from the flame front, 
due to the presence of the flame, was found to be of the order 
of 5 R per sec at a plate temperature of 1460 R. Because the 
period of time during which combustion occurred in the test appeared to be random within the extremes. These observa- 
chamber was very short namely, about 2 or 3 sec, the plate tions are in disagreement with those of Gross who reported 


temperature did not.change appreciably during each test. the upstream movement of the flame front to be four times as 
fast as the downstream motion. 


The proposed model for boundary layer flame stabilization 


Experimental Observations I 

: may be used as a basis to hypothesize a cause of the flame F 4, 
In addition to obtaining data for the stabilization perform- oscillations. Since, according to the proposed model, the § or... 
ance of propane-air flames on the heated plate, the experi- point of stabilization is defined by the intersection of two fy, g 
mental investigation allowed observation of several charac- curves, with very nearly equal slopes at the point of inter §, 4,,, 
teristics of boundary layer flame stabilization. Specifically, section, it is probable that the oscillation is a hunting motion J ,,,.,. 
it was observed that the position of the stabilized flame on the of the flame tip moving along the boundary layer in a range § 4, 5, 
heated plate was affected by the fuel-air ratio of the flowing where approximate conditions for the stabilization of the J ,...,, 
mixture. Also, a characteristic oscillation of the stabilized flame exist. If the stabilizing surface has a negative tem- Bj... 
flame front was noted in all the tests. perature gradient, the curves intersect at a greater angle. etal 
Although less pronounced at higher plate temperatures, the Physically, this means that the conditions necessary for the F pach 
equivalence ratio of the approach mixture had a considerable stabilization of a flame are more localized and definite. This belie 
effect on the position of the stabilized flame for any given should reduce the distance along which conditions may be § ,... 

plate temperature and approach velocity. The variation of favorable for the flame and reduce the flame oscillation. This 

flame attachment with equivalence ratio is shown in Fig. 12 effect is in agreement with experimental observations by 

for free stream velocity of 50 ft per sec. It may be noted that Dooley. 

the minima of the curves consistently fall within the range of Also it was observed that the rate of flame spreading across T" 
equivalence ratios of 1.15 to 1:25. This phenomen is in agree- the test chamber was low. This would be expected, since pl 
ment with the experimental results of Gross (2) and Gilbert, there is no obstruction in the stream to generate turbulence study 
Haddock and Metzler (4) in which the maximum blow-off as does a bluff body flameholder. The flame spreading results phenc 
velocities for a given stabilizing surface occurred at equiv- from essentially laminar flame propagation into the free Visg 
alence ratios of about 1.25. stream from the flame tip stabilized in the boundary layer. § “ste 
The effect of the equivalence ratio can be explained by the This is indicated by high speed flash schlieren pictures, of In the 
facts that the minimum required heated length would tend to which Fig. 13 is typical. The flash duration was approxi under 
occur at equivalence ratios at which the flame velocity is a mately 7 microsec. Viewing the slightly wrinkled flame sheet the th 
maximum, and that the flame velocity, according to the ex- from the side causes it to appear on the photograph as a clU* F Rec 
perimental work of Dugger (11), peaks at an equivalence ter of lines. The boundary layer suction slot is located at the J ' As 
ratio of about 1.1. Moreover, the dead space distance de- left edge of the photograph. *N 
Sepry 


598 


JET PROPULSION 


j 
sre 
nat 
~ 
1 
& 
SNE 
» 
: 
> 
= 
: 
> 


and 
lence 
rth Is 
is re- 
for a 
addi- 
tures 
3 the 

free 
ually 
‘orm- 
these 
ith a 


oscil. 

blur 
more 
arily 
cilla- 
layer 
oley 
more 
tem- 
ream 
sper, 
the 
‘able 


phe- 
and 
free 
0 of 
4 to 
tion. 
2 in. 
ol 
y of 
o be 
‘atio 
rva- 
rted 


S aS 


tion 
ame 
the 
two 
ter- 
tion 
nge 
the 
em- 
gle. 
the 
‘his 
be 
‘his 


by 


Under laminar conditions, the travel of the flame into the 
yburned mixture is relatively slow. However, in most 
practical applications of flame stabilization in high speed 
reams, the turbulence of the approach mixture would cause 
the rate of flame propagation across the combustion chamber 
toinerease. On the other hand, the rapid flame spreading as- 
geiated with bluff body flameholders would not normally be 
attained with boundary layer flame stabilization. “oa 
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n attempt is made to obtain a physical understanding 
oferosive burning of a solid propellant, taking into account 
the flame structure. At low pressure the burning rate of 
adouble base propellant is entirely determined by the heat 
transfer within the fizz zone with no heat exchanged from 
the high temperature flame zone. Therefore it has been 
assumed that erosive burning is due to the effects of turbu- 
lence within the fizz zone, the threshold velocity corre- 
sponding to the moment when incipient turbulence 
reaches the boundary of the fizz zone. By so doing it has 


agreement with the experimental values. 


Introduction and Basic Assumptions 


_— burning rate of a solid propellant is increased by the 
presence of a flow parallel to the burning surface. The 
‘tudy of experimental data suggests, moreover, that the 
bhenomenon takes place only when the flow mean velocity 
U is greater than a so-called “threshold veloc sity” Ui. (The 
existence of the threshold velocity is not clearly established.) 
lathe present paper, an attempt is made to obtain a physical 
inderstanding of erosive burning, assuming the existence of 
the the threshold velocity. 
Jen. 7, 1958. 
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hen possible to calculate an erosion coefficient in close 


1954, pp. 9-12. 


In this problem in which aerodynamics, heat exchanges — 
and chemical kinetics are simultaneously involved, it has 
been possible to show the influence of various parameters and 
even, despite numerous drastic assumptions, to find numerical 
results in accordance with the experimental values. The 
study of the combustion of a double base propellant is in 
itself a fairly difficult problem involving very complex 
theories. 

The present paper is restricted to double base propellants 
burning at low pressure. The numerical values used in the 
example are the results of elaborate calculations made by 
Rice and Ginell (1)? for the HES. 4016 propellant. 


Basic Assumption in Solid Propellant Combustion 


During combustion, the solid is transformed into a gas 
by chemical reactions (see Geckler (2)). Even if the trans- 
formation occurs through the formation of an intermediate 
liquid phase, there exists a certain position, called the burn- 
ing surface, where gas is being formed. The mechanism of 
the solid-phase decomposition is considered by the majority 
of authors as the rate-controlling step of propellant burning, 
and the gas is being formed at a rate 
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ot Fig. 1 Representation of the flame of a double base propellant 
a burning at fairly low pressure 
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where K is a constant, B the frequency factor, E the activa- 
tion energy, m the mass rate of burning per unit area 
(m = r pp), R the universal gas constant and 7’, the tempera- 
ture of the gas at the burning surface. 

The reaction occurring in the solid phase would have in- 
-_ ereased the temperature at the burning surface to a value 
ae T;’, but actually a temperature 7, greater than T',’ is ob- 
7 tained there because heat is conducted back to the solid from 
the gas. The essential theoretical problem, according to this 
view, is to discover how 7, depends upon the gas phase 
reactions and behavior. 

We have assumed that the problem of combustion without 
erosion has been satisfactorily solved (see (1, 2)). There- 
: fore the purpose of this study of erosive burning is to estab- 
>. lish the influence of a flow parallel to the burning surface on 
temperature T,,. 
ee? Fig. 1 represents a model of the flame of a double base pro- 


pellant, satisfactory from both the experimental and theo- 

retical points of view. After the solid phase decomposition, 

_ the gas phase reactions can be divided in three parts: Close 
to the burning surface, some exothermic reactions take place 
in the fizz zone. Then, in the preparation zone, activated 
products are formed without heat production. Finally, 
when a sufficient concentration of activated products is 
achieved, the final reaction occurs in the combustion zone at 
the end of which the isobaric combustion temperature is ob- 
tained. 

At pressures between 20 and 80 atm, the thickness of the 
preparation zone is much greater than the thickness of the 
other zones, and the heat transfered from the flame zone 
into the fizz zone can be neglected. Therefore the temper- 
ature 7’; at the end of the fizz zone is approximately equal to 
the temperature 7,’ obtained without any heat exchange 
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across the preparation zone, and temperature 7’, can be cal. 
culated considering only the heat exchange across the fiz; 
zone. 

This picture is probably fairly accurate when limited ty 
the case of an end-burning strand where the gas flow js 
laminar and one-dimensional, directed normally away from 
the burning surface. In the more complex case of an jm. 
portant velocity parallel to the surface, one can raise the 
question of the validity of this model. We have assumed 
that the presence of a strong velocity gradient parallel to the 
burning surface (boundary layer) does not alter the picture 
appreciably in the laminar sublayer where all transport of 
heat or chemical species is accomplished by movements on q 
molecular scale. However, in the turbulent layer where the 
eddies play an important role in the transport, it is clear that 
the picture of separate zones does not hold anymore. Mean. 
while one can expect that the flame structure still holds in 
the transition sublayer where turbulence is not yet fully de. 
veloped and where the dimensions of the eddies are much 
smaller than the thickness of a zone. 


Conduction Equation Across the Fizz Zone 
Without Erosion 


At a distance y; from the burning surface, the fizz zone re- 
action is completed, and therefore the temperature 7,’ is 
achieved if no heat is transferred back from the flame zone. 
A certain heat release pattern must be assumed in the fizz 
zone. Following Rice and Ginell (1) it was supposed that 
the heat evolved in the fizz zone is all released at the sam 
distance y;. Meanwhile provision must be made for diffusion 
and for distributed heat sources by means of a coefficient ( 
Therefore the equation giving the actual temperature 7, o! 
the burning surface can be written 


Where €, is the mean specific heat at constant pressure, and 
) the mean thermal conductivity across the fizz zone, m being 
the mass burning rate. 

For a given propellant the quantities 7, 7,’, 7,', m, ¢, 
X and y; can be computed as functions of one independent 
variable, the pressure, for instance. Therefore the pressur 
dependent value of C can be deduced from Equation [2]. 

Actually it is much simpler to consider as a whole the group 
Cz,/X, and for our present purpose it is useful to plot that 
term against y; from the results obtained by Rice and Ginel 
for HES. 4016 (see Fig. 2). 


Conduction Equation Across the Fizz Zone 
With Erosion 


In the presence of erosion due to the flow parallel to the 
burning surface, turbulence takes place in the combustion 
zone, increasing the thermal conductivity at each point by 
factor taking into account the turbulence level at that point 
Therefore the rate of heat transfer increases together wit! 
the temperature 7, at the burning surface and the mas 
burning rate m, when the thickness of the combustion zone 
shrinks. 

Indicating the values of 7’, and m without erosion by 7’,. and 
m,, Equation [1] can be written 


B 


Within the fizz zone the temperature 7’ is achieved at 4 
point yi’ < y: where the fizz zone reactions are possibly not 
yet completed. The conduction equation can be writtel 
directly between that point y,’ and the burning surface, acros 
a region in which the flow is still completely laminar 
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From Equation [4] the value of 7’, and consequently of 
wm, can be obtained for each value of y:’ < y: if the proper 
value of C or C is used. 

aan [4] is valid in the laminar sublayer, since the 
presence of a velocity gradient perpendicular to the burning 
wrface can hardly be expected to have a significant influence 
» the temperature distribution. That assumption is cor- 
wborated by the fact that no variation of the temperature 
iistribution (no erosion) seems to exist below the threshold 
elocity, which can be viewed as the velocity at which the 
incipient turbulence penetrates into the fizz zone. 

A difficulty arises with the choice of the proper value of 
('¢,/\ to be used in the laminar sublayer between the burn- 
ag surface and y,’ if the outer part of the fizz zone is in this 
turbulent region. The value of C ¢,/ taken at the point 
, without erosion would be too small and would not account 
jor the increase in the temperature level and for the asso- 
ited shrinkage of the fizz zone. It therefore seems rea- 
wnable to use the value of C ¢,/X corresponding to y;’ in Fig. 

That assumption is probably fairly accurate just above 
the threshold velocity but may be questionable at much 
higher velocity. 


Relation Between yi’ and U 


The dependence of the erosion ratio m/m, upon the dis- 
tance y’ at which the temperature 7'/’ is achieved, for a given 
pressure, will be of little help without a relation between that 
distance 7,’ and the mean flow velocity U parallel to the 
burning surface. 

The present study will be restricted to circular channels of 
ndius r. The flow pattern across the channel can be ex- 
ected to be somewhat similar to the flow inside a tube where 
the thickness of the laminar sublayer close to the wall de- 
reases when U increases. 

locity for which the Reynolds number of the flow is well 
hove the one corresponding to the transition, it seems rea- 
wnable to assume that, despite the presence of the flame, the 
fow pattern quantitatively follows the well-known relation 
giving the mean flow velocity U in terms of the so-called 

friction velocity u, (rough pipe) i 


(175 + 5.75 logio 
v 


where » = p/p is the cinematic viscosity. 
That approximation has also been used by Corner (3) and 
hay be explained by the fact that, near the threshold velocity 
. the velocity perpendicular to the burning surface repre- 
wits only 2 or 3 per cent of the mean flow velocity U. On 
e other hand, in the laminar sublayer, the flow pattern is 
aitirely determined by the viscous stresses. Therefore, neg- 
keting the velocity component perpendicular to the burning 
uirface, it can be written, by analogy with the classical flow 
ina tube 


: y* = [6] 
v 


here y* is the so-called friction distance parameter and y the 
stance from the burning surface. Equation [6] defines the 
iminar sublayer up to a critical value of the friction distance 
jarameter y.-* which is characteristic of the problem and 
presents the point where incipient turbulence appears. 
Above that limit, an increase in heat transfer is observed due 
toturbulence. 

Anormal approach will therefore be to work out completely 
the problem again with a variable thermal conductivity at 
ach point y, taking into account the turbulence level re- 
ated to the value of the mean flow velocity U. Such a 
uethod would lead to very elaborate calculations, as has 
en recognized by Geckler. We have tried to simplify the 
jtoblem by assuming that erosive burning is only related 
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Moreover, above the threshold 


Mo 


to the flow behavior inside the fizz zone. It is obvious that 
such an approach is valid only when, without erosion, the 
heat flux reaching the fizz zone from the flame zone is negli- 
gible (approximately below 80 kg/cm?). 

Moreover, we have supposed, as a working assumption, 
that inside the fizz zone, the temperature 7,’ is obtained at 
the limit of the laminar sublayer defined by the critical value 
Yer* of the friction distance parameter. In the neighborhood 
of the threshold velocity, the increase of the temperature level 
of the fizz zone is fairly small (T;’ ~ 7,), and turbulence is 
not yet fully developed at the beginning of the preparation 
zone. The influence of the turbulence is not only to increase 
the thermal conductivity by transporting hot gases toward the 
burning surface. Hitherto unreacted products are also 
transported back in the process in such a manner that they 
will react closer to the surface. Our working assumption 
seems to be a fair approximation, close to the threshold ve- 
locity. 

In the theory of flow in tubes, the critical value of the fric- 
tion distance parameter is about 5, but in the present case a 
lower value can be expected due to the flow perpendicular to 
the burning surface and to the strong temperature gradient. 
Therefore Equation [6] can be written 


where y;’ is the distance at which we have assumed that the 
temperature 7’;’ is achieved. 

The three unknowns to be determined are U, uz and y-r*, 
and it is therefore necessary, in addition to Equations [5, 7], 
to find a third condition. The use of the experimental value 
of U,, corresponding to y;’ = 4 gives that third condition and 
allows the computation of y.,*. As expected, in the follow- 
ing example, we found y.,* = 2.85 < 5. 


Computation of the Erosive Constant k,, 


We are now in a position to compute the erosive constant 
k, defined, for U greater than U.., by 


Consistent with the experimental data, the value of k, will 
be slightly dependent upon the value of the term (U — U,,) 
with its maximum when (U — U,,) is approaching zero. 

The sequence of the computation is the following: 

1 For different values of y:’ < y:, compute from Equa- 
tions [3 and 4] the corresponding values of the mass burning 
ratio m/mp, using the value of C ¢,/X corresponding to y:’. 

2 Fora given value of U,, known from experiments, com- 
pute from Equations [6 and 7] for y:’ = y: the corresponding 
value of u, and of the critical distance parameter y;. 

3 Use that critical value in the same equation for the 
computation of the u, and of the U corresponding to various 
y;’ smaller than y;. 

4 Finally from Equation [8] compute the 
erosion coefficient k.. 


values of the 


Numerical Example 


From Rice and Ginell the characteristics of the HES. 4016 
propellant can be summarized 


i 

T;'’ =1400K 

E/R = 8000 K cae 

v = 7.05/p cm?/see (p is the pressure given in atm) 


The results of these computations in the fizz zone are given 


in the first four columns of Table 1. 


We have also indi-. 


cated the corresponding values of Cé,/X calculated 
601. 


t 
Yor = a 
& 
= 
| 


Equation [2] from which Fig. 2 has been established. 
_ ealeulations are quite sensitive to the way these three points 
have been connected graphically.) The values of Cé,/X are 
in fairly good accord with the value of é)/X = 1140 em sec/ 


= gm estimated separately. See Table 1. 
Table 1 
Pp T's Mo "1 C 
atm K gm/cm? sec cm cm sec/gm 
42 770 0.42 0.0063 455 
22 840 1.01 0.00114 1400 
82 910 2.10 0.00031 1850 


Let us now compute the erosion coefficients k,, vs. U for the 
; "pressure of 22 atm assuming, in accordance with experimental 
results mentioned by Wimpress (4) and Green (5) for J.P.N. 
Ballistite, that U,, = 180 m/sec. (J.P.N.is a propellant simi- 
lar to HES. 4016.) We took r = 1.45 cm (where r = radius 
of the port area). See Table 2. 


Table 2 


Cap Ts 
em cmsec/gm K 


m/M Ur U ky 
em/sec m/sec sec/m 


0.0010 1468 845.3 1.062 912 205.3 0.00245 
0.0008 1570 855.5 1.188 1140 263.5 0.00225 
0.0006 1680 870.0 1.390 1520 362.0 0.00214 


a1 The curve m/m, vs. U is given in Fig. 3. It corresponds 
fairly well with the experimental curve given by Wimpress 
and Green for J.P.N. (5). These authors have found k, = 
0.00219 sec/m just above the threshold velocity. The correla- 
tion given by Geckler (2) and based on Green’s data, gives 
k,, = 0.00163 sec/m. 

At greater pressure (82 atm) we have found, using the same 
value of the threshold velocity, a slightly greater erosion 
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Fig. 3 Erosion ratio r/r) vs. mean flow velocity U 
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coefficient (about 0.003 sec/m). 


experimental and theoretical results suggests also an infly. 
ence of the pressure and possibly of the radius of the port 
area upon the value of threshold velocity U1, but the problem 
is not yet sufficiently advanced to allow definite conclusions 


te be drawn. 


~ Despite numerous drastic assumptions, the quan titative 
results obtained are in fairly good agreement with the ex. 
perimental data. Although some of the assumptions may ly 
questionable, the physical picture of the erosive proc:ss tak- 
ing into account the various zones of the flame repr: sents 
frontal attack upon a highly complex problem. The proposed 
method provides a physical explanation of the comp):tations 
made by Corner (3) who has assumed that the heat exchanges 
inside the flame zone can be characterized by the p::;ticular 
value of the thermal conductivity at the end of the fizz zone, 
(See also Geckler (2).) However Corner does not t:ke into 
account the possible existence of the threshold velocity and 
assumes that the rate-controlling process of the conibustion 
is in the gas phase. 

Lenoir and Robillard (6) have also studied the )roblem 
of erosion neglecting the threshold velocity but assuming 
that the rate-centrolling process of the combustion is the 
solid phase decomposition. According to them, the erosion 
is due solely to the convective heat transfer between the 
burning surface and the mean flow, and is governed by the 
same equations as a normal heat transfer problem includ- 
ing transpiration from the walls (as a matter of fact our 
method may be somewhat improved by the use of the bound- 
ary layer equations with transpiration). In our opinion 
such an approach is oversimplified since the temperature 
gradients, close to the burning surface, are mainly imposed 
by the chemical reactions inside the flame, the flow pattern 
having only a second order influence. Moreover their good 
numerical agreement has been obtained by matching the equa- 
tions with the experimental data in order to determine the 
values of the constants involved. 

A last attempt to study the erosive combustion has been 
made by Tavernier and Boisson (7). Their theoretical 
treatment assumes that the rate-controlling process of the 
combustion is in the gas phase and is therefore in opposition 
to the theory generally accepted. 

Meanwhile they have found an erosion coefficient about 
-_ twice greater than the one obtained by Green with a similar 
propellant. In our opinion, that interesting discrepancy can 
_ be explained by the experimental technique used. While 
_ Green measured the erosion of the main propellant grain in 
_ the chamber, Tavernier and Boisson used two grains, the main 


_ grain supplying the flow and the test sample on which erosion 


is measured being situated at the entrance of the nozzle and 
therefore being subjected to a higher turbulence. 
By assuming that the erosion is normally due to a process 
entirely localized in the fizz zone, our study enables one to 
- understand why the presence of instabilities can determine 
_ burning rates much higher than those that can normally be 
expected due to erosion alone. Indeed, with an abnormal 
re degree of turbulence, hot products from the flame zone may 
reach the fizz zone, correspondingly increasing the burning 
rate. That explanation is supported by the fact that, in a 
liquid propellant rocket, the temperature of the walls is 
_ strongly increased by instabilities (see (8)). 

The present theory is limited to double base propellants 
burning at fairly low pressures because we have been able 
only in that case to carry out quantitative calculations. 
Meanwhile, the fact that, after the disappearance of the dark 
layer of the flame, the pressure dependent burning rate law 
remains unchanged up to a much higher pressure suggests 


“ma a that the heat transfer toward the fizz zone remains small. 
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Such an increase seems to be 
in accordance with the data given by Green. The study of the 
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Therefore our method may possibly still hold above 80 atm. 

On the other hand, proposed theories for composite pro- 
pellant burning also allow for the decomposition of the flame 
into separate zones, and possibly the treatment may also be 
extended in that case if suitable characteristics of the different 
zones can be obtained. 
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Rhodes and Bloxsom 


Air at temperatures up to 9000 K and pressures to 40,000 
psi has been expanded in a deLaval supersonic nozzle. 
Experimental measurements of the air pressure indicated 
that the air was substantially in equilibrium during the 
entire isentropic expansion. After the supersonic jet 
was expanded to one atmosphere ambient pressure, flow 
patterns around various models were studied as to self 
luminosity of the flow, equilibrium state of the air, drag 
coefficients, and laminar flow stagnation heat transfer 
coefficients. The experimental drag and heat transfer 
coefficients agreed with commonly accepted theories in the 
field. Heat power inputs up to 363,000 w/cm? in the nozzle 
throat produced varying amounts of erosion of the throat 
material. This erosion is presented with the correspond- 
ing wall vaporization theory. 


Nomenclature 


a = speed of sound, cm/sec 

« = acceleration, cm/sec? 

A = surface area, cm? 

C, = gas specific heat at constant volume 

Cy = drag coefficient 

Cy = specific heat of solid wall, joules/gm K 

’, = specific heat of liquid wall material, joules/gm K 

d* = diameter of sonic flow throat, cm 

dy~= shock detachment distance, em 

F = force, dynes ‘ 

4H = gas enthalpy difference between stagnation and wall 
conditions, joules/gm 

Ky = Stanton number 

m = mass, gm 

Mp. = mass of wall material, gm 

= gas viscosity, poise * 

Q = heat, joules oi 

Qm = latent heat of melting of wall, joules/gm K 
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Supersonic Aerodynamic Experiments Using Very High 
Temperature Air Wind Tunnels 


D. E. BLOXSOM Jr.’ 


Canoga Park, Calif. 
agli 


« 
= latent heat of vaporiza 


Qr 

ge = equilibrium heat transfer coefficient, w/cm? 

ge = average equilibrium heat transfer coefficient, w/em? 
dep = black body heat transfer coefficient, w/cm? 

dsr = forced convection heat transfer coefficient, w/cm? 
Gne = nonequilibrium heat transfer coefficient, w/cm? 
R = sphere radius, cm ; 
R, = gas constant 

Re = Reynolds number based on throat diameter 

p = density, gm/em* 

pw = wall density, gm/em* 

S = displacement, em q 

T = temperature, K 

T = initial wall temperature, K 

Tm = melting temperature of wall, K 

T, = vaporization temperature of wall, K 

t = time, sec 

u = velocity, cm/sec 

y = erosion distance, cm 


subscript w refers to wall conditions 
superscript * refers to sonic flow conditions 


Introduction 


ERONAUTICAL and astronautical engineers at present 

need to measure aerodynamic parameters at elevated 
temperatures. These air temperatures are produced by high 
velocity flight in thé earth’s atmosphere. This high stagna- 
tion air temperature is plotted against flight Mach number in 
the isothermal stratosphere in Fig. 1. This figure was com- 
puted with the aid of the real gas air equation of state (1).* 


countered at earth satellite velocities. In addition to the high 
temperature requirement, high stagnation pressure is re- 
quired for wind tunnel operation due to the high total pressure 
drop encountered in very high Mach number steady flow isen- 
tropic air expansions. In order to meet these requirements, a 
method of heating air by means of condenser electrical sparks, 
was developed by the United States Air Force at the Arnold! 
Engineering Development Center (AEDC) at Tullahoma,’ 


; 7 ig. 1 shows that temperatures up to 10,000 K will be en- 


3 Numbers in parentheses indicate References at end of paper. 
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Fig. 2 Experimental static air jet pressures 


Fig. 3. Direct luminosity photograph of 20 deg steel cone in 
Mach number 4 flow 
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Tenn. This technique is described in (2) which includes the 
technique of operation and impurity level experimentally 
obtained. Bomb 2, described in (2), enabled air heating to 
5000 psi and 7200 K. Bomb 3, also described in (2), enabled 
air heating to 40,000 psi and 9000 K. The nozzles used for 
supersonic flow production opened out of the bombs. 


Supersonic Nozzle Air Expansion 


These experiments were performed in bomb 2 (2). ; 

The sudden rise of air pressure in the stagnation chamber 
of the wind tunnel blows a thin diaphragm at the wind tunnel 
throat. The air then rushes through the nozzle and starts the 
wind tunnel. The experimental tunnel starting times \ ith a 
back pressure of 0.97 atmosphere and an area ratio of 16.2 were 
measured with a 14,000 frames/sec Fastax camera. The start 
of the heating pulse could be easily seen on the film, ari! the 
fully started flow appeared about 130 microsec later. The 
testing described below occurred in the next 500 mic:osee 
following the tunnel opening. The stagnation pre-sure 
dropped less than 10 per cent of maximum pressure in the test- 
ing time (2). 

The static jet pressure was measured by inspection o! the 
first few frames showing the jet. The area ratio of the nozzle 
is chosen so that the static jet pressure is close to atmosplieric. 
The jet ambient pressure is usually slightly over the lalora- 
tory pressure (0.97 atmosphere), then matched, then slightly 
under the laboratory pressure as the stagnation pressure «!rops 
in bomb 2 as a function of time. The bomb pressure de- 
creases due to heat transfer and mass flow out of the tunnel 
throat. As little as 2 per cent change in static pressure c:n be 
detected in the jet configuration of this hot air near the labora- 
tory pressure matching point. Jet details were easily seen as 
changes in the self luminosity due to jet expansion waves and 
shock waves. The outline of the jet was not used in these ob- 
servations, which were made on the jet interior. The data on 
this static pressure are given in Fig. 2. If the flow is com- 
pletely in equilibrium, an isentropic real gas steady flow ex- 
pansion can be made, using a Mollier chart for air (3). If the 
gas is completely out of equilibrium, except for the translational 
and rotational degrees of freedom, an isentropic flow expan- 
sion using y = 1.4 can be made. Fig. 2 shows that the tunnel 
flow is almost in equilibrium. 

Assuming the flow to be in equilibrium, the tunnel calibra- 
tion (+15 per cent) can be completed. The Mach number is 
around 4 for these experiments. 


Flow Visualization by Gas Luminosity 


The heated gas is quite luminous. There are two types of 
luminosity. One, very brilliant, always occurs in the first 500 
microsec of flow from the tunnel, even when the ambient tem- 
perature is less than 1000 K and no luminosity is expected. 
The second type of luminosity is fairly steady in time, to 20 
millisee. For stagnation temperatures of 3000 K, this lumin- 
osity first appears in the boundary layer of blunt objects, such 


Fig. 4 Direct luminosity photograph of 10 deg wedge in Mach 
number 4 flow (14,000 frames/sec). Air stagnation tempera- 
ture: 3500 K. Jet ambient pressure: 1.0 atm 
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Fig. 5 Direct luminosity photograph of }-in. steel sphere in Mach number 4 flow (14,000 frames/sec) 


as spheres. As the stagnation temperature increases, this 
juminosity then appears in the shock heated gas on blunt 
object- (3350 K) and then appears finally in the shock heated 
gas of slender bodies (to 7200 K). This is probably the true 
gas luminosity. The first type of transient luminosity is 
proba! ly due to the dying down of the active nitrogen meta- 
stable atomic states excited by the heating spark. Reference 
(4) shows that the active nitrogen will not affect aerodynamic 
supersonic flow parameters significantly. Fig. 3 is a remark- 
ably clear time integrated photograph taken in the nitrogen 
aftergiow. Shock waves and expansion waves can be clearly 
seen around the 20 deg cone. 

Shock waves can also be easily seen in Fig. 4, which shows 
the development of a hot boundary layer on the upper surface 
of the 10 deg wedge. At low temperatures, this luminosity 
layer looks like a normal boundary layer. At high tempera- 
tures, however, the luminous layer gets very thick. It is not 
obvious why this luminous layer is so large at high tempera- 
tures. The luminous layer decreases in size aft of the wedge 
due to expansion waves in the jet. 


Equilibrium Studies of Laminar Sphere Flow 


When }-in. steel spheres were placed in the flow, it was 
possible to determine shock detachment distances by means 
of the 14,000 frames/sec Fastax camera. This camera was 
aimed perpendicular to the jet as is shown in Fig. 8. An f-2 
lens permitted reasonable exposures by means of the self 
luminosity of the flow. Fig. 5 shows such a Fastax camera 
frame of a sphere at 3550 K stagnation temperature. The 
shock wave, boundary layer, boundary layer separation and 
shock detachment distance may be clearly seen. The flow 
separates on all spheres used in this report at an angle of about 
73 deg measured from the stagnation point. This flow separa- 
tion is characteristic of laminar flow in the boundary layer. 
Recovered spheres were melted only to this angle, the rest of 
the sphere being intact. These shock detachment distances 
are plotted in Fig. 6. That the shock detachment distance 
would be sensitive to equilibrium state was first suggested by 
Schwartz and Eckerman (5). Their formula for Mach num- 
ber 4 for the shock detachment distance is given also in (5), 
page 170. The values for C,/R, = 3/2 (monatomic gas) and 
C,/R, = 5/2 are plotted in Fig. 6. Two other curves are 
plotted i in Fig.6. One, entitled ‘Tunnel Flow in Equilibrium, 
Sphere Flow Out of Equilibrium” uses the C,/R, characteristic 
of the equilibrium ambient tunnel flow; whereas the curve en- 
titled ‘Tunnel Flow and Sphere Flow in Equilibrium” uses 
the value of C,/R, characteristic of complete equilibrium 
sphere flow conditions. The experimental data of Fig. 6 in- 
dicate that only translational and rotational degrees of free- 
dom are activated in the sphere flow. 

The air equilibrium data discussed in this paper do not dis- 
agree with Logan’s published data (6). Logan has collected 
theoretical and experimental relaxation times for dissociation 
at | amagat* density for air, oxygen and nitrogen, using shock 
tube techniques. AVCO Research Laboratory has also re- 
ported on air dissociation rates in this temperature region, 
using radiation shock tube techniques (15). These data are 


‘| amagat is the density at 300 K and 1 atmosphere pressure. 
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summarized in Fig. 7. The stagnation gas is in equilibrium 
by the time the pressure is measured (in 100 microsec), since 
the characteristic gas holding time referenced to 1 amagat 
density is over 1 millisec. ‘The tunnel gas is marginal in equi- 
librium, as it has about 16 microsee to achieve equilibrium. 
The sphere flow is completely out of equilibrium, as it has a 
characteristic time of only 10-8 see to achieve equilibrium. 
Hence, the experimental equilibrium data of this report are not 
in conflict with previous theoretical and experimental determi- 
nations of the relaxation time of air. Inreduction to 1 amagat 
density, the relaxation times for air were assumed to scale ac- 
cording to the following law 


ta 
Data is currently available (15) which indicates the air dis- 
sociation density dependence is tap~-™. a density cor- 
rection of Fig. 7 is small on the sphere data, however. 


Drag Coefficients of Spheres in Laminar Flow 


In order to demonstrate the feasibility of measuring forces 
in a blowdown tunnel of millisee time duration, a steel sphere 
+ in. in diameter was mounted ballistically in the tunnel, as 
shown in Fig. 8. A small wire was used to suspend the sphere 
in the center of the wind tunnel. While the wind tunnel was 
running, the sphere moved from position 1 to position 2 in 500 
microsec, due to the drag forces on the sphere. Gravitational 
force during this time was negligible compared with the drag 
force. 

By matching seven displacement points, obtained by the 
14,000 frames,sec Fastax camera, with a curve of the type 


S = ut + dat? 


is is possible to determine the accelerations to 3 per cent. 
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Fig.6 Experimental shock detachment distances on }-in. steel 
sphere 
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Experimental relaxation times as 


Therefore, the forces are determined to the same accuracy 


per unit time 


suitable hypersonic boundary layer corrections. 


Laminar Flow Heat Transfer Coefficients 


from work on high temperature tunnels is realistic heat trans- found in (12). 


AIR STAGNATION TEMPERATURE, °K —> 


Knowing this mass we compute 


Q = mo[Cu(Tm — Tw) + gm] joules 


Fig. 9 Experimental drag coefficients of spheres 


fer coefficients, defined as average heat transferred per unit 


One of tl awe ase , 5 This early crude technique was intended to indicate feasibility 
ne of the most important types of information to come of heat transfer measurement only; more refined techniques are 
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and with the help of the tunnel calibration, a drag coefficient 
may be obtained based on the projected sphere area ae a - 
= = 2/ 
Our method of measuring the heat transfer coefficients ex- 

These drag coefficients, plotted in Fig. 9, are only accurate to perimentally involves the melting of models. This technique 
15 per cent due to uncertainties in the tunnel calibration. has been used successfully at 19 w/cem? (9). We use this 
Fair agreement is obtained with the values of drag coefficients technique to measure heat transfer coefficients in this tunnel 
for hypervelocity pellets, in this velocity and Reynolds num- by the melting of solid steel models. The steel cones exposed 
ber range, as given by Allen (7). to the hot flow eroded to a surface of minimum heat transfer 

Associated with the curve matching technique was an for a 20 deg cone. This technique can be used to find such 
initial velocity about one quarter to one half the final velocity minimum heat transfer surfaces for any configuration. For 
as the sphere left the test section. This means that the start- sealing purposes, the mode of heat transfer, laminar skin fric- 
ing processes put a force about five to ten times the steady tion, turbulent skin friction or radiation must be duplicated. 
drag force on the sphere for 50 microsec. This ballistic tech- When a steel model of 0.003-in. radius hemispherical tip 
nique looks promising, even at very low drag forces, as nylon is placed in the flow, the model acts as a heat sink for about 
models may be used with correspondingly higher accelera- 500 microsec. Then suddenly the tip erodes, spreading high 
tions. luminosity in the flow around the tip. Melted and vaporized 

The dynamic force }pu2A changed about 48 per cent material can be seen to flow away from the tip. This melting 
as the sphere moved in the source flow test section. For greater time can be read to 70 microsec on the Fastax pictures. The 
accuracy an axisymmetric Foelsch (8) type nozzle should be amount of melted material can also be seen by means of the 
used. This nozzle produces a parallel and uniform jet with model’s change in shape in the frame that the melting occurs. 
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About 5 micrograms were melted in this time, requiring about 
4 millijoules of heat. Since the surface area exposed to the 
jaminar flow is known to be 365 X 10~* cm?, an average heat 
transfer coefficient may be computed. This average heat 
transfer coefficient is plotted in Fig. 10. 

Lees has published a theory on laminar heat transfer over 
spheres at hypersonic flight speeds (10). The stagnation 
point heat transfer for equilibrium flow is given in (10). 
The average sphere equilibrium heat transfer for laminar flow 
can be computed from the heat transfer rate distribution as 


de = 0.4 de 


Using the tunnel calibration and sphere radius, the average — 


heat transfer equilibrium coefficients can be computed. 
These computed values are plotted in Fig. 10 as ‘‘Lees’ Equi- 
librium Limit.” 

According to the rate processes given by Lees, the maximum 
difference for a given flight speed in surface heat transfer 
rate between the equilibrium and completely out of equilib- 
rium flow with a catalytic wall is given by 


Gne = 1.3 de 


This is also plotted in Fig. 10 as “Lees’ Nonequilibrium 
Limit.’ It will be seen that there is very good agreement be- 
twee: the theory and the experimental data. The maximum 
possible black body radiation is also plotted to show that this 
form of heat transfer is not important in this case. 

Observation of the models in the flow indicate that the flow 
is separating about 85 deg from the stagnation point and 
hence is laminar. 

Further work on the theory of stagnation point heat trans- 
fer in dissociated air has been done by Fay and Riddell (11). 
This work indicates that Lees’ values are slightly low. Further 
experimental work in this field has been done by Rose and 
Stark (12). The heat transfer coefficients of Fig. 10 agreed 
with the published data in (12). 


Wind Tunnel Throat Erosion 


If bomb 3 (2) is converted into a wind tunnel, heat power 
rates over 100,000 w/cm? are obtained in the sonic flow 
throat. A picture of this bomb running as a Mach number 5 
wind tunnel is given as Fig. 11. The throat is cooled by the 
vaporization of the throat material. The total energy input 
into the throat is given by 


Q = Ate = 
Ayp(Cu(Tm — Tw) + — Tm) + + Qr] joules 


This equation assumes one-dimensional heat transfer and 
neglects the unsteady heat conduction through the unmelted 
wall. If the values from Table 1 are used in the above equa- 
tion, it will be seen that the latent heat of vaporization is 
much greater than the other terms of the wall heat equation 
combined. 

If oxidation occurs at the point of vaporization, the oxida- 


-— 
Table 1 Properties of wall materials 
Ma- Car- Tung- Molyb- Units 
terial bon sten denum 
K 0.24 1.99 1.45 joules/K-cm sec 
Cy 1.90 0.163 0.314 joules/gm K 
1,» 300 300 300 K 
Tm \Sublimes 3670 2920 K 
T, 1 at 3950 6227 5103 K 
Q, 59,800 4605 6823 joules/gm 
pQ, 134,600 87 , 360 69,400 joules/em’ 
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Fig. 10 Experimental heat transfer coefficients 


Fig. 11 Mach number 5 flow. Air stagnation temperature: 
9000 K. Air stagnation pressure: 40,000 psi. Jet ambient 
pressure: 1 atm 


tion process will produce heat at a rate in excess of the cooling 
by vaporization, since all these oxidation processes are exo- 
thermic. This would make the vaporization heat sink capa- 
bilities of various substances useless. However, we will show 
below that oxidation close to the surface does not occur in the 
boundary layer in this particular wind tunnel throat. This 
oxidation undoubtedly occurs in wakes or supersonic portions 
of a wind tunnel, whose throat is eroding in this manner (13). 
The steady state heat vaporization condition will only be 
reached after a certain amount of time—the time required to 
heat the surface to the melting point (14). This time is 


Using Table 1, we obtain for = 363,000 w/em? 
carbon t = 0.098 X 1073 sec _ 
tungsten t = 0.49 X 1073 sec 
molybdenum t = 0.24 X 1073 sec 


The characteristic testing time in the wind tunnel is 0.00227 
sec. 

After the time required to melt the wall, the wall is cooled 
by the vaporization of its surface, and the heat balance equa- 
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Fig. 12 Sonic throat erosion by means of high temperature air 
flow 


tion may be approximated for our heat power input regime 
ypQy = joules/em? 

The following conditions are obtained in the sonic wind 
tunnel throat of diameter d* (as computed from real air isen- 
tropic expansion from measured stagnation conditions) 


= 3.61 X 107? poise 

d* = 0.0762 cm 

T'stag.* = 9000 K 

p* = 0.0501 gm/cm* 

a* = 2.13 X 10° cm/sec 

(AH)* = 19,000 joules/gm 

Re* = 5.2 X 10° (based on tunnel throat diameter) 


This leads, by Reynolds’ analogy for turbulent flow, to a 
Stanton number of 0.0016 = Ky, (Prandtl number 1), for 
smooth walled pipes. This is probably a minimum Stanton 
number for this kind of flow. Three small corrections which 
all raise the heat transfer rate are: 

1 As the gas is black, the heat input rate per unit area for 


black body radiation is 
dep = oT* = 5.67 X 10712(9000)* = 37,600 w/cm? 


2 Due to the fact that we only have three diameters of 
pipe flow, the inlet heat transfer rate is about 20 per cent 
higher over these three diameters. 


3 Oxidation can change the heat input rate, raising it by 
as much as a factor of three for some materials; however, it 


does not occur here, as will be shown. 
The expected heat input rate for the minimum Stanton 
number is 


Gsr = Kyp*u*(AH)* = 3.25 X 10° w/cm? 
Adding the black body radiation 
Ge = Gps + Gsr = 363,000 w/cm? 


This is taken as a minimum total. 

The bomb has been observed to run for at least 0.020 sec in 
supersonic flow. However, the primary heat transfer occurs 
in the first 5 millisec (observed by means of Fastax camera). 
If we assume a characteristic running time of 0.00227 sec, we 
get a constant expected energy input to all tested throat 


materials of 823 joules/cm? = get. 
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The quantity pQ, is governed by the material and is called 
the volumetric heat of vaporizatiou. For our materials, this 
quantity is tabulated in Tables 1 and 2. 

Expected values of y as radius enlargement of initial radiys 
= 0.038 cm throat are found in Fig. 12. 

Graphite, tungsten carbide, molybdenum, lava (AISiO;), 
steel and beryllium copper were run in the sonic throat, 
These results are plotted in Fig. 12 along with the theory as 
worked out in the first part of this paper. Results are plotted 
as a function of melting temperature. The clear superiority 
of graphite over all other materials at 40,000 psi stagnation 
pressure seems to indicate that it is the preferred materi:! for 
wind tunnel throat under these conditions. 

The lava broke in pressure fracture. Improperly supported 
carbon also broke in pressure fracture. The other matvrials 
withstood the pressure but not the high heat power rate. If 
the melting point of the material gets much below 2500 K, 
the simple theory proposed above fails, due to the fact that 
the wall softens and the pressure blows the molten throat 
out in drops. This problem also occurs in electrical arv dis- 
charges at the electrodes, and it is recommended that such 
materials be avoided in this heat power range. 

Since the experimental erosion of the best matvrials 
(graphite, tungsten carbide and molybdenum) followe: the 
simple vaporization theory without oxidation, this is evidence 
that no additional heat input due to oxidation occurred in the 
sonic throat. If the wall had oxidized, the experimental! y de- 
termined erosion depths would have varied with the differing 
heats of oxidation and not with the heats of vaporization as 


was determined. q 
Conclusions 


A series of typical aerodynamic experiments have been pre- 
sented to demonstrate that data of interest to aeronautical 
engineers can be obtained in wind tunnels operating in millisec 
intervals. 

Supersonic nozzles were constructed and tested by means 
of jet pressure observations, proving that the isentropic 
equilibrium flow computations are valid for such high tem- 
perature, high pressure nozzle expansions. The flow about 
models can be observed in the self luminosity of the air. The 
flow equilibrium studies that have been carried out agree with 
other experimental shock tube results. Drag and heat trans- 
fer coefficients have also been made. For spheres, these re- 
sults agree with commonly accepted present theories in the 
field. 

Throat erosion is a very serious problem in this type of 
wind tunnel, due to the high temperatures and high heat flux 
rates. An experimental study is presented which demon- 
strates theoretically and experimentally that tungsten and 
carbon throats are to be preferred if wall vaporization is the 
throat cooling technique. 
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Table 2 Heat of vaporization and melting temperature 


Material Heat of Vaporization Melting — 
K 


joules/cm* 
Graphite 134,600 3950 7 
Tungsten 87,360 3670 
Molybdenum 69, 400 29200 
Iron 56,850 1835 
Copper 


JET PROPULSION 


~ 
Tl 
en 
- 

na 

- 
du 
vo 

of 
Ju 

Fr 

: 

of 

cor 

SIZ 

fro 
4 per 
pre 
cen 

pre 
the 

Hy 

I 

H 
m 
i 
; 
— 
po heer 
pres 

% 


Without these facilities the work would not have been done. 
The author would especially like to acknowledge the assistance 
given him by B. V. Rhodes and R. Smelt, whose continual 
encouragement and ideas were most valuable. 


References 


1 Hilsenrath, J. and Beckett, C. W., “Tables of Thermody- 
namic Properties of Argon-Free Air to 15, 000°K,’’ AEDC-TN- 
56-12, Sept. 1956. 

2 Bloxsom, D. E., Jr., 
duce High Temperature, High Pressure Air,’’ 
vol, 28, no. 9, ry 1958. 

3 F ‘eldman, $ “Hypersonic Gas Dynamic Charts for Equi- 
librium Air,’’ AV CO Research Laboratory, Jan. 1957. 

4 Benson, J. M., “Measurements of the Physical Properties 
of Active Nitrogen,’’ Journal of Applied Physics, vol. 23, no. 7, 
July 1952, pp. 757-763. 

5 Schwartz, R. N. and Eckerman, J., “Shock Location in 
Front of a Sphere as a Measure of Real Gas Effects,’’ Journal of 
Applied Physics, vol. 27, no, 2, Feb. 1956, pp. 169-174. 

6 Logan, J. G., “Relaxation Phenomena in Hypersonic Aero- 
dynamics,’’ IAS Pre sprint 728, Jan. 1957. 

7 Allen, W. A., Rinehart, J. S. and White W. C., 


“Use of Capacitor Discharges to Pro- 
JET PROPULSION, 


“Phe- 


An electrical spark gas discharge technique, using con- 
densers, is described, which has been used to heat air to 
9000 K and 40,000 psi at close to 100 per cent energy transfer 
efficiency. The heat losses from the air at near stagnation 
conditions have been measured. These heat losses restrict 
the holding time of the hot air to milliseconds in laboratory 
size apparatus. The impurities introduced into the air 
from the heating and electrode erosion have been ex- 
perimentally measured, and these measurements are also 
presented. The total of these impurities is about | to 5 per 
cent of the total heated air mass. Other experiments are 


presented to show that the airisinequilibrium. Various 


theoretical explanations are given for the experiments. 2] 


Nomenclature ry 
4 = surface area, cm? 
D bomb diameter, em 
H, = specific enthalpy of air at stagnation temperature condi- 
ditions, joules/gm 
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of Stagnation 
Aero- 


Temperature, High Pressure Air 


dsr = forced convection heat transfer coefficient, w/cm? — 
Qvy = heat of vaporization, joules/gm 

p = air density, gm/cm$ 

Re = Reynolds number of air a> 

= air temperature, K 

t = time, sec 

u = air velocity, cm/sec 


AV = voltage change, at electrodes ~30 v 


Introduction 


ECEDING work in high pressure, high temperature 
electric condenser discharges has been done by J. A. 
sarly 1920’s (1, 2). This work involved the 
sudden release of electrical energy stored originally in capaci- 
tors. The air heating technique described in this paper is 
similar to Anderson’s technique except that the electrical 
explosion is confined in a closed bomb, and the data are taken 
by means of fast pressure probes. These modifications of 
Anderson’s technique extend the time scale of the experiment 
to millisecond time intervals and make possible gas dynamic 
and heat transfer studies at these elevated air temperatures. 
Roth (3) has done work indicating very high air electrical 
heating efficiency from smaller ignition-type sparks. 


Air Heating by Capacitor Discharge 


The spark heating of air considered here deals with a con- 
stant volume and constant density condition. No details of 
heating during spark growth are available from the instrumen- 


3 Numbers in parentheses indicate References at end of paper. 
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tation used, as the minimum time of response of the instrumen- 
tation was longer than the entire time of spark discharge. 

The air was heated in a closed chamber called a bomb. 
Bomb 1 was a glass tube 18 in. long with an inside diameter of 
in. The bomb was used to heat air to a maximum pressure 
of 300 psi. Bomb 2 was 2 in. long with a 1.6 in. inside diame- 
ter. This bomb was made of steel and permitted holding of 
air pressures to 5000 psi. Bomb 3 was also fabricated of steel, 
with a 2 in. inside diameter and was 1 in. long, which per- 
mitted holding of air pressures to 40,000 psi. 

The procedure for using these bombs was to first fill the 
bomb with high pressure air at laboratory temperature. The 
high pressure air was restrained with a diaphragm and valves. 
Then, a bank of capacitors was discharged into the high 
pressure air in the bomb. A small wire was used to initiate 
the spark between the electrodes, as the gaps were too long to 
spark over. A low resistance spark gap, which could be trig- 
gered by external means, was used to permit 3 microsec high 
voltage rise across the tease wire. The usual discharge circuit 
for bombs 2 and 3 may be found in Fig.1. The total discharge 
time of the circuit was 50 microsec. It is very important, for 
this type of electrical heating, to keep the discharge time of the 
electrical circuits small relative to the cooling time of the hot 
air. 

As the air is heated to higher and higher specific enthalpies, 
more and more of the energy is transferred to other than 
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_ Fig. 1 Schematic of electric circuit for bombs 2 and 3 


La. 
BOMB 3 
a EXPERIMENTAL POIN 
—— 
c 
> 


0 2000 4000 6000 8000 10,000 


AIR TEMPERATURE °K —> 


Fig. 2 Overall efficiency of energy transfer from condensers to 
usable air 


translational degrees of freedom. The amount of this energy 
for thermodynamic equilibrium may be found in the equation 
of state for air (4). By means of these equations of state, the 
specific internal energy, specific enthalpy, etc., may be plotted 
as a function of density, pressure and temperature. If ther- 
modynamic equilibrium is assumed, then knowledge of two 
of the three air parameters, pressure, temperature and 
density, permits the other to be computed. Since the bomb is 
closed, the density remains constant in the absence of |arge 
disturbances in the air. This density is measured by measur- 
ing the initial air pressure in the bombs at 300 K. Gas density 
is measured here in amagat units, the density of a gas at one 
atmosphere pressure and 300 K temperature—for air, one 
amagat is 1.17 K 10-3 gm/cm’. 

The pressure rise in the air due to spark heating was !.eas- 
ured by means of a pressure pickup in the bomb wall. This 
was an Electro Products Laboratory, Inc. Model 37(0-C 
Pressograph capacitance pressure pickup with a critically 
damped rise time less than 100 microsec. The high current of 
the spark discharge produced a signal as large as the full pres- 
sure signal, but this inductive pickup was confined to the first 
50 microsec of the pressure trace, i.e., the first half of the 
pickup’s rise time. The pressure pickup was calibrated 
statically by means of registering oscilloscope deflections pro- 
duced by dead weight pressure on the pickup. This calibra- 
tion was made before and after a series of runs with a given 
pressure pickup. The calibrations agreed to within 5 per cent. 
The oscilloscope used in these experiments was the Tektronix 
513D. 

Using the equation of state for air from (4), equilibrium air 
temperatures were computed from the pressure data. Ex- 
periments involving the sudden cooling of the air (10° K ‘sec) 
in a short isentropic nozzle discharging into a cold air mass 
indicated that the air is in equilibrium in bomb 2 in a time less 
than 1 millisee. After 0.1 millisec, the pressure in the bombs 
was slowly decreasing due to heat transfer only. Finally, 
shock tube measurements of air equilibrium (5) indicate that 
the air is in equilibrium in essentially the time required for the 
pressure pickup to measure the pressure rise in the spark 
heated air. A comparison of these data is found in Fig. 7 of 
(6). 

Using the equilibrium air internal energy indicated by the 
pressure and density measurements, an overall efficiency of 
energy transfer from the condensers to the air can be ob- 
tained. These data are found in Fig. 2. 

It will be seen from Fig. 2 that overall efficiencies around 
50 per cent of the initial capacitor energy are obtained. 
Voltage and current measurements of the electrical power de- 
livered to the air spark were made as a function of time, using 
the circuit of Fig. 1. Care was taken not to introduce spurious 
signals into the measuring circuits during discharge by induc- 
tive loops responding to the 100,000 amp current of the spark 
discharge. To prevent these inductive loops, the resistance 
of the voltage circuit must be kept high, but not too high so 
as to prohibit high frequency performance with the distributed 
capacitance of the lines and oscilloscope plates. Around 5000 
ohms are suitable for a circuit given in Fig. 1. The inductance 
of the current measuring device must be kept as low as 
possible. If an M-shaped foil loop is used as a resistor, prop- 
erly insulated and supported, the inductive component of this 
current test resistor can be made as low as 10-!! henries. 
Both the current and voltage circuits had 0.01 microsec rise 
times to measure 50 microsec pulses. Proof that spurious 
signals were not entering the circuits was obtained by dis- 
connecting the high signal lead but leaving the lead close by 
with its ground connected. Neither circuit had enough 
spurious pickup to trigger the oscilloscope, thus proving the 
spurious signals were below 5 per cent of the peak signals. 
A study of these current and voltage oscilloscope traces of the 
spark discharge indicated that almost the full condenser 
energy enters the ohmic spark resistance. Hence the lower 
efficiencies measured by the pressure pickup (Fig. 2) are the 
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Fig.3 Pressure vs. time oscillogram of air heat transfer in 

bomb 1 


result of heat losses during the actual spark formation. No 
information concerning these losses is given in this paper, as 
the actual variable density high pressure spark formation 
could not be followed by the pressure instrumentation. 


Experimental Measurement of the Heat Losses 


The pressure measurements indicated a slow drop in pres- 
sure with time. A typical example of these data is found in 
Fig. 3. The useful experimentation time in bombs 1-3 at 
high temperature is defined as the time for the pressure to 
drop to 0.90 of maximum pressure. This time is found in Fig. 
4as a function of temperature. Fig. 4 shows that at all tem- 
peratures there is about 0.5 millisee available for such hot air 
experimentation. The maximum time of electrical heating is 
0.05 millisee. 

From the drop in average air pressure with time, at con- 
stant air density, it is possible to compute a heat transfer 
coefficient for a given time interval 


These heat power rates are found in Fig. 5 for bomb 1 heat 
losses, and in Fig. 6 for bombs 2 and 3. 

The actual heat power loss mechanism is not known at the 
present time. The black body limit to radiation heat power 
loss is given by 


dan = 6.67 X 107 [2] 


This black body limit is plotted in Figs. 5 and 6. 

Another possible mechanism of heat transfer from the hot 
air is forced convection. The heat power given by this 
mechanism is 


dsr = Kuy(Ho — Hy)pu w/em?............ [3] 


These Stanton numbers have been measured by means of 
pipe flow friction factors and are given by Reynolds’ analogy 


[4] 


These Stanton numbers are a slowly varying function of the 
Reynolds numbers for turbulent flow 


The viscosity of high temperature air is not very well 
known. Viscosity values for air used here may be found in (7). 
The friction factors may be found in (8) as a function of 
Renolds number and pipe surface roughness. 

In order to compute the forced convection heat transfer ex- 
pected in the three bombs, the average air velocities must be 
experimentally measured as a function of time. All other 
parameters of the forced convection heat transfer Equation [3] 
are known from other measurements and computations from 
the average pressure oscilloscope traces (using an equilibrium 
equation of state). 
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Fig. 5 Experimental heat transfer coefficients in bomb 1 
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Fig.6 Experimental heat transfer coefficients in bombs 2 and 3 


A measurement of the air velocities may be obtained from 
the actual pressure fluctuations measured at one end of bomb 
1. A typical oscilloscope pressure trace is given in Fig. 3. 
The N-shaped wave imposed on the steady average pressure 
permits the Mach number of the near sonic shock to be 
uniquely determined. The air velocity driving this near 
sonic shock is also uniquely determined from the Mach num- 
ber, since the speed of sound is known in the hot air. The air 
velocity-speed of sound ratio so determined as a function of 
time is found in Fig. 7. This shock wave decay is to be ex- 
pected in the bombs due to the stability of weak shocks, near 
Mach number 1, in the time scale of the experiment. It must 
be emphasized that there was no net air flow out of bomb 1, 
since the shock wave simply moved the air around in the bomb. 
Additional optical evidence of the existence of the shock wave 
is found in a series of smear camera photographs, which 
record luminosity perpendicular to the bomb 1 axis as a func- 
tion of time. A representative photograph is given in Fig. 8. 
The initial high intensity spark may be seen, and later, the 
glow of the hot gas cooling. The decaying shock wave is also 
easily seen as a sharp Z-shaped division between brighter and 
darker luminosity. Particle paths may also be seen moving in 
the tube. Gas temperature measurements from this shock 
wave velocity agree with the temperature measurements from 
the constant volume pressure rise to 10 per cent, so the gas 
temperature is believed accurate to this percentage. 
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Fig. 7 Longitudinal air shock decay as a function of time in 


bomb 1 


Figs. 5 and 6 contain estimates of the forced convection heat 
transfer, using the measured air velocities of Fig. 7 and the 
turbulent flow smooth pipe friction factors of (8). 

The experimental heat transfer coefficients of Figs. 5 and 6 
seem to follow the turbulent flow forced convection theory for 
low temperatures and the black body radiation theory for high 
temperatures. It should be emphasized again, however, 
that additional experimentation is necessary to determine the 


heat transfer mechanism. | 

Impurities can come from several sources: 

1 Condensable impurities due to vaporization of the tease 
wire and electrode surfaces. These impurities are found ex- 
perimentally as dust in the flow gas. 

2 Permanent gases eroded from the electrode shields, such 
as carbon dioxide. 

3 Permanent impurities frozen into the air ry high cool- 
ing rates, such as nitric oxide. oo, | 


Impurities 


Fig. 8 Smear camera view of longitudinal shock wave in argon luminosity, argon temperature 4150 K, argon density 1 amagat, 143 
microsec time interval between film sprocket holes on y axis (time mark, 1000 cps), 18 in. distance between film sprocket holes on x axis 


mie JET PROPULSION 


re 


5 
» 
er 
= 
: 
wey — 
= 
aed 
f 
| 
| 
3 
| 
e > « 
= 


Table 1 Dust impurity measured in bomb 2 
Initial Initial Dust mass, Dust per cent of > Expected dust 
tempera- density, gm gas mass, based per cent of gas 
ture, amagat — on mean dust mass from tease 
K ~ mass wire and elec- 
7200 4,37 0.0391 5.9 4.6 
3550 0.0146 1.8 1.3 
3200 28.2 0.0246 0.8 0.7 
: Table 2 Composition of air samples after heating in bomb 2 and cooling at constant volume 
Initial Initial 
temperature, density, —Composition by volume percentage of-————--—-——-—-~ 
amagat Ne H.O A Co, NO 
\ir sample be- 
fore heating, “4 
300 K 1 76.7 20.0 2.4 0.93 0.03 no trace < 0.1 
7200 4.37 78.1 19.0 0.9 0.96 0.41 no trace <0.1 
5350 a 1.40 79.1 18.7 0.5 0.95 0.74 no trace <0.1 
3500 _ : 14.6 81.3 15.1 0.5 0.95 2.14 no trace <0.1 
3100 28.2 80.8 15.1 0.6 0.95 2.50 no trace <0.1 | 


The mass of electrode material eroded from very high cur- 
rent are electrodes is given by (9) 


The effects of the impurities of the vaporized tease wire and 
electrode surfaces were determined by collecting and weigh- 
ing the dust after a closed bomb 2 run. In these runs, all 
openings into the bomb were closed, and the hot air permitted 
to cool by heat transfer to the cold walls. The cooled air did 
not change in density from the initial density before heating. 
These dust impurities, which were weighed with a chemical 
balance, are found in Table 1. The expected dust per cent of 
gas mass from the known weight of tease wire and electrode 
erosion (from Equation [6]) are also given in Table 1. 
The greater contribution came from the vaporized electrode. 
There is substantial agreement between the measured and ex- 
pected dust impurities. The maximum dust impurity comes 
at low gas masses and amounts to about 6 per cent. Fastax 
camera pictures at 14,000 frames/sec indicated that it takes 
over 10 millisee for the dust impurities to diffuse from the 
electrodes into the body of the heated air and affect the hot 
air experiments. 

The cooled gaseous products from the closed bomb 2 runs 
were then analyzed with a mass spectrometer (Model 20-610 
Consolidated Engineering Corporation) to determine the 
composition of the air after spark heating. This technique 
will detect any impurity present as permanent gases. The 
results of the mass spectrometer analysis are found in Table 2. 
It will be seen that there is a rise in carbon dioxide level in the 
air, to a maximum of 2.5 per cent by volume. This mass 
spectrometer did not respond to small traces of nitric oxide. 

According to equilibrium theory, some nitric oxide will be 
formed at the temperatures achieved in the hot air. Then if 
the air experiment involves cooling the air, the kinetics of the 
nitric oxide reaction (10) indicate that, at a finite cooling 
rate, a part of the nitric oxide will be frozen in and will re- 
main as a permanent impurity. 

Calculations of the equilibrium concentrations from Hirsch- 
felder (11), of nitric oxide at high air temperatures, will be 
found in graphical form in Fig. 9. These theoretical results 
are compared with data showing the result of bubbling the 
cooled air through several water-collecting towers in series 
(Bubble Tower) and analyzing the nitric acid produced with a 
pH meter. The nitric oxide gas was also collected on other 
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occasions by scrubbing with water (Scrubbing Tower). From 
the known pressure decay in the closed bomb, a rate of 6000 
K/sec was measured at 2000 K, the most probable freeze-in 
temperature. The work of Steinmetz (10) indicates a yield of 
about 1.3 per cent at this rate, and the measured quantity of 
nitric oxide approximates this value. Using cooling rates of 
500,000 K/sec, Hertzberg (12) reports that 51 per cent of 
Steinmetz’s yield was realized. 

The determination of temporary impurities was also made 
with the help of bomb 2. A question arises concerning the 
state of the air after sudden cooling (to 106 K/sec). These 
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Table 3 Relaxation time upper limits and production of 
nitric oxide produced by 10° K/sec cooling in bomb 2 
Initial Initial Cooled Time to Per cent 

tempera- density, gas tem- reach final nitric 
ture, amagat perature, pressure, oxide by 
c K millisee volume 
7200 4.37 1170 1 0.65 
5350 7.76 579 1 0.67 
3550 14.6 702 1 0.93 
3200 28.2 507 3 0.89 


data are important if rapid isentropic expansions are to be 
studied. To investigate this condition, the hot gas was 
dumped into a chamber with ten times the initial volume of 
hot gas and filled with cold gas (300 K) at the same density. 
The 0.1 millisec rise time pressure probe was used to record 
the pressure rise in the large chamber. Data resulting from 
tests of the cooled gas are presented in Table 3. Also shown 
is the time necessary to equilibrate the pressure to final pres- 
sure (of the large chamber heated gases) within 3 per cent, 
and the per cent of nitric oxide determined by water absorp- 
tion. At this cooling rate, the expected yield of nitric oxide 
from (10) is 2.5 per cent, and Table 3 shows about one third 
this yield is achieved. 

A technique has been described by which any gas can be 
heated electrically in a constant volume and constant density 
process to temperatures up to 9000 K and pressures up to 40,- 
000 psi. Overall heating efficiencies of air are as much as 
100 per cent even at the highest temperatures. The natural 
cooling times of these gases in laboratory size vessels are on 
the order of milliseconds, permitting various gasdynamic ex- 
periments on this time scale. The possibilities of raising the 
gas temperature are only limited by the size of the electrical 
power supplies available as compared with the losses from the 
heated gases. Due to the fact that the heat power inputs to 
the electrodes are so high in these times, some impurity con- 
taminates the heated gas. The impurity rises as the gas is 
held for longer times, appearing in various experiments in 
about 10 millisec. The design and use of these hot gas 
sources is then a compromise between sufficient testing time to 


Conclusions 


the extreme temperature and heat power inputs. 


perform a given experiment and the rise in impurity level due 
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Measurement of Turbulent Heat Transfer Rates on the 
Aft Portion and Blunt Base of a Hemisphere Cylinder 


Turbulent heat transfer rates on the aft portion and on 
the blunt base of a hemisphere cylinder were measured in 
the 2; X 2] in. GALCIT shock tube over a range of shock 
Mach numbers between 3.25 and 5.1 and initial pressures 
bet. een 3 and 17 em Hg. The local Reynolds numbers on 
the cylindrical afterbody varied between 3.5 < 10‘ and 3.0 
xX 10° per em. A side support was used for the model to 
eliminate the disturbing effect of a rear sting support. 
The measured turbulent heat transfer rates on the cylin- 
drical portion agreed very well with previous flat plate 
measurements for small temperature differences, although 
the ratio of stagnation to surface enthalpy varied between 
3to 8 in the present tests. Only a slight effect of this large 
variation in h,/h,, was detected in this range of local Mach 
numbers, i.e., 1.25 < M, < 1.5. The measured heat transfer 
rate on the base indicated that at the center of the base the 
heat transfer rate is comparable to that on the surface just 
ahead of the base, while the heat transfer rate falls off to 
one-half to one-third of this value toward the rim of the 
base. This unexpected distribution of heat transfer rate 
over the base, and particularly the high value at the center, 
shows the necessity for a careful study of wake phenomena. 


Nomenclature 

Cy = Stanton number = 

haw — hw) 
h = enthalpy, cal/gr 
M = Mach number 
Nuz = Nusselt number = Pr+ Rex 
Rez = Reynolds number = pete* 4 

Ke 
Pr = Prandtl number = “2#¢ 

R = body radius, em 
Pp = pressure, cm Hg abs 4 
= heat transfer rate, cal/em?/sec 
r = recovery factor 
u = velocity, cm/sec 
f) = boundary layer thickness, em 
p = density, gr/em? 
m = viscosity coefficient, gr/em-sec 
Subscripts 
e = value at the outer edge of the boundary layer 
1 = ahead of shock front ear 
2 = behind shock front 7 
8 = stagnation 7 
aw = recovery 


w = wall ; 
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in the Shock Tube 
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HE use of the shock tube to simulate various aspects of 

hypersonic flight is now widely accepted. Recent develop- 
ment of the film resistance thermometer has made possible 
accurate heat transfer studies in the extremely short testing 
times available in the shock tube. Measurements of laminar 
heat transfer rates at and near the stagnation region of a blunt- 
nosed body (1 and 2)* indicated a very good agreement with 
the equilibrium laminar boundar layer theory for hyper- 
sonic flight speeds as developed by Lees (3) and Rose and 
Riddell (4). The present work extends these measurements 
to the aft region of a blunt body, i.e., to the aft surface and to 
the model base. For reasons which will be explained in detail 
in section 3, it is very difficult to obtain undisturbed laminar 
boundary layer flow on the aft portion of a long model in the 
present straight section of the 27 X 23 in. shock tube. Thus 
the work was confined to tarbalent ae at transfer studies. 

Published experimental data on turbulent heat transfer 
rates with large temperature differences across the boundary 
layer are scarce. In a recent report, Libby and Cresci (5) 
compare their turbulent heat transfer rate data measured on 
a blunt-nosed body in a wind tunnel with several analyses and 
find the data to agree best with the local flat plate approxima- 
tion. However, they suggest that their result should be sub- 
stantiated at a higher stagnation enthalpy, since their wind 
tunnel tests were limited to ratios of stagnation enthalpy to 
surface enthalpy of 1.3 to 2.0. The shock tube enables us 
to realize much higher ratios of stagnation to surface enthalpy, 
and in the present work this ratio lies between 3 and 8. The 
results of the present investigation tend to confirm the pre- 
vious conclusion of the validity of the flat plate analysis. This 
result, which is at first somewhat surprising, fits well into the 
overall trend indicating no startling or unexpected phe- 
nomena associated with high temperatures at ordinary pres- 
ures. 

The heat transfer problem on a blunt base is considerably 
more difficult. Much experimental and theoretical work has 
been done on the base pressure problem from the point of 
view of an overall “mixing” analysis. However, the heat 
transfer rate probably depends also on local conditions, so 
that a knowledge of the flow field near the base surface is re- 
quired before local heat transfer rates can be predicted. The 
-measurement of base heat transfer rates in the shock tube is 
also complicated by the question of the time required to es- 
tablish a steady wake behind a model. Thus the present ex- 
periment has two main objectives: (a) To 


- established in the time interval vail ible in the heed k tube 
(about 400-600 sec). (b) To investigate heat transfer rates 
on a blunt base when the body is supported from the side, in 
order to eliminate the usual disturbing effect of a rear sting 


support. 


3 Numbers in parentheses indicate References at end of paper 
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The shock tube and the instrumentation used in this investi- 


io r 4 gation are described in (1); therefore, only a brief summary 
ie —; will be included in section 2. Section 3 presents the experi- 
“| mental results, including a discussion and indication of future 
| work on this problem. 
Rubber Observation x The author would like to thank Prof. L. Lees for suggesting 
Section Windows Sein etna the problem and for his many helpful discussions and sugges- 
And Model Supports tions. 


Fig. X in, GALCIT shock tube 
Experimental Technique 


Shock Tube and Instrumentation 


This investigation was performed in the 23 X 2} in. GAL- 
CIT shock tube (Fig. 1). Heat transfer rate and schlieren 
studies were carried out with the model placed at the 19-ft sta- 
tion of the straight section. The instrumentation block (ia- 


Start Stop t gram is shown in Fig. 2 (see (1)). Shock speed is determined 
O x by measuring the time for the shock to travel between two 
reg monitoring stations spaced 2 ft apart. The heat trans/er 
Y gages are used for detection, and the time is measured by the 

7360 Berkeley counter. The schlieren spark is triggered 
through a time delay, which is activated by the downstre:m 
wave-speed monitoring gage. 
ia The heat transfer gages are sputtered directly on the 
| A model. Their initial resistance was measured by a Wheat- 
— stone bridge and the current was measured by a milliammeter. 
|. D.C. Power Supply The gage output is fed through a Tektronix 121 wide band 
2. Heating Ponels amplifier to a 585 Tektronix oscilloscope, and the records are 
3. Resistance Bridge 
obtained by a polaroid camera. 
Instrumentation diagram The gages are calibrated for quantitative heat trans/er 
measurements utilizing the method described in (1). The 
7 overall accuracy of the heat transfer instrumentation is 


a is about +5 per cent as demonstrated in (1). 
The Hemisphere Cylinder Model 

The model consists of a brass cylinder ? in. in diameter and 
22 in. long, with a hemispherical nose. A 3%; in. spike ? in. 
long is attached to the nose in order to produce a separated 
region and thus insure a fully developed turbulent boundary 
layer upon reattachment of the flow to the cylinder surface 
(Fig. 3). The aft portion of the model consists of a hollow 
Pyrex glass “cup” 2 in. in diameter and 1 in. long, which is 
waxed on the mating brass cylinder. The outside diameter 
of the glass and brass parts are matched to provide a con- 
tinuous smooth surface, and the platinum film resistance 
gages are sputtered on the glass cup (Fig. 3). Four gages are 
positioned on the flat base—gage 1 at the center and gages 2, 
3 and 4 on the circumference of a 3 in. circle as shown in Fig. 4 
Three additional gages are positioned on the cylindrical sur- 
face 0.082, 0.442 and 0.772 in., respectively, from the base, 
and are denoted as gages 5, 6 and 7 on Fig. 4. 

The model is side-supported by a 45 deg swept-back, double 
wedge that is 2 in. wide and 0.100 in. thick. This support 
enables schlieren observation of the aft part of the model and 
portions of the wake, and eliminates the need for a rear-sting. 

The leads on the glass cup are provided by a thick layer of 
silver paint to which wires are soldered on the unexposed 
edge of the glass cup (Fig. 3). The wires are then led 
through machined grooves in the model surface, then along 
the side support to the outside of the shock tube. The cup 
is aligned so that gages 3, 5, 6 and 7 are opposite the wake 
; of the side support in order to minimize the effect of this 
a wake on the readings of these gages (Figs. 3 and 4). 


7 Side view (b) Back view 


Fig. 3 Hemisphere cylinder model a. 


GAGE 


3 Experimental Results 


Range of Experiments 


| 
a A preliminary investigation was undertaken to determine 
La . the maximum size of the heat transfer model permitted by 
Ras Fig.4 Diagram of gage positions on the model the blockage effects—specifically by the effect of the reflection 
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Fig. 5 Schlieren studies of flow over the aft portion of the model. 


(a) (b) (ec) 


(a) Ms = 3.8, p, = 11.7 cm Hg, At = 40, sec. (b) Ms = 3.8, 


pi = 11.3 cm Hg, At = 140 yu sec. 


by the bow shock wave from the shock tube wall. It was de- 
termined that a 2 in. diam body can be accommodated in the 
22 x 22 in. section without choking the flow. For this model 
the ratio of blocked area to the tube cross section area is 
0.07, including the projected area of the side support. How- 
ever, the flow Mach number is only 1.3 to 2.0, which presents 
a more difficult problem so far as model length is concerned. 
The bow wave angle is relatively steep, and the first shock re- 
flection from the shock tube wall reaches the model at a sta- 
tion between ? and 13 in., a measured from the model nose. 
Thus it is impractical to attempt to test a model which is 
shorter than the distance of the first shock reflection from the 
nose. The model was designed to extend about 1 in. behind 
the first shock reflection, thus anticipating a range of shock 
strengths within which the second shock reflection from the 
tube wall will be weak enough and will strike the wake far 
enough downstream of the base so as to have little or no effect 
on the wake structure and the base heat transfer rate. 

The flow around the actual heat transfer model is shown in 
the schlieren photographs in Fig. 5. Figs. 5(a) and 5(b) show 
the developing flow at At = 40 yw sec and 140 u sec, respec- 
tively, for a shock Mach number M, of 3.8. Fig. 5(c) shows 
the flow at At = 350 usec for WM, = 3.0. Fig. 5(b) indicates a 
fully developed wake and a relatively weak second reflection, 
while at the lower Mach number, Fig. 5(c), the flow is ob- 
viously choked. At high Mach number, /, > 4.5, the re- 
flected shocks again increase in strength, and a choking effect 
again appears. Thus the base heat rate measurements were 
limited to a narrow range of shock Mach numbers (3.5 < M,< 
4.0) because of the particular combination of flow Mach 
number and test section geometry. The initial pressures 
were varied between 3 and 17 cm Hg corresponding to 
Reynolds numbers/em in the range 3.5 X 104 to 3 X 10°. 
The lower limit is imposed by the minimum pressure required 
for a fully developed turbulent boundary layer, and the upper 
limit is set by structural considerations. The advantages of 
larger tube cross section dimensions and of higher flow Mach 
number are obvious. These considerations indicate the ad- 
vantages of using an expansion nozzle as the test area for this 
type of investigation. 

The limitations on the range of shock Mach numbers for the 
turbulent heat rate measurements on the cylindrical surface 
are not as severe, because the expansion around the shoulder 
limits the feedback from the wake. Good results are obtained 
for shock Mach numbers in the range 3.25 to 5.1. 
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(c) M, = 3.0, p, = 12.5 cm Hg, At = 360 uw sec 


The disturbance caused by the shock wave reflection on the 
model surface ahead of the heat transfer gages causes dis- 
turbances in the laminar boundary layer when one attempts to 
obtain laminar heat rate measurements. However, the effect 
is negligible for turbulent boundary layers, so the present in- 
vestigation was confined to the turbulent case. 


Turbulent Heat Transfer Rate Measurements on _ the 
Cylindrical Surface 


Representative gage outputs at various positions on the 
model are shown in Fig. 6. Fig. 6(d) is a record of gage 7 on 
the surface of the cylinder. Measured heat transfer rates from 
gages 5, 6 and 7 at .W/, = 5.8 are shown in Fig. 7 as a function of 
initial pressure p;. The lines drawn‘ show a relation g~ 
The heat transfer rate measured by gage 1 on the base is also 
shown for comparison. These base measurements will be dis- 
cussed in the next section. 

In Fig. 8 the heat transfer rates measured at the various 
locations on the cylindrical surface of the model are plotted as 
g/p:°8 vs. M, (were p; is measured in em Hg abs and ¢g is 
cal/em?/sec). For comparison the stagnation point heat 
transfer rate for p, = 1 cm Hg abs calculated for the equi- 
librium laminar boundary layer according to (3), is also 
shown. Turbulent heat transfer rates on the cylindrical after- 
body are one quarter to one third of the laminar stagnation 
point heat transfer rate for comparable free stream conditions. 

Turbulent heat transfer rates are generally expressed in 
terms of the Stanton number Cy where 


Semi-empirical correlations of the available experimental data 


where a = —0.2, 6 = —2, and B depends on the local Mach 
number and the temperature (or enthalpy) ratio across the 
boundary layer. Also 


hew = hal — ( 


4The actual Mach number varies between 3.5 and 4.0, and the 
heat transfer rates are normalized to M, = 3.8 by the factor 


[M/(M,)|*. 
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effective distance ¥ is determined directly from the measured 


(a) (b) 
Fig. 6 Representative gage output at various positions on the model. (a) Gage No.1: 1/; = 3.5, p, = 8.0 cm Hg, sweep = 100 y 
sec/div, Ry) = 24.9 2, Jo = 20 ma, sensitivity = 10 mv/div. (b) Gage No.3: Ms; = 3.6, p,; = 6.6 cm Hg, sweep = 100 u sec div, 
Ro = 19.0 2, Io = 20 ma, sensitivity = 5 mv/div. (c) Gage No.4: M; = 3.5, p, = 6.9 cm Hg, sweep = 100 yu sec/div, Ry = 20.5 o 
Io = 20 ma, sensitivity = 5mv/div. (d) Gage No.7: M; = 3.81, p, = 6.86 cm Hg, sweep = 100 u sec/div, Ro = 19.32, 1, = 20 ma, 
sensitivity = 10 mv/div 
100 where the recovery factor r is approximately Pr’/*; i.e., = 
SS as | | | = 0.89 for Pr = 0.7. According to Equations [1, 2] 
GAGE | NO7 
60] GAGE NO.5 | NO6___ = — hy) [3] 
0 GAGE NO.6 
Vv GAGE NO.7 re NO.1,5 where Z is some effective length of run for the fully developed 
40 re ee Oa oz turbulent boundary layer. In other words in the present ex- 
eZ | periments the quantity g¥°-?/p,°* should be the same for all 
30 | three gage locations at the same shock Mach number. ‘This 
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Fig. 7 Turbulent heat transfer rate variation with pressure at 
(Ms) = 3.8 
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Fig. 8 Turbulent heat transfer rate variation with shock Mach 
number 
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values of g and turns out to be remarkably constant for all 
three gages over the range 3 < M, . 5.1, with the following 
values: ¥; = 0.50 in., %, = 0.82 in., ¥; = 1.19 in. within +5 
per cent. In Fig. 9 the quantity pay 2/p,0 8 is plotted against 
shock Mach number; evidently the “spread” between the 
three gages shown in Fig. 8 is practically eliminated. 
Previous measurements of turbulent heat transfer rates 
with zero pressure gradient are correlated very well by Equa- 
tion [3], where B = 0.03 for M, = 0 with T,,/T,=1; B= 
0.0266 for M, = 0.87, and B = 0.021 for M, = 1.62, with 
T/T aw = 1, i.e., small temperature differences. These data 
are obtained from (6, 7, 8, 9, 10) and are reproduced in Fig. 
10. In order to compare the present results with this earlier 


5 This effective distance < for any gage is very nearly equal to 
the distance measured from the point on the model at which the 
bow shock reflected from the shock tube wall strikes the model 
surface (Fig. 5). 


40 
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Fig. 9 Comparison of turbulent heat transfer rate with fiat 
plate measurements 
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work the “external” flow quantities (pu), and yw, for the 
present experiments were computed by assuming a normal 
shock ahead of the model nose and an isentropic expansion 
along the body surface from the stagnation point to the free 
stream pressure po. The turbulent heat transfer data ob- 
tained in the shock tube agree remarkably well with the earlier 
data of (6, 7, 8) in which 7.,/Taw ~ 1. In Fig. 9 the curves 
representing the average of these data (obtained from Fig. 10) 
at Mach 0, 0.87 and 1.62 are shown. The shock tube data 
which correspond to local Mach numbers between 1.25 and 
1.5 fall between the curves for M = 0.87 and M = 1.62. 
The remarkable fact here is that the corresponding ratio of 
stagnation to surface enthalpy increases from 3 to 8 in the 
present tests. 

When one plots the data in the form Nu; vs. Rez disregard- 
ing the variation in enthalpy level (Fig. 11), the scatter in the 
data increases slightly as compared to Fig. 9. This scatter 
indicates a slight effect caused by the variation of h,/hw. One 
can conclude that in this Mach number range‘ the effect of 
enthalpy level is surprisingly small. This conclusion is in 
agreeinent with the previous results of Libby and Cresci (5). 

The experiments of Libby and Cresci (5) can also be used 
to observe the effect of falling pressure gradients on turbulent 
heat transfer rates. These experiments were performed in a 
blowdown wind tunnel with varying stagnation pressures and 
temperatures during the run. In order to correlate the accu- 
mulated data, modified Reynolds and Nusselt numbers, based 
on stagnation conditions, were used. It is difficult without 
recomputing the data to obtain local heat transfer values. 
However, a comparison with flat plate analysis is presented 
in (5), from which the agreement with the local flat plate data 
can be examined.? The geometrical distance from the stag- 
nation point has been used in evaluating the Reynolds and 
Nusselt numbers. This procedure is consistent with most 
turbulent boundary layer studies on flat plates. The fact that 
it seems to give good results in the case of an axially sym- 
metric body with a falling pressure gradient is worth noting. 
Fig. 7 of (5) shows that the measured heat transfer rates fall 
within +10 per cent of the curve described by Equation [6] of 
(5). Translating the particular definition of the modified 
Reynolds and Nusselt numbers into their conventional form 
results in the relation 


Nuz = [4] 


where p’ and y’ are defined at reference conditions determined 
by the local pressure and reference enthalpy 


h’ = 0.5hy + 0.22h,, + [5] 


Now for the Libby and Cresei experiments h,/h» varies be- _ 
tween 1.3 and 1.7; thush’ =h,. The local Mach number at 
the thermocouple positions in these experiments varied be- 
tween 0.5 to 0.8. Thus one observes that these data agree 
with the flat plate data with no pressure gradients at the 
corresponding Mach numbers. 

One can conclude from the collection of data presented here 
that for subsonic and low supersonic flow Mach numbers, the 
turbulent heat transfer rates seem to be to a surprising degree 
independent of falling pressure gradients, stagnation to wall 
enthalpy ratio, and radius of curvature when 6/R <1. The 
only appreciable effect is found to be the local Mach number, 
which is responsible for about 15 per cent reduction in the 
heat transfer coefficient by the increase of local Mach number 
from 0 to about 1.5. 


‘This result does not imply that the effect of hs/hy is small 
also at high local Mach numbers. 

7 Tt may be that some of the scatter of the experimental results 
may be due to the correlation with respect to stagnation rather 
than local conditions. A replotting of the data in terms of the 
conventional local quantities would be very helpful in generalizing 
these data to other heat transfer applications. 
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Fig. 10 Measurements of turbulent heat transfer rates on a 
flat plate, 7) > Tw (6, 7, 8, 9, 10) 
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rates with the flat plate results 


Turbulent Base Heat Transfer Measurements 


The arrangement of the four heat transfer gages on the 
blunt base is shown in Fig. 4, and representative outputs of 
gages 1, 3 and 4 are shown in Figs. 6(a), 6(b) and 6(c). The 
output of these gages is not as smooth as that of the surface 
gages 5,6 and 7. However, the output of gages 1 and 3 indi- 
cates a steady heat transfer rate after a “building-up”’ 
period of about 100-150 uw sec. Gage 4 lies in the wake of the 
side support, and the output of this gage never gave a steady 
heat transfer rate. The comparison of the base heat rate at 
the center and at a location } in. from center as measured with 
gages 1 and 3 is shown in Fig. 8. Gage 1 at the center of the 
base indicates a heat rate comparable with that of gage 5 on 
the surface just ahead of the base. At the gage 3 position the 
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heat rate is only one half to one third of this value. One 
should note that the output of gage 4 (though disturbed) is 
comparable to that of gage 3, which increases one’s confidence 
in this result. 

The experimental data indicate that a “steady” heat trans- 
fer rate is established after about 100-200 u sec. However, 
the total testing time is about 400-600 yu sec, which may not 
be long enough to establish steady equilibrium wake flow such 
as may be experienced in a wind tunnel or in free flight. Until 
this question is settled the present base heat transfer data are 
considered to be preliminary and are reported mainly to call 
attention to this result and stimulate further thought and re- 
search into this effect. However, the fact that these measure- 
ments are repeatable, that the heat rate is steady and uniform 
with time after about 100-150 u sec, and that the heat rate at 
the base center varies like p;°- may be used as arguments to 
indicate that these effects are real. 

The unexpected distribution of heat transfer rate over the 
base, particularly the high value at the center of the base, 
shows the necessity for careful study of wake phenomena. 
One way to attack this problem is to perform similar experi- 
ments in the shock tube with a wide range of model sizes. 
If the observed effect is caused by the unsteady flow condi- 
tions, the size of the model should affect the results. Because 
of the limitation in size of the present shock tube, such a 
program could not be undertaken at this time. In order to 
shed more light on this problem, an investigation of base 
and wake flows is now under way in the GALCIT 5 in. con- 
tinuous-flow hypersonic tunnel. 
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Technical Notes 


A Comparison of Turbojets and Ram- 
i jets for High Speed Flight 


M. ARENS! 


Israel Institute of Technology, Haifa, Israel 


Nomenclature 
F, = thrust, Ib 
f/a = fuel-air ratio 
g = gravitational acceleration 
total pressure 
p = static pressure a? 
SFC = specific fuel consumption, lb/hr/Ib a 
T = total temperature 
V = velocity 
w = airflow, lb/sec 
y = ratio of specific heats 
n = efficiency 
7 = total pressure ratio 
T = ratio of total temperature to ambient temperature 
Subscripts 
b = burner 
c = compressor 
e = nozzle exit 
n = nozzleentry 
0 = ambient 
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t = turbine 

T.J. = turbojet 

d = turbine entry 
r = ram 


HE relative advantage of ramjet powerplants over turbo- 
jet powerplants in high speed applications has taken on 
renewed significance with the rapidly increasing flight speed 
range of current airplanes. If we limit our attention to the 
propulsive performance of the two powerplants, in terms of 
thrust per pound of airflow or specific fuel consumption, it is 
possible to clearly delineate an area of superiority for each 
powerplant. Since current applications of turbojets to high 
speed vehicles require the use of an afterburner, the com- 
parison is based on an afterburning turbojet, although the 
method can be applied to non-afterburning engines as well. 
The specific fuel consumption of an air-breathing engine 
can be expressed as 


If we assume that the effective combustion efficiency of the 
two-stage combustion process in the turbojet is equivalent 
to that of the ramjet combustion process, and neglect con- 
sideration of any air-leakage or bleed in the turbojet, then the 
fuel-air ratio to attain a desired nozzle entry temperature will 
be the same for both engines. Their relative superiority can 
therefore be measured in terms of thrust per pound of airflow 


w 


Eprror’s Nore: The Technical Notes and Technical Comments sections of Jer PROPULSION are open to short manuscripts describing new 
developments or offering comments on papers previously published. Such manuscripts are published without editorial review, usually 
within two months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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Fig. 1 Optimum turbojet compression ratio 
orn terms of nozzle exit velocity, since the powerplants are — 
compared at the same flight speed. For a fully expanded 


nozzle? 


and therefore, for any desired afterburner or ramjet temper- 
ature, the engine providing the highest nozzle pressure ratio 
isthe more efficient. If it is assumed that the total pressure 
losses in the turbojet afterburner are equivalent to those of 
the ramjet combustion chamber, the requirement for turbojet 
superiority is 


or that the turbojet “gas generator” produce an increase in 
pressure. 

It is obvious that for any flight speed and turbine-entry 
temperature there is an optimum compression ratio subject 
to assumed component efficiencies. Equation [5] gives the 
gas generator pressure ratio 


By introducing the equivalence of compressor and turbine 
work, as well as their respective efficiencies, one obtains 
¥/(v¥-1)) 
(y-1)/ 
Wr.j. = TH, | 1 — —— (2-‘7 — 1) [6] 


The optimum compression ratio can then be found by differ- 
entiating Equation [6] with respect to 7. and setting the re- 
sult equal to zero 
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Fig. 2 Optimum turbojet gas generator pressure ratio 
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Fig. 3 Regions of turbojet and ramjet superiority 


TANCHe 


1/(y—1) 


OT. \ 
in which case the optimum compression ratio is 


1 manent \y/(y—1) 
om.(optimum) = + [8] 
and the condition for turbojet superiority is that Equation 


. Sutton, G. P., “Rocket Propulsion Elements,” Wiley, 2nd 
edit., p. 53. 
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[6] be greater than one when the optimum compression ratio 
is used. Figs. 1 and 2 show the optimum compression ratio 
and gas generator pressure ratio as a function of flight speed 
and turbine-entry temperature, while Fig. 3 indicates the re- 
gion of turbojet superiority. It may be concluded from Fig. 
3 that the ramjet is the appropriate powerplant for a signifi- 
cant portion of the flight spectrum. 
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test stand 


Simulation of Free Flight Structural 
Characteristics by Captive Missile 


I. Computer Analysis 

HE limitations of conventional test stands for preparing 

missiles to withstand the vibration environment o/ free 
flight are revealed by a computer study (1)? of a missile 
structure. A test stand mount that is soft and has «2 low 
damping characteristic is quantitatively shown to reduc: the 
error between free and captive flight tests; however, struc- 
tural considerations limit the degree of softness that cin be 
obtained for a particular sized missile. 

The purpose of this paper, which was inspired by a previous 
report by Lee and Miesse (3), is to determine the exteit of 
difference between the vibration environment of missili's in 
captive and free flight tests during single frequency excit:tion 
and to suggest a method of obtaining captive test data that 
will accurately predict free flight performance. At present, 
until the nature of the vibration environment is better de- 
fined, the advantages of the classical steady-state vibration 
theory can be used. Although aerodynamic and rocket 
forces are random in nature, they are filtered by structural 


P. LIEBERMAN! 


Armour Research Foundation, Chicago, III. 


Received April 11, 1958. 
1 Associate Research Engineer, Combustion Research Section. 
? Numbers in parentheses indicate References at end of paper. 


4 
’ 
10? 


Ix to! 
WwW 
\ 
5 
WJ \ 
\ 
IxlO” 
\ 
ixtor® 
3 
Oo Ix 
2 \ 
\ 


N 
0.1 1.0 10 


FREQUENCY OF VIBRATION 


100 1000 


CYCLES PER SECOND 
——-— Free flight Ks = 1.16 X 106 Fs 
lb 
Captive flight C; = — 
fps 
Fig. 3 Frequency response of typical missile in soft mount test 


stand 


JET PROPULSION 


men 
quel 
Thu: 
valid 
vibré 
mod 


resul 


missi 
inter 
Fig. 


respo 
in Fi 
mour 
appa. 
siders 
opera 
violet 


In 


ROCKET MOTOR DISPLACEMENT — FreEGrT 


Fig. 4 


SEPTED 


Ma 
degre 
* 
: test 
; 
; 
+. 
: 
j | 
— 


e 


members and emerge as vibration confined to a narrow fre- 
quency band with an apparent fluctuating amplitude (2). 
Thus the classical single frequency excitation technique is 
valid for any degree of damping. Although the missile can 
vibrate in axial, bending and torsional modes, only the axial 
modes were considered. However, the test stand restraint 
resulting from an axial analysis is also applicable to the other 
degrees of freedom. 

The variational portion of the rocket motor thrust causes a 
missile to vibrate as though it were a multimass spectrum 
interconnected by springs and dash-pots, as illustrated by 
Fig. 1. Restraint of the missile vibratory motion by a static 
test stand causes thrust energy to be reflected and absorbed 
ina manner which produces considerable difference between 
free flight and test stand frequency response. 

A comparison of free flight and captive missile frequency 
response to an oscillatory force of 1000 lb amplitude is shown 
in Fig. 2, for which the spring constant of the stiff missile 
mounts was set up at 1.16 X 10° lb/ft. It is immediately 
apparent that missile mounts of this stiffness introduce con- 
sideral le undertesting of the system, so that a missile which 
operates satisfactorily on the test stand is apt to experience 
violent destructive vibrations in free flight. 


In an effort to avoid the errors shown by Fig. 2, various 
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agencies have used nylon cord, cable and springs as means of 
restraint for small missiles in order to simulate the air cushion 
experienced in free flight. For large missiles, the steel torsion 
bar is thought to be a practical approach to soft mounting. 
The beneficial effects of soft mounting are vividly portrayed 
by Fig. 3, which assumes a spring constant of 1.16 X 10° lb/ft. 
For frequencies above 6 cps, it is apparent that the soft 
mount gives excellent simulation. If the missile vibrates at 
frequencies lower than 6 cps, however, then the high ampli- 
tudes must be restrained by auxillary damping and servo- 
control systems that will inhibit the simulation but save the 
test stand for future tests. 

The soft mount test stand is adequate for a range of 
missile sizes, as well as under the condition of changing missile 
mass during firing. Fig.4 shows that the placing of an IRBM- 
type missile on a soft mount test stand designed to be optimum 
for an ICBM-type missile, still permits acceptable simula- 
tion for the former. Another extreme condition occurs when 
the missile empties 85 per cent of its mass by ejecting propel- 
lant. As shown in Fig. 5, the soft mount designed for full 
tanks continues to yield adequate simulation at the end of the 
firing. Apparently the smaller missile masses on the flexible 
mount raise the missile stand resonant frequency; however, at 
the higher frequency, inertia reduces the amplitude discrep- 
ancy between free and captive frequency response. 

As a result of this investigation, the following conclusions 
‘an be drawn: 

1 Ground testing of missiles in conventional test stands is 
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an inadequate method for preparing missiles to withstand the 
free flight vibration environment. 

2 For frequencies above 6 cps, the soft mount gives exeel- 
lent simulation of free flight. 

3 Auxillary dampers and a control system need to ae. 
company a soft mount if low frequency, between 0.01 and 
6 eps, excitations occur. An uncontrolled low frequency 
excitation can tear the missile away from the test stand. 

4 A soft mount test stand designed for a large missile with 
full tanks provides adequate free flight simulation for smaller 
missiles, as well as the large missile with empty tanks. 
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Aircraft Gas Turbine Engineers 


MECHANICAL DESIGN | 
AERO-THERMODYNAMICS 


These are two with 


lines of responsibility. 

: You serve as a consultant on major (aerodynamic 
and thermodynamic) (mechanical design) prob- 
lems encountered in conceiving and developing 
aircraft gas turbines. 


In this capacity, you provide technical guidance i. 


on (aero-thermodynamic) (mechanical design) oa 
programs, contributing your broadest and most 


advanced engineering knowledge and under- 


standing. 


Specifically, you initiate and integrate the long- 
range planning required to exploit broad 
(aero-thermodynamic) (mechanical design) 


advances. You also promote acceptance of 
these plans by related engineering components, 


to assure coordinated engineering effort. 


xperience: 


You should have a minimum of 5 years’ re- ‘s 
+4 sponsible (aero-thermodynamics) (mechanical 
: “4 design) experience in aircraft gas turbines, plus 
other work in the field, from preliminary design ah 
through production. 
‘ 


in complete confidence to: 


Mr. William Merrill 


_ SMALL AIRCRAFT ENGINE DEPARTMENT 


1119 Western Avenue 
aa West Lynn, Mass. 4 
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New Patents__ 


George F. McLaughlin, Contributo; 


Convertible aircraft with jet-driven lifting 
rotors (2,823,875). A.C. Peterson, !dina, 
Minn. 

Means are provided for rotational pro- 
pulsion of the rotor unit or alternative 
translation propulsion of the aircraft. 
Fabricated gas turbine structures (2,824, 
359). E. C. Rhodes and D. W. Rhys, 
Hounslow, England, assignors to The 
International Nickel Co., Inc. 

Hollow blade of chromium alloy adapted 
to be subject to corrosive conditious pre- 
vailing at temperatures of about ‘00 to 
850 C. 


External turbine jet engine (2,831,320), 
W. E. Duncan, Miami, Fla. 

Casing integral with an end wail and 
rotatably mounted on a stationary sup- 
port. A fluid stream is directed :gainst 
the periphery of the wall and through fluid 
guiding passages which interact with the 
stream to rotate the casing. 


Devices for releasably holding the rear 
ends of rockets (2,831,400). H. Hosli, 
Geneva, Switzerland, assignor to Lrevets 
Aero-Mechaniques, S. A. 

Pressure from the rocket nozzle actuates 
a central piston connected to a pair 
pivoted jaws holding the rear of a rocket, 
opening the jaws and releasing the rocket. 


Automatic parachute release (2,831,721). 
R. Gross, Dayton, Ohio, assignor | 
USAF. 

Load release for shroud lines, consisting 
of an explosive charge in the lower end of a 
cylinder. Slow-burning power charges 
ignite the main charge upon actuation of a 
trigger, blowing a closure member away 
from a hook holding a ring attached to the 
load. 


Convertible turbojet and ramjet engine 


(2,832,192). N. N.  Budish, Seattle, 
W ‘ash., assignor to Boeing Air lane Co. 

Bladed gas turbine mounted for shifting 
in the axial direction between two posi- 
tions. In first position, blades project 
within the gas flow passage; in second 
position, blades are withdrawn, leaving the 
gas flow passage unobstructed. 
Gas turbine fuel control system responsive 
to — and compressor load (2,832,193). 
H. J. Wood, Sherman Oaks, Calif., as- 
signor to Garrett Corp. 

Means for subjecting to differential 
pressure a movable wall connected to op- 
erate a fuel valve. A valve in the line 


Eprtor’s Norte: Patents listed above 
were selected. from the Official Gazette of 
the U.S. Patent Office. Printed copies of 


patents may be obtained from the Com- 
missioner of Patents, Washington 25. 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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ASTRODYNE 


Astrodyne was founded early this year to bring to- 
gether in one productive group all the skills, experi- 
ence, and facilities it takes to design and manufacture 
America’s major solid propellant missile systems. 

From Phillips Petroleum Company came men 
with impressive backgrounds in research, design, 
and manufacture of superior solid propellants. From 
North American Aviation came men who know mis- 
sile systems management, and who have designed 


—aimed for the future in solid rocketry 7 


ASTRODYNE, INC. 


McGREGOR, TEXAS 


and built the largest rocket engines in use today. 

Today Astrodyne has the experience to design, 
develop, and manufacture complete propulsion sys- 
tems, extruded and cast propellants, solid propellant 
rocket motors, and boosters, and gas generator 
charges for auxiliary power units. 

Inquiries are welcomed on any phase of the solid 
propellant field—from preliminary design to quan- 


tity production. 
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and imagination to carry out individual studies and participate in group efforts on propulsion systems and 
installations for advanced aircraft and missiles. 


PROPULSION DEPARTMENT 


Propulsion Aerodynamicist 


MS in AE with at least five years’ experience in aerodynamics of fluid flow. Emphasis is on aero- 
dynamics of propulsion systems—particularly external flow effects on induction system and exit 
system design of high-speed vehicles. Knowledge of supersonic and hypersonic flow theory man- 
datory. Experience in the propulsion aerodynamics of VTOL and STOL aircraft desirable. Will be 
responsible for setting up applied research projects: and performing design and — work i in 
support of advanced design effects. 

BS or MS in AE, ME, Physics or Chemistry with 3-5 years’ experience in rocket and/or ramjet 
propulsion system analysis and system application. Additional experience involving supervision, 
military requirements, customer liaison, research, engine development, or testing is desirable. 


Engine Performance Analyst 


BS or MS in AE or ME with 3-5 years in turbojet and ramjet performance analysis. Experience re- 
quired with IBM 704 programing application to engine performance. 


Heat Transfer Engineer 

BS or MS in AE, ME or Physics, with 3-5 years of research or analytical experience in applied heat 
transfer or thermodynamics. Work involves aircraft air conditioning, aerodynamic heating, engine 
installation cooling, and ice removal. 


Applied Mathematician 

Master’s Degree in mathematics with experience in applied mathematics or mechanics. Work 
involves studies and investigations in the fields of heat transfer, aerodynamics, thermodynamics, 
propulsion and statistical computational techniques. 


EQUIPMENT INSTALLATION DEPT. 


Fuel System Engineer 
BS in ME, CE or AE or equivalent with 2-5 years’ experience in fuel system design for recent military 

S- jet aircraft. Technical knowledge and experience in the following fields required: 
- Basic jet aircraft fuel system design and analysis. 

a. Design of fuel system components and accessories and attendant vendor and military service 
liaison problems. 
4 .. Design and installation engineering of internal self-sealing, bladder and integral fuel tanks. 
_= ye Design of pressure fueling and flight refueling systems and flight refueling packages. 


Power Piant Installation Engineers 


Graduate engineers with 5 or more years’ experience in the design and installations of turbo-jet, 
turbo-prop, ramjet and rocket engines for aircraft and missile propulsion systems. Additional 
experience in supervision, manufacturing, or test development would be beneficial. 


U. S. Citizenship required. We pay moving expenses. 
For immediate interview, send your resume to: 


MR. A. WILDER, Engineering Personnel Director 


New projects at GRUMMAN AIRCRAFT have created responsible positions for graduate engineers with initiative 


connected to the downstream side of { 
throttle modulates the pressure level aet 
ing on the wall. 


Airplane design (182,597). A. Kurtveli, 


Huntington, N. Y., assignor to Republ 
Aviation Corp. 

Twin-engine single-seater jet with ¢ 
hausts close to the fuselage, below 1 
wing roots. 


Ramjet wing system for jet propelled 
aircraft (2,836,379). B. T. Salmon, San 
Diego, Calif., assignor to General Dynam- 
ics Corp. 

Ramjet propulsion integrated with tl 
wings of an aircraft as to constitute : 
power plant which is compact in size al 
efficient in operation at supersonic speeds 
The configuration propels and serves 
sustain the aircraft in flight. 


Missile launcher (2,837,970). H. A. 
Helstrom Jr., Dallas, Tex., assignor to 
Chance Vought Aircraft, Inc. 

Device for successively launching roe 
ets. Missiles are retractably mounted 
guides, and positioned for launching | 
chains and sprockets actuated by & 
electrical motor. 


Hydraulic ejection equipment for missiles 
(2,837,971). R. Wosak, Quaker Hi 
Conn., assignor to General Dynamics 
Corp. 

Balancing system for ejecting missiles 
from a submarine or other installatiot 
under the sea. A passageway connects 
a water cylinder with the sea wherel 
the hydraulic pressure from the sea is tl 
same in the cylinder as it is at the muzz 
end. 


Automatic emergency radio control system 
(2,838,255). J. J. Hagopian, San Jo 
and s. R. Hennies, Los Altos, Cali 
assignors to Northrop Aircraft Corp. 
Control circuit for an aircraft to | 
initially launched and connected to ope 
ate a control member in accordance with 
radio signals from a manually operated 
remote transmitter. A time delay relay 
operates a reset switch after the sw itch 
is energized for a predetermined time in- 
terval. The relay is energized simultane- 
ously with launching of the aircraft. 
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research and 


advanced 
development 


Development of a re-entry simulation facility for 
the Air Force Titan ICBM nose cone required 

a method for heating air to extremely high tem- 
peratures. A special type of electrical discharge, 
sometimes called a ‘‘plasma jet,” was considered. 
Both fundamental research—to obtain a 
quantitative understanding of the are properties 
and mechanism—and applied research in actual 
aerodynamic and materials testing were needed. 
Advanced development—the actual design 

and construction of improved test facilities— 
was also an inherent part of the problem. a 


Theoretical and experimental physicists, 
aerodynamicists and physical chemists began 

a simultaneous attack on the fundamental 
processes. At the same time, and in close 
association with the others, instrumentation 
and mechanical engineers began to build and 
use a series of continually improved facilities for 
re-entry testing. Those inclined toward either 
fundamental or applied research in their particular 
specialty were able to maintain their integrity 
and their own research approach. At the 

same time they were able to benefit from the 
close working association with one another. 


The results were especially fruitful. The plasma 
jet has proven an excellent tool for re-entry 
studies. Investigations of important applications 
in high-temperature chemistry, space propulsion 
and refractories are following as an added bonus. 


Research and Advanced Development is more 
than a descriptive title at Avco. It is a concept. 
that promotes creativity. 


For information on unusual career opportunities 
for exceptionally qualified scientists and engineers, 
write to: Dr. R.W. Johnston, 

Scientific and Technical Relations, 

Avco Research and Advanced Development Div., 
201 Lowell Street, Wilmington, Mass. 


— 


Sk. PERFORMANCE ANALYST 


TURBOJET MAC HINERY 


(AE or ME) 


. a unique application of your skills on the 


AIRCRAFT NUCLEAR PROPULSION PROGRAM 


Now in an advanced development stage 
at General Electric, the ANP program 
has an important new assignment for an 
engineer with 6 to 10 years of pro- 
gressively more advanced experience in 


formance of turbomachinery with all- 
over power plant design for nuclear 
powered aircraft. 


He will formulate and direct perform- 
ance analysis projects to further these 


one or more of the following fields: cafes: es : 
c. objectives; coordinate the work of con- 


tributing groups: assess the effectiveness 
of analysis, design and development 
work involved. 


Propulsion System Performance 
Analysis; Aero-Thermodynamics of 
Turbomachinery: Power Plant In- 


stallation (Aero-Thermodynamics) Qualified Applicants, 


challenge of this pioneering assignment, 


are requested lo write lo: 


GENERAL @@ ELECTRIC 


Cincinnati 15, Ohio 


This engineer's primary responsibility 
will be to integrate, cycle-wise, the per- 


Mr. P. W. Christos, Div. 34-MU 


AIRCRAFT NUCLEAR PROPULSION DEPT. 


P.O. Box 132 


7 


“MONOBALL” 


Self-Aligning Bearings 


PLAIN TYPES 


ANALYSIS 


Stainless Steel 
Ball and Race 


RECOMMENDED USE 


For types operating under “ae temper- 
ature (800-1200 degrees F.). “ 


Chrome Alloy For types operating under high radial 
Steel Ball and Race {ter loads (3000- 893,000 ibs. 


Bronze Race and For types operating under normal loads 
Chrome Steel Ball with minimum friction requirements. 
Thousands in use. Backed by years of service life. Wide variety 
Bre Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
_ similar size range with externally or internally threaded shanks. 
mg Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
‘types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
_ have been used to meet specific needs. Write for revised Engi- 
_ neering Manual describing complete line. Dept. JP-58. 


SOUTHWEST PRODUCTS CO. 


MOUNTAIN AVE., MONROVIA, CALIFORNIA 


Book Reviews __ 


Ali Bulent Cambel, 


- 
Northwestern University, 
* ox Associate Editor 


Microwave Measurements, by Edwiid L, 
Giuyton, McGraw-Hill, New_York, 
1957, 476 pp. $12. 


Reviewed by ARTHUR BRONWEL| 
Worcester Polytechnic Institut: 


“Microwave Measurements”’ serve~ the 
very valuable function of bringin: to- 
gether in one volume ‘much of the siznifi- 
sant material on this subject, The 
author’s research experience in this field 
enables him to appraise the material 
accurately and also to bring this up to date 
with a wealth of information on recent 
developments, many of which have «ome 
out of researches of the laboratories which 
he directs at Stanford University. This 
book presumes a knowledge of basic ma- 
terial in electromagnetic theory. The 
derivations of the equations are not given, 
but most beginning texts in microwaves 
will supply this information. 

The treatment of klystrons goes far 
beyond that of other texts in supplying 
both practical information regarding de- 
sign and theoretical formulations for 
analyzing performance. It considers such 
aspects as klystrons used as oscillators, 
converters; amplifiers and frequency multi- 


COMPLETE 


TESTING 


FACILITIES 


* Qualification Tests 
Evaluation Tests 

* Performance Tests 
* Environment Tests 


Manufacturers of Metal Boss 
Seals to Military Specifications 


TESTING CO. 
1812 Fleet St., Baltimore 31, Md. 
ORleans 5-8337 ORleans 5-2222 
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pliers, as well as methods of obtaining 
tunable klystrons. There is a section on 
reentrant cavity oscillators and a short 
section on traveling wave tubes. The 
absence of information on the magnetron 
js unfortunate. 

The sections on the generation and de- 
tection of laboratory signals will be ex- 
tremely helpful to those engaged in 
jaboratory practice. The chapter on 
measurement of microwave power reflects 
the excellent work done on balometer 
bridges and calorimeters during the -war, 
which has not been much improved on in 
basic principle since. In the chapter on 
impedance concepts there is a rather shy 
treatment of Maxwell’s equations, fol- 
lowed by an extensive analysis of standing 
waves on transmission lines, a subject 
whicl: is further elaborated upon to good 
effect in the following chapter on imped- 
ance measurement. 

Subsequent chapters deal with the 
meastirement of wave length and fre- 
quency, as well as a rather thorough 
treatment of the measurement of Q of 
cavities. An interesting innovation is the 
inclusion of a section on the use of molecu- 
lar absorption as microwave frequency 
standards. 

One gets the impression that the author 
is much more at home with circuit equa- 
tions than with Maxwell’s field equations. 
Perhaps it is in the nature of the text, deal- 
ing as it does with the entire gamut of 
microwave measurements, that circuit 
concepts should predominate. There is a 
dearth of information on waveguides and 
microwave filters using waveguide ele- 
ments. However, the book is thorough in 
ifs treatment and contains a great deal of 
material of value to those who are work- 
ing in the microwave field. ) 


Automatic Feedback Control System 
Synthesis, by Prof. John G. Truxall, 
McGraw-Hill, New York, 1955, 675pp. 

Reviewed by Scorr H. CAMERON 
Armour Research Foundation 

In spite of the fact that the years since 
World War II have seen the development 
of many new and powerful techniques for 
the synthesis of automatic feedback 
systems, much of the design of such 
systems is still carried out in our develop- 
ment laboratories by means of inefficient 
cut and try procedures. This book by 
Dr. Truxall should do much to correct 
this situation. While intended as a text 
for course work at the graduate level, the 
clarity of style and the well chosen ex- 
amples render the book suitable for self 
study. 

The book begins with a discussion of the 
LaPlace transform which is unusually good 
for a servomechanism text. The emphasis 
in this chapter is placed on the time do- 
main-frequency domain relationships which 
are so essential to a firm grasp of the 
material to follow. The LaPlace trans- 
form discussion is followed by a chap- 
ter on signal flow diagrams as applied to 
feedback theory in which the signal flow 
diagram is used to derive and prove some 
of the more important Bode feedback 
theorems with astonishing simplicity. 
Other chapters treat the synthesis of RC 
het works, the synthesis of complete feed- 
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NOW IN 


VOLUME 
... the TC-106 


(USAF TYPE MA-IA) 


PORTABLE 
COMPRESSOR 


for 
GROUND 
SUPPORT 
APPLICATION 


d Typical of the finest results of Continental development is 
the TC-106 portable compressor unit for ground support 
applications. This advanced new model, with a high per- 
formance turbine compressor as its heart, weighs one-third 
less than its predecessor, yet has 17 per cent higher output, 
and in addition, other important qualities: greater mobility, 
less noise, and a completely automatic-'control system. 
. . « It is now in volume production at*the Continental 
Aviation and Engineering Toledo plant. 


ODUCTION 
AT C.A.E. 


MODEL 141 
TURBO-COMPRESSOR 


ENGINE 


C.A.E. gas turbine models—the J69-T-9, the J69-T-2, and the J69- 
T-19A are being built for Cessna’s T-37A twin jet trainer, Temco's 
TT-1 Navy jet trainer, the Beech jet Mentor trainer, and the Ryan 

Q-2A Fire Bee target drone. 


12700 KERCHEVAL AVENUE, DETROIT 15, MICHIGAN 


SUBSIDIARY 


OF CONTINENTAL MOTORS CORPORATION 


631 
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i 
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When you need an aluminum fuel 
pump this big for your rocket to the 
moon, Morris Bean & Company will 
cast it for you. 
Actually, no such casting exists. 
It could. 


We have developed unique foundry 
techniques which we now use in 
production on castings weighing 
thousands of pounds. These castings 
have extensive areas 9 to 12 inches 
thick which must be porosity free. 

We are also casting solid and hol- 
low mandrels used to form solid pro- 
pellant. 

The outstanding contribution of our 
large aluminum castings is drastically 
reduced first and ultimate costs. 

Telephone or write for information 
now. 
Morris Bean & Company 
Yellow Springs 6, Ohio 


| 


back control systems through zero-pole 
configurations, the root-locus method of 
Evans, the application of statistical design 
principles together with the appropriate 
background material, sampled data sys. 
tems, and the analysis of nonlinear systems 
by means of the describing function and 
phase plane techniques. Much of the 
material in these later chapters appears 
here for the first time in book form, hiaving 
appeared previously in the form of tech- 
nical papers and unpublished Dvctora] 
theses. 

The book is entirely theoretical, and 
Dr. Truxall seems to be capable of ex- 
tremely lucid theoretical discussions. H¢ 
first presents the broad theoretical icture 
and then fills in the details by me:ns of 
appropriate examples and well ex cuted 
drawings. For those with a desire for 
still more detail the references are |ibera 
and complete. This reviewer regards 
Dr. Truxall’s book as essential for anyon 
concerned with the design or analysis of 
automatic feedback systems. The singk 
criticism might be the omission o/ an 
problems which, it is noted, might have 
made this text, already long, somewhat 
unwieldy. 


| 

ay Book Notices 

Elements of Engineering Thermody- 
namics, by Rolf H. Sabersky, Mc(Graw- 
Hill, 1957, 318 pp. $9.50. This book has 
been prepared for the undergraduate 
course in engineering thermodynamics. 
The author has succeeded in writing a 
concise and up-to-date book. The many 
illustrations related to jet propulsion 
should make it particularly suitable for 
students interested in this field. 


Engineering Thermodynamics, by ©. 0. 
Mackey, W. N. Barnard and F. O. Ellen- 
wood, John Wiley & Sons, 1957, 428 pp. 
$7.95. Although based on the well- 
known text, “Heat Power Engineering,” 
by Barnard, Ellenwood and Hirshfeld, this 
book actually is an entirely new book. It 
is aimed at undergraduate engineering 
students. 


Books Received 


Strengthening Science Education for 
Youth and Industry, by T. A. Edison Foun- 
dation Institute, New York University 
Press, New York, 1956, 162 pp. $5. This 
constitutes the Proceedings of the Seventh 
Thomas Alva Edison Foundation Insti- 
tute, held Nov. 19-20, 1956. 


Fatigue in Aircraft Structures, edited 
by A. M. Freudenthal, Academic Press, 
1956, 456 pp. $12. This constitutes the 
Proceedings of the International Confer- 
ence on Fatigue in Flight Structures which 
convened at Columbia University in 1956. 
It includes the contributions of the fol- 
lowing authors: R. J. Atkinson, P. D. 
Brooks, J. Cornillon, S. Eggwertz, P. J. 
E. Forsyth, A. M. Freudenthal, E. Gass- 
ner, H. Giddings, R. F. Hanstock, R. A. 
Heller, M. Hempel, H. T. Jensen, W. W. 
Johnstone, P. Kuhn, B. Lundberg, A. 0. 
Payne, C. E. Phillips, R. L. Schleicher, N. 
Thompson, F. Turner, W. Weibull, W. A. 
Wood. 
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These welder phase-shift heat-control patterns were 
directly recorded with a Honeywell 906 Visicorder at _ 
Bristol Aircraft (Western) Limited in Winnipeg. . 

Since the welding heat generated is proportional to 
the square of the current value, phase shift must be accu- 
rately controlled in order to determine the heat value. 
If the phase shift dial is not accurately calibrated, the 
result is too much or too little heat, and a poor weld. 

In this application, the Visicorder is an essential 
guide to accurate calibration, since ink-type recorders 
do not cover the sensitivities and frequencies needed and 
an oscilloscope would present a continually changing 
pattern since most recording periods are less than 10 
cycles. The directly-recorded Visicorder patterns allow 
a convenient study of the exact time when the current 
wave form was being cut off. 

Here is the circuit used in this test. 

Input Circuit for Oscillogram of D. C. Welding ——y 


Lower Arm R 
of Welder series 


To V450-55B 
Galvo & shunt 
No amplification!) 


Welder Arm 
Calibrated 2.16 x 106 
Ohms between 
Connections 


Reference Data: Write for Visicorder Bulletin 
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Honeywell 


O 
across Ton TUBES, 


Model 906A Visicorder 
pictured with 

Record Takeup and — 
Latensifier Unit. 


The Honeywell Visicorder is the first high-fre- 
quency, high-sensitivity direct recording oscillo- 
graph. In laboratories and in the field everywhere, 
instantly-readable Visicorder records are pointing 
the way to new advances in product design, 
rocketry, computing, control, nucleonics ...in any 
field where high speed variables are under study. 

The new Model 906A Visicorder, now avail- 
able in 8- and 14-channel models, produces 
longitudinal grid lines simultaneously with the 
dynamic traces, time lines, and trace identification 
by means of new Accessory units. 

To record high frequency variables—and moni- 
tor them as they are recorded—use the Visicorder 
Oscillograph. Call your nearest Minneapolis- 
Honeywell Industrial Sales Office for a demon- 


Minneapolis Honeywell Regulator Co., Industrial Products Group, Heiland Division, 5200 E. Evans Avenue, Denver 22, Colorado 
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Technical Literature Digest___ 


Fluid Flow 


Approximate Methods for Predicting 
Separation Properties of Laminar Bound- 
ary Layers, by N. Curle and S. W. Skan, 
Aeron. Quart., Aug. 1957, pp. 257-268. 


_ M.H. Smith, Associate Editor, and M. H. Fisher, Contributor 
ca The James Forrestal Research Center, Princeton University 


Jr. and J. W. Westwater, Ind. & Engng 
Chem., vol. 49, Sept. 1957, pp. 1445 
1448. 

On the Steady-state Temperature of a 
Flat Plate Placed in a Flow of Reacting 
Gas Mixture (Russian—Ob_ ustanoviy- 
sheisia temperature ploskoi plastiny, ob- 


‘ tekaemoi reagiruiushchei gasovoi smes’ju), 
On the Calculation of Heat and Mass py E. P. Vaulin, Doklady Akad. Nauk 
Transfer in Laminar and Turbulent SSR, vol. 113, Apr. 21, 1957, pp. 1235- 
[Ds Layers, by D. R. Davies and 1238. 
D. E. Bourne, Quart. J. Mech. Appl. 
Math., vol. 9, Dec. 1956, pp. 457-488. 

A Low-speed Experimental Investiga- 


CONTROLS 


Heat Transfer to Flowing Gas-solids 
Mixtures in a Circular Tube, by Leonard 
Farber and Morgan J. Morley, Jnd. & 


AIR VALVES tion of the Effect of a Sandpaper Type of Engng. Chem., vol. 49, July 1957, pp. 
Boundary-layer Transi- 1143-1150. 
tion, by Albert EF. von Doenhoff and -di i P 
Elmer A. Horton, NACA TN 3858, Oct. Pai 
1956, 45 pp. Univ. Maryland, Inst. Fluid Dynamies 

ACTUATORS Investigation of Transient Pool Boiling Appl. Math., TN BN-86 (AFOSR TN 
Due to Sudden Large Power Surge, by 56-537), (AD 110356), Nov. 195, 14 
Robert Cole, NACA TN 3885, Dec. pp. 

GROUND TEST 1956, 44 pp. Heat Transfer in Fully-developed Flow 

EQUIPMENT Two-phase, Steam-water Critical Flow, between Parallel Plates with Variable 
by H. S. Isbin, J. E. Moy and A. J. R. Heat Sources, by Albert L. Loeffler Jr., 

Write for complete data | DaCruz, AIChE J., vol. 3, Sept. 1957, p Nucl. Sci. & Engng., vol. 2, Sept. 1957. 

or consult the Barber- 61 T Air T 

Colman engineering sales : ransient Air Temperatures in a Duct, 

office nearest you: Heat-Transfer to Boiling Methanol— by S. E. Rea and C. M. Ablow, Trans. 


Los Angeles, Seattle, 

Fort Worth, New York, 
Baltimore, Montreal, 

/ Rockford. 


ASME, vol. 79, no. 7, Oct. 1957, p. 1536. 
Effect of Wall Roughness on Convec- 
tive Heat Transfer in Commercial Pipes, 


Effect of Added Agents, by A. J. Lowery 


Epitor’s Note: Contributions from E. R. 


G. Eckert, J. P. Hartnett, T. F. Irvine Jr. by J. W.- Smith and Norman Epstein, 
BARBER-COLMAN COMPANY and P. J. Schneider of the Heat Transfer Alc hE J., vol. 3, June 1957, pp. 242- 
Dept. | — 1470 Rock Street, Rockford, Illinois. | Laboratory, University of Minnesota, are 248. 

gratefully acknowledged. New 


Finite-difference Technique for 


_ The General Electric Jet Engine Dept. Announces — 
the establishment of 2 new 


Senior Engineering Design Analysis Positions __ 

Because of “expanding activity we are broaden- __ problems relating to engines in operation, and v 
ing our organization to include two additional have provided valid engineering solutions. 
Senior Engineers, and we are inviting inquiries 
from qualified men who wish to join our organiza- 
tion at this level. 


cre 


al 


Briefly, you must be able to accept any hard- 
ware problem that turns up in the field and stay 
with it until you get the answer. 

In keeping with General Electric’s policy of ; 
parallel opportunities for specialists and man- 
agers, salaries established for these new positions 
are equivalent to those of high management level. 


The Engineers who fill these positions will be 
making important contributions to jet engine 
design and utilization. You will be responsible 

We are seeking design specialists with 12 to 20 for both overall and detailed analysis of com- 
years of practical hardware experience. To ponents and systems of the newest designs and 
qualify, you must be thoroughly familiar with all developments to assure maximum utilization of 
mechanical, aerodynamic, hydraulic, control, and special knowledge gained during factory, spe 
installation areas involved in aireraft gas turbine 224. operations. As a member of 
design, with heavy specialized practical experience Senior Review Panel, you will incorporate results 
in at least one of those areas. You will have been ° Your practical experience in our new designs. 
responsible for successfully putting engines And, as a senior specialist in the field, you will 
through their 150 hour test and you will have had represent the Jet Engine Department oe 
many years of experience in analyzing design Government and Industry technical panels. 


GENERAL ELECTRIC 


JET ENGINE DEPT. CINCINNATI 15, OHIO 


If you qualify, call or write 

M. L. Umble, Supervisor 

Professional Placement, POplar 1-4100, 
Extension 3828, Bldg. 501. 
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The Avco Research Laboratory was 
founded a little more than three years © 

ago for the purpose of examining high- — 

temperature gas problems associated with © 

ICBM re-entry. The success of this | 
research led to the birth of a new cor- — 

porate enterprise, Avco’s Research and b 

Advanced Development Division. 


The Research Laboratory, now estab- — 
lished as a separate Avco division, has — 
expanded to embrace all aspects of physi- © 
cal gas dynamics. We are currently gravid — 
with several embryonic projects which we | 
anticipate will likewise grow into new © 
corporate enterprises. Our work in the > 
physics, aerodynamics and chemistry of 
high-temperature gases is growing in the 
following areas: 


Magnetohydrodynamics— 
Flight and industrial power- 
generation applications 


Space flight— 4 
Manned satellites 
Electromagnetic propulsion 


NO MAGNETIC FIELD. This shock tube photograph, taken by 
emitted light only, shows the typical shock wave configuration 


formed by high-velocity gas flowing around a pointed cone. — 


These developments have created a num- 
ber of openings for physicists, aerody- 
namicists and physical chemists. If your 
background qualifies you to work in any 
of these areas, we would be pleased to 


hear from you. 


Dr. Arthur Kanirowitz, Director 
Avco Research Laboratory 


WITH MAGNETIC FIELD. Here is shown the magnetohydrodynamic 
displacement of the shock wave. The magnetic field is caused 
by electric current flowing through a coil of wire within the cone. 
This experiment qualitatively demonstrates the interaction of 
a high-temperature gas with a magnetic field. This effect would 
be expected to produce drag and reduce heat transfer to the body. 


P, §, A listing of laboratory research re- 
ports indicative of the scope and depth 
of our activities is available. Address 
your request: Attention: Librarian, Avco 
Research Laboratory, 2385 Revere Beach 
Parkway, Everett, Massachusetts. 


ies of electrically conducting fluids interacting 
RESEARCH with magnetic fields. 
LABORATORY 


Crosley Broadcasting Corporation 
Research and Advanced Development Division 
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A Division of Avco Manufacturing Corporation/Everett, Mass. 


Other divisions and subsidiaries are: 


Ezee Flow Division 
Lycoming Division 


AK Division 
Crosley Division 


New Idea Division © 


= 


Moffats Limited 
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OGRESS REPORT FROM AVCO ESEARCH LABORATO ae. 
PR R RY 
: 


A 


broad creative mission at Chance Vought 


MILITARY ASTRONAUTICS: Vought studies toward space research 
vehicles and manned spacecraft include multistaging, space 
communications, nuclear and ionic propulsion, celestial navi- 
gation, and membership on Boeing’s Dyna Soar space glider 
development team. 


ASW RESPONSIBILITY: Under the Office of Naval Research, 
Vought antisubmarine specialists are conducting extensive 
research and development work in ASW detection and clas- 
sification phases. Goal is to bring detection abreast of 
destruction capabilities. 


Advanced Weapons Staff Engineer. Ph.D. preferred, with at 
least 10 years background in guidance or navigation and control 
systems. To develop completely new concepts in guidance, naviga- 
tion, or control systems. 

Electro-mechanical Systems Engineer or Specialist. A.E., E.E., 
or M.E., (advanced degree preferred) with at least 7 years expe- 
rience in autopilot, flight control, stability systems and inertial guid- 
ance systems and design work. To make high-level technical studies 
of various control and stabilization systems for advanced weapons. 
Guidance Design Engineer. E.E. or Physics Degree, plus 2 or 
more years experience. To design various active and self-contained 
missile guidance systems, and to design and develop radar beacons. 
Lead Instrumentation Engineer. B.S. or M.S. in E.E., plus 5 
years electronics experience in telemeter applications of flight test 
instrumentation data. To guide a Research and Development group 
in the design and application of electronic gathering and reduction 
systems. 

Hydrodynamicist. B.S. or M.S. in Engineering plus 5-7 years 
experience in hydrodynamics and cavitation of torpedoes and other 
fully submerged vehicles. 

Qualified engineers and scientists who would like to join Vought’s 
development of advanced weapons are invited to — 7 


ml A. L. Jarrett, Manager, Advanced 
ees Weapons Engineering, Dept. JP-3. 
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SEES 


Solution of the Heat-conduction Equa 
tions, Especially Near Surfaces with Con. 
vective Heat Transfer, by H. G. Elrod 
Jr., Trans. ASME, vol. 79, Oct. 1957, p, 
1519. 

On the Use of Matrices for Solvi 
Periodic Heat Flow Problems, by G. W. T, 
White, Appl. Sci. Res., vol. 6, Sec. & 
Nos. 5-6, 1957, pp. 433-444. 

Influence of Turbulence on Transfer 
Heat from Cylinders, by J. Kestin 
P. F. Meader, NACA TN 4018 
1957, 78 pp. 

Measurements of Heat Transfer from 
Fine Wires in Supersonic Flows, by John 
Laufer and Robert McClellan, J. Fluid 
Mech., vol. 3, Sept. 1956, pp. 276-289. 

Heat Transfer and Flow-Friction Char. 
acteristics of Woven-screen and Crvssed- 
rod Matrixes, by L. 8S. Tong and A. L, 
London, Trans. ASME, vol. 79, Oct, 
1957, p. 1558. 

Heat Transfer to Lead-Bismuth in 
Turbulent Flow in an Annulus, by R. A. 
Seban and D. F. Casey, Trans. ASME, 
vol. 79, Oct. 1957, p. 1514. 

Experiments on the Inhibition of Ther- 
mal Convection by a Magnetic Field, by 
Y. Nakagawa, Proc. Roy. Soc., Ser. A, 
no. 1220, vol. 2401, 1957, pp. 108-1 

Cross Sections of Dissociative Re- 
combination, by Ernest Bauer and Ta. 
You Wu, Can. J. Phys., vol. 34, Dee, 
1956, pp. 1436-1447. 


Combustion, Fuels and 
Propellants 


Final Technical Report, Contract AF 
18(600)-430, Texas Univ. (AFOSR-TN- 
56-63; ASTIA AD 115028), Dec. 1956, 
18 pp. 

Thermodynamic Calculations of the 
Performance of the Hydrogen Fluorine 
System, by Rudolph Edse and William L. 
Doyle, Wright Air Dev. Center, Tech. Rep. 
53-427, March 1956, 46 pp. 

Theoretical Combustion Performance of 


| Ramjet Fuels: Hydrogen, by Waldo T. 


Renich, Johns Hopkins Univ. Appl. 
Phys. Lab., CF 2601, Dec. 1956, 3 pp. 5 fig. 

Research on Some Properties of Certain 
Propellants, by Jacques Francon, France, 
Ministere de Air, Publications Scien- 
tifiques et Techniques, Notes Techniques 66, 
1957, 46 pp. (in French). 

The Diesel Powered Ballistic Centrifuge 
by W. Hendricks and W. M. Bogart, 
Hercules Powder Co., Allegany Ball. Lab., 
Rep. ABL/X-17, Aug. 1957, 13 pp. 

Flame Stabilizing Effects of Inclined 
Air Jets, by D. P. Duclos, Allan Schaffer 
and A. B. Cambel, Indust. & Engng. 
Chem., vol. 49, Dec. 1957, pp. 2063-2066. 

Investigations of the Detonation Prop- 
erties of Condensed Explosives with 
Equations of State Based on Intermolec- 
ular Potentials, I: RDX with Fixed 
Product Composition, by W. Fickett, 
W. W. Wood and Z. W. Salsburg, J. 
Chem. Phys., vol. 27, Dec. 1957, pp. 1324 
1329. 

The Thickness of Laminar Flames im 
Pre-mixed Reactants: Optical Considera- 
tions, by F. J. Weinberg, Proc. Royal Soc. 
vol. 243, Dec. 10, 1957, pp. 107-118. 

Progress Report on Study of Combus- 
tion for Supersonic Ram Jet (ODS Pro- 
ject) for Period Apr. 1—Sept. 30, 1945, CM 
174, (Report PDN-2980), by Esso Labs., 
Oct. 1945, 24 pp. 3 tab. (Declassified 
from Confidential by authority of Johns 
Hopkins Univ., Appl. Phys. Lab., letter 
dated 2/5/57, signed Paul E. Clark.) 
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